[gradsusr] WAFS GRIB2: Cannot contour grid - all undefined values
matt Grzych
mattgrzych at yahoo.com
Sun Dec 15 14:37:14 EST 2013
Ryan,
Thanks for the advice about the -verf flag. I did set it up and tested with -verf, but unfortunately I'm still getting the "all undefined values" warning.
thanks,
Matt
Sent from my iPad
On Dec 14, 2013, at 7:47 PM, "Maue, Ryan" <rnm02c at my.fsu.edu> wrote:
> g2ctl -verf
>
> gribmap -i
>
>
> f06 = 6-hr forecast on its own would require the -verf flag
>
>
> Cheers,
>
> RYAN
>
>
> From: gradsusr-bounces at gradsusr.org <gradsusr-bounces at gradsusr.org> on behalf of Jeff Duda <jeffduda319 at gmail.com>
> Sent: Saturday, December 14, 2013 3:57 PM
> To: Matthew Grzych; GrADS Users Forum
> Subject: Re: [gradsusr] WAFS GRIB2: Cannot contour grid - all undefined values
>
> Hey, Matt,
> Great to hear from you!
>
> I don't have a ton of experience working with GRIB2 files, but I do see in the description of the CAT*, CTP*, and ICIP* variables different pressure levels listed. Is it possible that these arrays are valid on different sets of pressure fields that are not covered in your ZDEF entry? You only have 40000 35000 30000 25000 20000 15000 listed, but I wonder if it will work if you add 80000 70000 60000 and 50000 to the ZDEF entry.
>
> Jeff
>
>
> On Sat, Dec 14, 2013 at 4:43 PM, Matthew Grzych <mattgrzych at yahoo.com> wrote:
>> I am having no luck getting WAFS GRIB2 files to plot in GrADS and get the "Cannot contour grid - all undefined values" error after many hours of trying seemingly everything. I am downloading the files from:
>>
>> http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.2013121212/WAFS_blended_2013121212f06.grib2
>>
>> Then I generate the control files using g2ctl.pl and index files using gribmap. Here's an example ctl file:
>> _______________________________________________________________________________________
>> dset /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2
>> index /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2.idx
>> undef -0.5
>> title /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2
>> * produced by g2ctl v0.0.8.4
>> * command line options: /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2
>> * griddef=1:0:(288 x 145):grid_template=0:winds(N/S): lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.250000 lon 0.000000 to 358.750000 by 1.250000 #points=41760:winds(N/S)
>> dtype grib2
>> ydef 145 linear -90.000000 1.25
>> xdef 288 linear 0.000000 1.250000
>> tdef 1 linear 18Z14dec2013 6hr
>> * PROFILE hPa
>> zdef 6 levels 40000 35000 30000 25000 20000 15000
>> options pascals
>> vars 9
>> CATaveprs 6,100 0,19,22,0 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]
>> CATmaxprs 6,100 0,19,22,2 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]
>> CBHEclm 0,10,0 0,6,25 ** entire atmosphere Horizontal Extent of Cumulonimbus (CB) [%]
>> CTPaveprs 6,100 0,19,21,0 ** (700 600 500 400 300) In-Cloud Turbulence [%]
>> CTPmaxprs 6,100 0,19,21,2 ** (700 600 500 400 300) In-Cloud Turbulence [%]
>> ICAHTcbb 0,11,0 0,3,3 ** cumulonimbus base ICAO Standard Atmosphere Reference Height [m]
>> ICAHTcbt 0,12,0 0,3,3 ** cumulonimbus top ICAO Standard Atmosphere Reference Height [m]
>> ICIPaveprs 6,100 0,19,20,0 ** (800 700 600 500 400 300) Icing [%]
>> ICIPmaxprs 6,100 0,19,20,2 ** (800 700 600 500 400 300) Icing [%]
>> ENDVARS
>> ____________________________________________________________________________________
>>
>> Here's the output from my .gs file:
>> ___________________________________________________________________________________
>> Running in Batch mode
>> open WAFS_blended_2013121412f06.grib2.ctl
>> Scanning description file: WAFS_blended_2013121412f06.grib2.ctl
>> Data file /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2 is open as file 1
>> LON set to 0 360
>> LAT set to -90 90
>> LEV set to 400 400
>> Time values set: 2013:12:14:18 2013:12:14:18
>> E set to 1 1
>> set LON 0 358.75
>> LON set to 0 358.75
>> set gxout shaded
>> d CATmaxprs
>> Cannot contour grid - all undefined values
>> q dims
>> Default file number is: 1
>> X is varying Lon = 0 to 358.75 X = 1 to 288
>> Y is varying Lat = -90 to 90 Y = 1 to 145
>> Z is fixed Lev = 400 Z = 1
>> T is fixed Time = 18Z14DEC2013 T = 1
>> E is fixed Ens = 1 E = 1
>> quit
>> EXEC error: error in turb_dev.gs. EXEC stopped.
>> No hardcopy metafile open
>> GX package terminated
>>
>> __________________________________________________________________________________
>>
>> Any ideas?
>> Thanks,
>> Matt
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> --
> Jeff Duda
> Graduate research assistant
> University of Oklahoma School of Meteorology
> Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131215/325bb071/attachment-0001.html
More information about the gradsusr
mailing list