<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>I grabbed a sample file: </div><div><a href="ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/RD.20120119/PT.grid_DF.bb/fh.0006_tl.press">ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/RD.20120119/PT.grid_DF.bb/fh.0006_tl.press</a></div><div><br></div><div>There are 3 different variables (grib codes 6, 7, and 102). Each variable has 8 records, each with a different grid, except for variable 102, which has 9. The ninth grid for variable 102 is a latlon grid with 1813 points, all the others are thinned grids with nxny=3447. Grib2ctl creates a pdef file for the thinned grid, presumably based on the grid specs it finds in record #1. I ran wgrib -V on this file and then grepped on 'grid=' and got this:</div><div><br></div><div><div>rec 1:00040:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=44</font> tropopause 6hr fcst:</div><div>rec 2:05915:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=41</font> tropopause 6hr fcst:</div><div>rec 3:11790:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=43</font> tropopause 6hr fcst:</div><div>rec 4:17233:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=37</font> tropopause 6hr fcst:</div><div>rec 5:23108:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=39</font> tropopause 6hr fcst:</div><div>rec 6:28983:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=42</font> tropopause 6hr fcst:</div><div>rec 7:34858:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=38</font> tropopause 6hr fcst:</div><div>rec 8:40733:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) <font class="Apple-style-span" color="#ff221b">grid=40</font> tropopause 6hr fcst:</div></div><div><br></div><div>When checking for matches for PDEF'd grids, gribmap will check the nxny value to make sure it matches what the PDEF entry says, but doesn't go beyond that. The one latlon grid in the file does spawn an error from gribmap, "GRIB grid size does not match descriptor". Each of the 8 thinned grids will appear identical to gribmap, so they will all match, and the one that the index file will point to will be the last record in the file to match, which will not likely correspond to the information used to create the pdef file. Wesley's suggestion to break up the file and separate out the nine different grids is the only way this data can be read with GrADS. </div><div>--Jennifer</div><div><br></div><br><div><div>On Dec 16, 2011, at 5:06 PM, Wesley Ebisuzaki wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Brian,<br><br>Looks like a problem with gribmap in GrADS v2.<br><br>The files have several different grids in them. GrADS should check that <br>the grid<br>matches the grid in the ctl file (it does with grib1). It doesn't so all <br>the grids match<br>and you can plot the wrong data. :(<br><br>Options: grib1, delete all but one of the grids.<br><br><blockquote type="cite">18.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">19.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">20.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">21.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">22.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">23.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">24.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">25.1: lev1=101,0 var=0,3,1 MATCH<br></blockquote><blockquote type="cite">grib2map: Writing out the index file<br></blockquote>Wesley<br><br><br><br>Walsh, Brian J wrote:<br><blockquote type="cite"><br></blockquote><blockquote type="cite">Hi,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I’m wondering if anyone has been able to plot the UKMET grib data <br></blockquote><blockquote type="cite">available on the TGFTP ftp site <br></blockquote><blockquote type="cite">(<a href="ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/">ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/</a>) using <br></blockquote><blockquote type="cite">GrADS. It is a thinned lat-long grid and I am able to convert it to a <br></blockquote><blockquote type="cite">GRIB2 file using cnvgrib and run wgrib2 on the files, but it’s not <br></blockquote><blockquote type="cite">plotting. Just wanted to see if anyone had success plotting these <br></blockquote><blockquote type="cite">files and how it was accomplished.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Brian<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The information contained in this e-mail is intended only for the <br></blockquote><blockquote type="cite">individual or entity to whom it is addressed.<br></blockquote><blockquote type="cite">Its contents (including any attachments) may contain confidential <br></blockquote><blockquote type="cite">and/or privileged information.<br></blockquote><blockquote type="cite">If you are not an intended recipient you must not use, disclose, <br></blockquote><blockquote type="cite">disseminate, copy or print its contents.<br></blockquote><blockquote type="cite">If you receive this e-mail in error, please notify the sender by reply <br></blockquote><blockquote type="cite">e-mail and delete and destroy the message.<br></blockquote><blockquote type="cite">------------------------------------------------------------------------<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">_______________________________________________<br></blockquote><blockquote type="cite">gradsusr mailing list<br></blockquote><blockquote type="cite"><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br></blockquote><blockquote type="cite"><a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a><br></blockquote><blockquote type="cite"><br></blockquote><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></blockquote></div><br><div apple-content-edited="true">
<span class="Apple-style-span" style="font-size: 12px; "><div>--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span>
</div>
<br></body></html>