[gradsusr] WAFS GRIB2: Cannot contour grid - all undefined values

[matt Grzych]

Hey Jeff,

Good to hear from you too!

Good point on the pressure levels. I changed the zdef line in my ctl file to:

zdef 10 levels 80000 70000 60000 50000 40000 35000 30000 25000 20000 15000

Thought for sure this would fix it, but nope! still getting "Cannot contour grid - all undefined values". Ugh. 

I also changed g2ctl.pl to have "-verf" as suggested by Ryan, which didn't fix this issue.  What else is there? 

Thanks,
Matt


Sent from my iPad

On Dec 14, 2013, at 4:57 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:

> 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/11e02d69/attachment.html