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

[Jennifer Adams]

Hi, Matt -- 
A few changes need to be made to your descriptor file. There are more than 6 pressure levels, so change your ZDEF statement to:

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

Not all variables are defined at all these levels, but this set of levels is the superset to catch all the records in the file. 

Also, you need to add one more grib code to the units field in your variable declarations:

CBHEclm     0,10,0  0,6,25      ** entire atmosphere Horizontal Extent of Cumulonimbus (CB) [%]9,
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]
CATaveprs   10,100  0,19,22,0,3 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]
CATmaxprs   10,100  0,19,22,2,3 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]
CTPaveprs   10,100  0,19,21,0,3 ** (700 600 500 400 300) In-Cloud Turbulence [%]
CTPmaxprs   10,100  0,19,21,2,3 ** (700 600 500 400 300) In-Cloud Turbulence [%]
ICIPaveprs  10,100  0,19,20,0,3 ** (800 700 600 500 400 300) Icing [%]
ICIPmaxprs  10,100  0,19,20,2,3 ** (800 700 600 500 400 300) Icing [%]

These changes resulted in a "MATCH" for all records in your file, using gribmap (in verbose mode) in GrADS 2.0.2. 
--Jennifer




On Dec 14, 2013, at 5:43 PM, Matthew Grzych 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

--
Jennifer M. Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
111 Research Hall, Mail Stop 2B3
George Mason University
4400 University Drive
Fairfax, VA 22030 





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131216/6ad6717a/attachment-0001.html