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

Jennifer Adams jma at cola.iges.org
Sun Dec 22 16:58:29 EST 2013


Matt,
You have ZDEF 10, but "6,100" in the levels field of your variable  
declaration. GrADS assumes that the 6 levels for those variables are  
the first six in the ZDEF list, which is why you're not getting any  
plots for 300mb and above. Change your "6,100" to "10,100" and then  
you will see the grids for the variables that are defined at higher  
levels.
--Jennifer

On Dec 22, 2013, at 1:16 PM, Matthew Grzych wrote:

> Jennifer,
>
> This fix worked great for me when testing lev 400, but I have  
> scripted to make images for all levels and only levels 350 and 400  
> work. For levels 300, 250, and 200 I am still getting this error:
>
> "Data Request Warning:  Request is completely outside file limits
>   Entire grid contents are set to missing data
>   Grid limits of file:     X = 1 288  Y = 1 145  Z = 1 10  T = 1 1   
> E = 1 1
>   Grid limits of request:  X = 1 288  Y = 1 145  Z = 9 9  T = 1 1  E  
> = 1 1
>   Warning issued for variable = catmaxprs
> Cannot contour grid - all undefined values"
>
> it says im out of the range, but it clearly shows I am not. Any  
> ideas here?
>
> My control file:
> dset /home/matt/gfs_wafs/12z/WAFS_f18.grib2
> index /home/matt/gfs_wafs/12z/WAFS_f18.grib2.idx
> undef 9.999E+20
> title /home/matt/gfs_wafs/12z/WAFS_f18.grib2
> * produced by g2ctl v0.0.8.4
> * command line options: -verf /home/matt/gfs_wafs/12z/WAFS_f18.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 06Z22dec2013 1mo
> * PROFILE hPa
> zdef 10 levels 80000 70000 60000 50000 40000 35000 30000 25000 20000  
> 15000
> options pascals
> vars 9
> CATaveprs   6,100  0,19,22,0,3 ** (400 350 300 250 200 150) Clear  
> Air Turbulence (CAT) [%]
> CATmaxprs   6,100  0,19,22,2,3 ** (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,3 ** (700 600 500 400 300) In-Cloud  
> Turbulence [%]
> CTPmaxprs   6,100  0,19,21,2,3 ** (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,3 ** (800 700 600 500 400 300) Icing [%]
> ICIPmaxprs   6,100  0,19,20,2,3 ** (800 700 600 500 400 300) Icing [%]
> ENDVARS
>
>
> My script:
> open /home/matt/gfs_wafs/12z/WAFS_f18.grib2.ctl
>
> *set mpdraw off
> *set grid off
> set lon 0 358.75
> set lev 200
>
> set cmin 5
> set cmax 50
> set rgb 21 255 51 255
> set rgb 22 255 0 25
> set rgb 23 204 0 204
> set rgb 24 153 0 153
> set rgb 25 102 0 102
> set rgb 26 204 0 102
> set clevs 4 8 12 16 20 24
> set ccols 0 21 23 24 25 26 22
> set gxout shade2b
> d CATmaxprs
> printim /home/matt/gfs_wafs/images/turb_f18_200.png x3200 y2400
> quit
>
>
> Thanks!
> Matt
>
>
>
>
> On Monday, December 16, 2013 10:59 AM, Matthew Grzych <mattgrzych at yahoo.com 
> > wrote:
> Jennifer,
>
> That works for me! Thank you so much for the help as I have been  
> banging my head over this for days.
>
> I added the additional zdef levels (to make a total of 10) and then  
> added the additional "3" to the grib codes and it works great now.
>
> Thanks again,
> Matt
>
>
>
> On Monday, December 16, 2013 7: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
>
>
>
>
>
>
>
>
>

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131222/f1fdceb0/attachment-0001.html 


More information about the gradsusr mailing list