[gradsusr] WAFS GRIB2: Cannot contour grid - all undefined values
Jeff Duda
jeffduda319 at gmail.com
Sat Dec 14 18:57:19 EST 2013
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/20131214/3b20c57e/attachment.html
More information about the gradsusr
mailing list