<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>I think the data file is corrupted. I can duplicate the error, but the other diagnostic tools (wgrib2 and grib2scan) throw errors when trying to read that file. </div><div>--Jennifer</div><br><div><div>On Oct 21, 2011, at 1:17 PM, Wesley Ebisuzaki wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Bo,<br><br> I compared your ctl file with mine and they were the same except for the location of the files (dset/index/comments).<br><br> Maybe it's a problem with gribmap. Perhaps gribmap is confusing the 30mb TMP with another field. Try running<br><br> gribmap -v -i test.ctl<br><br><br> Wesley<br><br><br>Bo Zhang wrote:<br><blockquote type="cite">Hi Wesley,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Since John tried and found the same problem. Perhaps, some issue happens during control files generation. Could you please send me the control files you create and may be g2ctl also.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> I have tried both GrADS2.0 and 2.0.a9 on my machine which all have this issue. My OS: Red Hat Enterprise Linux Server release 5.6 (Tikanga). <br></blockquote><blockquote type="cite">Bo Zhang<br></blockquote><blockquote type="cite">Dow 824,Civil and Environmental Engineering<br></blockquote><blockquote type="cite">Michigan Tech University<br></blockquote><blockquote type="cite">1400 Townsends Drive. Houghton, MI, 49931<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">----- Original Message -----<br></blockquote><blockquote type="cite">From: "Wesley Ebisuzaki" <<a href="mailto:Wesley.Ebisuzaki@noaa.gov">Wesley.Ebisuzaki@noaa.gov</a>><br></blockquote><blockquote type="cite">To: "Bo Zhang" <<a href="mailto:bzhang3@mtu.edu">bzhang3@mtu.edu</a>><br></blockquote><blockquote type="cite">Sent: Thursday, October 20, 2011 4:40:13 PM (GMT-0500) Auto-Detected<br></blockquote><blockquote type="cite">Subject: Re: Different output values from Grads and wgrib2<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Bo,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> With your file, I get the right number for 30 mb with GrADS installed on my<br></blockquote><blockquote type="cite">machine. (64-bit RedHat 5 - linux, GrADS 2.0a9 from the COLA site).<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> The problem that you are finding may be version / machine specific. You<br></blockquote><blockquote type="cite">should send to the GrADS list<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> (1) OS<br></blockquote><blockquote type="cite"> (2) version of GrADS<br></blockquote><blockquote type="cite"> (3) who compiled the GrADS or where obtained from<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Wesley<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Bo Zhang wrote:<br></blockquote><blockquote type="cite"> <br></blockquote><blockquote type="cite"><blockquote type="cite">Hi Wesley,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">It seems my space is not working correctly, the file broken when John downloaded it. So, excuse me for sending you the data file directly. Here is my results of running GrADS on this data file. <br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">With wgrib2:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">[bzhang3@flexruns wgrib2]$ ./wgrib2 /home/campus03/bzhang3/Grads/fnl_20100508_06_00 -match TMP -s -ijlat 181 91<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">3:120296:d=2010050806:TMP:10 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=239<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">8:393851:d=2010050806:TMP:20 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=220.5<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">13:666758:d=2010050806:TMP:30 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=215.5<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">18:961202:d=2010050806:TMP:50 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=204.5<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">23:1271647:d=2010050806:TMP:70 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=193.7<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">With GrADS<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">[bzhang3@flexruns bin]$ ./grads<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">ga-> set lat 0<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">LAT set to 0 0 ga-> set lon 180<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">LON set to 180 180<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">ga-> d tmpprs (lev=10)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Result value = 239 ga-> d tmpprs (lev=20)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Result value = 220.5 ga-> d tmpprs (lev=30)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Result value = 300.9 ga-> d tmpprs (lev=50)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Result value = 204.5 <br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">All numbers match except at level of 30mb. I don't think the problem is from the data file since wgrib2 gives correct numbers. I use g2ctl and gribmap packed in GrADS to generate .ctl and .idx files. I am not sure my GrADS version, probably CentOS 5.6. Information on GrADS:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Grid Analysis and Display System (GrADS) Version 2.0.0<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Copyright (c) 1988-2011 by Brian Doty and the<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Institute for Global Environment and Society (IGES)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">GrADS comes with ABSOLUTELY NO WARRANTY<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">See file COPYRIGHT for more information<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Config: v2.0.0 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Thanks.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Bo Zhang<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Dow 824,Civil and Environmental Engineering<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Michigan Tech University<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">1400 Townsends Drive. Houghton, MI, 49931<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"> <br></blockquote></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> <br></blockquote><br>dset ^fnl_20100508_06_00<br>index ^fnl_20100508_06_00.idx<br>undef 9.999E+20<br>title fnl_20100508_06_00<br>* produced by g2ctl v0.0.8.3<br>* command line options: fnl_20100508_06_00<br>* griddef=1:0:(360 x 181):grid_template=0:winds(N/S): <br>dtype grib2<br>ydef 181 linear -90.000000 1<br>xdef 360 linear 0.000000 1.000000<br>tdef 1 linear 06Z08may2010 1mo<br>* PROFILE hPa<br>zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000<br>options pascals<br>vars 1<br>TMPprs 26,100 0,0,0 ** (1000 975 950 925 900.. 70 50 30 20 10) Temperature [K]<br>ENDVARS<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="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><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></span></span>
</div>
<br></body></html>