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

Maue, Ryan rnm02c at my.fsu.edu
Sat Dec 14 21:47:42 EST 2013


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<mailto: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<http://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<http://turb_dev.gs>.  EXEC stopped.
No hardcopy metafile open
GX package terminated

__________________________________________________________________________________

Any ideas?
Thanks,
Matt




_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto: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/9f160df8/attachment-0001.html 


More information about the gradsusr mailing list