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

Sun Dec 15 18:16:21 EST 2013

I'm not following your problem real closely...but shot in the dark here....
try:

g2ctl -no_profile -verf

If the problem you're having is with variables at a height, this will
eliminate that so you don't have to set the level, you get a unique
variable for each one... IE instead of set lev 500 and then pull tmp you
would have a new variable called t500mb or something like that.

Be sure to do a 'q file' to see what the vars are (or just open the ctl in
a text editor)

Andrew Revering - Follow
@andrewrevering<https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>on
Twitter
Convective Development, Inc.
PO Box 511, Cambridge, MN 55008
http://www.convectivedevelopment.com/

On Sun, Dec 15, 2013 at 12:39 PM, matt Grzych <mattgrzych at yahoo.com> wrote:

> 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
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131215/f50214ed/attachment.html 