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

Mon Dec 16 11:06:31 EST 2013

Jennifer,

  Next version of wgrib2 will add the missing "3" to the ctl file.

   Wesley

On Mon, Dec 16, 2013 at 9:16 AM, Jennifer Adams <jma at cola.iges.org> wrote:

> Hi, Matt --
> A few changes need to be made to your descriptor file. There are more than
> 6 pressure levels, so change your ZDEF statement to:
>
> zdef 10 levels 80000 70000 60000 50000 40000 35000 30000 25000 20000 15000
>
> Not all variables are defined at all these levels, but this set of levels
> is the superset to catch all the records in the file.
>
> Also, you need to add one more grib code to the units field in your
> variable declarations:
>
> CBHEclm     0,10,0  0,6,25      ** entire atmosphere Horizontal Extent of
> Cumulonimbus (CB) [%]9,
> 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]
> CATaveprs   10,100  0,19,22,0,3 ** (400 350 300 250 200 150) Clear Air
> Turbulence (CAT) [%]
> CATmaxprs   10,100  0,19,22,2,3 ** (400 350 300 250 200 150) Clear Air
> Turbulence (CAT) [%]
> CTPaveprs   10,100  0,19,21,0,3 ** (700 600 500 400 300) In-Cloud
> Turbulence [%]
> CTPmaxprs   10,100  0,19,21,2,3 ** (700 600 500 400 300) In-Cloud
> Turbulence [%]
> ICIPaveprs  10,100  0,19,20,0,3 ** (800 700 600 500 400 300) Icing [%]
> ICIPmaxprs  10,100  0,19,20,2,3 ** (800 700 600 500 400 300) Icing [%]
>
> These changes resulted in a "MATCH" for all records in your file, using
> gribmap (in verbose mode) in GrADS 2.0.2.
> --Jennifer
>
>
>
>
> On Dec 14, 2013, at 5:43 PM, Matthew Grzych 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
>
>
> --
> Jennifer M. Adams
> Center for Ocean-Land-Atmosphere Studies (COLA)
> 111 Research Hall, Mail Stop 2B3
> George Mason University
> 4400 University Drive
> Fairfax, VA 22030
>
>
>
>
>
>
> _______________________________________________
> 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/b7e64a44/attachment.html 