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

Matthew Grzych mattgrzych at yahoo.com
Sat Dec 14 17:43:59 EST 2013


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131214/681b5034/attachment.html 


More information about the gradsusr mailing list