[gradsusr] WAFS GRIB2: Cannot contour grid - all undefined values
Matthew Grzych
mattgrzych at yahoo.com
Sun Dec 22 13:16:30 EST 2013
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131222/efb8573f/attachment-0001.html
More information about the gradsusr
mailing list