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

Wesley Ebisuzaki - NOAA Federal wesley.ebisuzaki at noaa.gov
Mon Dec 16 08:55:16 EST 2013 

Matt,

  I couldn't get the g2ctl/gribmap combo to work either.  When I did a
"gribmap -v -i junk.ctl",
most of the fields didn't match.  Not sure why.    It may be a version
issue as I am only
using GrADS  v 2.0.1.

   I was able to get the plots to work with the alt_g2ctl/alt_gmb
combination.

  152  alt_g2ctl WAFS_blended_2013121212f06.grib2 >junk.ctl
  153  vi junk.ctl
  154  alt_gmp -i junk.ctl

Line 152 made a ctl file with too long variable names.  I had to change the
variables
names to something shorter.  Changed "dspatial" to "s".

The status of g2ctl and alt_g2ctl is

   g2ctl/gribmap:  stable, handles ensemble dimension and files > 2GB
          Can create ctl/idx files that alt_g2ctl/alt_gmp can't.

   alt_g2ctl/alt_gmp:  late beta, no ensemble dimension, files < 2GB,
           faster for large templated data sets, Can create ctl/idx files
that
           g2ctl/gribmap can't.  g2ctl v1.0 is getting closer to release.

Wesley







On Sun, Dec 15, 2013 at 6:16 PM, Andrew Revering <andy at f5data.com> wrote:

> 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
>>
>>
>
> _______________________________________________
> 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/20131216/e3457cd8/attachment-0001.html

More information about the gradsusr mailing list