[gradsusr] Different output values from Grads and wgrib2

[Wesley Ebisuzaki]

Bo,

   I compared your ctl file with mine and they were the same except for 
the location of the files (dset/index/comments).

   Maybe it's a problem with gribmap.  Perhaps gribmap is confusing the 
30mb TMP with another field.  Try running

   gribmap -v -i test.ctl


              Wesley


Bo Zhang wrote:
> Hi Wesley,
>
>     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.
>
>     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). 
>   
>
> Bo Zhang
> Dow 824,Civil and Environmental Engineering
> Michigan Tech University
> 1400 Townsends Drive. 
> Houghton, MI, 49931
>
> ----- Original Message -----
> From: "Wesley Ebisuzaki" <Wesley.Ebisuzaki at noaa.gov>
> To: "Bo Zhang" <bzhang3 at mtu.edu>
> Sent: Thursday, October 20, 2011 4:40:13 PM (GMT-0500) Auto-Detected
> Subject: Re: Different output values from Grads and wgrib2
>
> Bo,
>
>     With your file, I get the right number for 30 mb with GrADS 
> installed on my
> machine. (64-bit RedHat 5 - linux, GrADS 2.0a9 from the COLA site).
>
>     The problem that you are finding may be version / machine specific.  You
> should send to the GrADS list
>
>    (1) OS
>    (2) version of GrADS
>    (3) who compiled the GrADS or where obtained from
>
>              Wesley
>
>
>
> Bo Zhang wrote:
>   
>> Hi Wesley,
>>
>> 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. 
>>
>> With wgrib2:
>> [bzhang3 at flexruns wgrib2]$ ./wgrib2 /home/campus03/bzhang3/Grads/fnl_20100508_06_00 -match TMP -s -ijlat 181 91
>> 3:120296:d=2010050806:TMP:10 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=239
>> 8:393851:d=2010050806:TMP:20 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=220.5
>> 13:666758:d=2010050806:TMP:30 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=215.5
>> 18:961202:d=2010050806:TMP:50 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=204.5
>> 23:1271647:d=2010050806:TMP:70 mb:anl::(181,91),lon=180.000000,lat=0.000000,val=193.7
>>
>> With GrADS
>> [bzhang3 at flexruns bin]$ ./grads
>> ga-> set lat 0
>> LAT set to 0 0 
>> ga-> set lon 180
>> LON set to 180 180
>> ga-> d tmpprs (lev=10)
>> Result value = 239  
>> ga-> d tmpprs (lev=20)
>> Result value = 220.5 
>> ga-> d tmpprs (lev=30)
>> Result value = 300.9 
>> ga-> d tmpprs (lev=50)
>> Result value = 204.5 
>>
>> 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:
>>
>> Grid Analysis and Display System (GrADS) Version 2.0.0
>> Copyright (c) 1988-2011 by Brian Doty and the
>> Institute for Global Environment and Society (IGES)
>> GrADS comes with ABSOLUTELY NO WARRANTY
>> See file COPYRIGHT for more information
>> Config: v2.0.0 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile
>>
>> Thanks.
>>
>> Bo Zhang
>> Dow 824,Civil and Environmental Engineering
>> Michigan Tech University
>> 1400 Townsends Drive. 
>> Houghton, MI, 49931
>>     
>
>   

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: test.ctl
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20111021/bf1a61b0/attachment-0003.pl