Replacing line 115 of dppgrb.f:
IF ( qdif .ne. 0.0 ) THEN
with:
IF ( ABS(qdif) .gt. TINY(qdif) * imax ) THEN
and rebuilding should fix this bug.
Steve
On 6/20/19 7:51 PM, elowery@xxxxxxxxxxxx wrote:
> Hello all,
>
> I have a legacy processes which is still using 2.50 degree GFS data
> (being phased out). NOAA stopped creating the files several years ago,
> and I've been using NAGRIB2 and GDBIINT to convert 1.00, 0.50, or 0.25
> degree grib2 -> 2.50 degree grid data.
>
> Occasionally, GDBIINT will hang at ~98% CPU when converting a file.
> Testing this across both older and newer versions of GEMPAK results in
> the same hung GDBIINT process. Attached is a link to some test files
> which GDBIINT handles flawlessly, and one file which it does not.
>
> *Source:*
> https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.20190616/00/
>
> *Files (GDBIINT status)*
> gfs.t00z.pgrb2.1p00.f159 (ok)
> gfs.t00z.pgrb2.1p00.f162 (hangs)
> gfs.t00z.pgrb2.1p00.f357 (ok)
> gfs.t00z.pgrb2.0p50.f162 (ok)
>
> Google drive link (gfs_sample_files.zip
> <https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdrive.google.com%2Ffile%2Fd%2F1rXC8Y8UaRa7FHzUe9-Hgiv2QJVBhtnUG%2Fview%3Fusp%3Dsharing&data=02%7C01%7Csgdecker%40envsci.rutgers.edu%7C7074cde7102a4190266a08d6f5da39cb%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C1%7C636966714966523229&sdata=sIahh0S62h0lornH62id9qAMkH1RWvjU1PLk%2BqhioEU%3D&reserved=0>)
>
> *Sample script*
> convert_1p00_2p50.sh {res} {hhh}
>
> *Sample command (convert 1p00 159 file to 2p50)*
> convert_1p00_2p50.sh 1p00 159
>
> Has anyone come across this? My assumption is that it's a corrupt grib2
> file from NOAA and should be brought to their attention as well, but
> maybe it could be handled gracefully in GDBIINT?
>
> This error randomly occurs and is not consistent with a specific file or
> resolution.
>
> Thanks in advance,
>
> Evan Lowery
>
>
> _______________________________________________
> NOTE: All exchanges posted to Unidata maintained email lists are
> recorded in the Unidata inquiry tracking system and made publicly
> available through the web. Users who post to any of the lists we
> maintain are reminded to remove any personal information that they
> do not want to be made public.
>
>
> gembud mailing list
> gembud@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe, visit:
> https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.unidata.ucar.edu%2Fmailing_lists%2F&data=02%7C01%7Csgdecker%40envsci.rutgers.edu%7C7074cde7102a4190266a08d6f5da39cb%7Cb92d2b234d35447093ff69aca6632ffe%7C1%7C0%7C636966714966543256&sdata=eRp2O%2FIGjv0rfEDAFMQ4%2FIwcdVUlvsFxN3C3rEni0IQ%3D&reserved=0
>
--
Steve Decker, Associate Teaching Professor
Director, Meteorology Undergraduate Program
Department of Environmental Sciences Phone: (848) 932-5750
Rutgers University Fax: (732) 932-8644
14 College Farm Rd Email: decker@xxxxxxxxxxxxxxxxxx
New Brunswick, NJ 08901-8551
