<div dir="ltr"><div><div>Jennifer,<br><br></div> Next version of wgrib2 will add the missing "3" to the ctl file.<br><br></div> Wesley<br><br><div><div><div><div class="gmail_extra"><div class="gmail_quote">On Mon, Dec 16, 2013 at 9:16 AM, Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hi, Matt -- </div><div>A few changes need to be made to your descriptor file. There are more than 6 pressure levels, so change your ZDEF statement to:</div>
<div><br></div><div><div><font face="Courier"><span style="font-size:13px">zdef 10 levels 80000 70000 60000 50000 40000 35000 30000 25000 20000 15000</span></font></div></div><div><font face="Courier"><span style="font-size:13px"><br>
</span></font></div><div>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. </div><div><br></div><div>Also, you need to add one more grib code to the units field in your variable declarations:</div>
<div><br></div><div><div><font face="Courier"><span style="font-size:13px">CBHEclm 0,10,0 0,6,25 ** entire atmosphere Horizontal Extent of Cumulonimbus (CB) [%]9,</span></font></div><div><font face="Courier"><span style="font-size:13px">ICAHTcbb 0,11,0 0,3,3 ** cumulonimbus base ICAO Standard Atmosphere Reference Height [m]</span></font></div>
<div><font face="Courier"><span style="font-size:13px">ICAHTcbt 0,12,0 0,3,3 ** cumulonimbus top ICAO Standard Atmosphere Reference Height [m]</span></font></div><div><font face="Courier"><span style="font-size:13px">CATaveprs 10,100 0,19,22,0,3 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]</span></font></div>
<div><font face="Courier"><span style="font-size:13px">CATmaxprs 10,100 0,19,22,2,3 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]</span></font></div><div><font face="Courier"><span style="font-size:13px">CTPaveprs 10,100 0,19,21,0,3 ** (700 600 500 400 300) In-Cloud Turbulence [%]</span></font></div>
<div><font face="Courier"><span style="font-size:13px">CTPmaxprs 10,100 0,19,21,2,3 ** (700 600 500 400 300) In-Cloud Turbulence [%]</span></font></div><div><font face="Courier"><span style="font-size:13px">ICIPaveprs 10,100 0,19,20,0,3 ** (800 700 600 500 400 300) Icing [%]</span></font></div>
<div><font face="Courier"><span style="font-size:13px">ICIPmaxprs 10,100 0,19,20,2,3 ** (800 700 600 500 400 300) Icing [%]</span></font></div></div><div><br></div><div>These changes resulted in a "MATCH" for all records in your file, using gribmap (in verbose mode) in GrADS 2.0.2. </div>
<div>--Jennifer</div><div><br></div><div><br></div><div><br></div><br><div><div>On Dec 14, 2013, at 5:43 PM, Matthew Grzych wrote:</div><br><blockquote type="cite"><div><div style="font-size:10pt;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
<div>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:<br>
</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif"><br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
<a href="http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.2013121212/WAFS_blended_2013121212f06.grib2" target="_blank">http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.2013121212/WAFS_blended_2013121212f06.grib2</a></div>
<div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif"><br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
Then I generate the control files using <a href="http://g2ctl.pl" target="_blank">g2ctl.pl</a> and index files using gribmap. Here's an example ctl file:</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
_______________________________________________________________________________________<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
dset
/home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2<br>index /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2.idx<br>undef -0.5<br>title /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2<br>* produced by g2ctl v0.0.8.4<br>
* command line options: /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2<br>* 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)<br>
dtype grib2<br>ydef 145 linear -90.000000 1.25<br>xdef 288 linear 0.000000 1.250000<br>tdef 1 linear 18Z14dec2013 6hr<br>* PROFILE hPa<br>zdef 6 levels 40000 35000 30000 25000 20000 15000<br>options pascals<br>vars 9<br>CATaveprs 6,100 0,19,22,0 ** (400 350 300 250 200 150) Clear Air Turbulence (CAT) [%]<br>
CATmaxprs 6,100 0,19,22,2 ** (400 350 300 250 200 150) Clear Air
Turbulence (CAT) [%]<br>CBHEclm 0,10,0 0,6,25 ** entire atmosphere Horizontal Extent of Cumulonimbus (CB) [%]<br>CTPaveprs 6,100 0,19,21,0 ** (700 600 500 400 300) In-Cloud Turbulence [%]<br>CTPmaxprs 6,100 0,19,21,2 ** (700 600 500 400 300) In-Cloud Turbulence [%]<br>
ICAHTcbb 0,11,0 0,3,3 ** cumulonimbus base ICAO Standard Atmosphere Reference Height [m]<br>ICAHTcbt 0,12,0 0,3,3 ** cumulonimbus top ICAO Standard Atmosphere Reference Height [m]<br>ICIPaveprs 6,100 0,19,20,0 ** (800 700 600 500 400 300) Icing [%]<br>
ICIPmaxprs 6,100 0,19,20,2 ** (800 700 600 500 400 300) Icing [%]<br>ENDVARS</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
____________________________________________________________________________________</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">Here's the output from my .gs file:</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
___________________________________________________________________________________<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
Running in Batch mode<br>open WAFS_blended_2013121412f06.grib2.ctl<br>Scanning description file: WAFS_blended_2013121412f06.grib2.ctl<br>Data file /home/matt/gfs_wafs/12z/WAFS_blended_2013121412f06.grib2 is open as file 1<br>
LON set to 0 360<br>LAT set to -90 90<br>LEV set to 400 400<br>Time values set: 2013:12:14:18 2013:12:14:18<br>E set to 1 1<br>set LON 0 358.75<br>LON set to 0 358.75<br>set gxout shaded<br>d CATmaxprs<br>Cannot contour grid - all undefined values<br>
q dims<br>Default file number is: 1<br>X is varying Lon = 0 to 358.75 X = 1 to 288<br>Y is varying Lat = -90 to 90 Y = 1 to 145<br>Z is fixed Lev = 400 Z = 1<br>T is fixed Time = 18Z14DEC2013 T = 1<br>E is fixed Ens = 1 E = 1<br>
quit<br>EXEC error: error in <a href="http://turb_dev.gs" target="_blank">turb_dev.gs</a>. EXEC
stopped.<br>No hardcopy metafile open<br>GX package terminated<br><br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
__________________________________________________________________________________</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">Any ideas?</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
Thanks,</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">Matt<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif">
<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif"><br><br></div></div></div>_______________________________________________<br>
gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div><br><div>
<div style="word-wrap:break-word;font-size:12px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><div>
--</div><div>Jennifer M. Adams<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>111 Research Hall, Mail Stop 2B3<br>George Mason University<br>4400 University Drive<br>Fairfax, VA 22030 <br><br></div><div><br></div><br>
</span></span></div><br>
</div>
<br></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div></div></div></div></div>