On Tue, Feb 15, 2011 at 8:01 PM, Charles Seman <span dir="ltr"><<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear GrADS Users,<br>
<br>
We're finding that GrADS 2.0.a7.oga.3 (OpenGrADS) is giving better results on our x86_64 analysis cluster than on an i686 workstation...<br>
<br>
The script <a href="http://climatology.gs" target="_blank">climatology.gs</a> (attached) works on the analysis cluster, but we get a "box_ave error." for one of the regridding commands on the workstation... Please find attached the GrADS script <a href="http://climatology.gs" target="_blank">climatology.gs</a>, the netCDF dataset used by the script, and output files from running the script on the two computers: climatology.output.i686 is from running on the workstation and contains the "box_ave error." for one of the regridding commands, </blockquote>
<div><br></div><div>Could it be a memory problem? The 32-bit binaries have some limitation on the amount of memory it can use.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
while climatology.output.x86_64 is from running on the analysis cluster and does not give this error (but gives "/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads: /opt/ifort/11.1.073/lib/intel64/libimf.so: no version information available (required by /app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads)"; is this something to be concerned about?).<br>
<br></blockquote><div><br></div><div>You may be using a version of libimf that is not quite compatible (most likely, it is getting it from your intel compiler while the opengrads extensions are built with gfortran.). To be sure, move the libimf.so from</div>
<div><br></div><div> /app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/libs</div><div><br></div><div>into</div><div><br></div><div>/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/gex</div><div>
<br></div><div>and this message should disappear.</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks for all your help, and please let me know if you need more information and/or have any questions about this issue.<br>
<br></blockquote><div><br></div><div>Let me know if this fixes your problem.</div><div><br></div><div> Arlindo</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks,<br>
Chuck<br>
<br>
-- <br>
<br>
Please note that <a href="mailto:Charles.Seman@noaa.gov" target="_blank">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>
email address, not <a href="mailto:cjs@gfdl.noaa.gov" target="_blank">cjs@gfdl.noaa.gov</a>.<br>
<br>
********************************************************************<br>
Charles Seman <a href="mailto:Charles.Seman@noaa.gov" target="_blank">Charles.Seman@noaa.gov</a><br>
U.S. Department of Commerce / NOAA / OAR<br>
Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547<br>
201 Forrestal Road fax: (609) 987-5063<br>
Princeton, NJ 08540-6649 <a href="http://www.gfdl.noaa.gov/~cjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
********************************************************************<br>
<br>
"The contents of this message are mine personally and do not reflect any<br>
official or unofficial position of the United States Federal Government,<br>
the United States Department of Commerce, or NOAA."<br>
<br>
<br>'reinit'<br>
'q config'<br>
say<br>
say '************************************************************************'<br>
say result<br>
say '************************************************************************'<br>
say<br>
'sdfopen <a href="http://climatology.nc" target="_blank">climatology.nc</a>'<br>
'q file'<br>
say result<br>
say '************************************************************************'<br>
say ' lat/lon unspecified'<br>
say '************************************************************************'<br>
'q dims'<br>
say result<br>
'set z 1'<br>
'set t 1'<br>
'define fieldtemp1 = ave(olr,t=1,t=12)'<br>
'define avemodann = aave(fieldtemp1,g)'<br>
'd avemodann'<br>
say 'global annual averages'<br>
say 'no regrid'<br>
say result<br>
'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,0.,-90)'<br>
'd aave(modelolrann,g)'<br>
say result<br>
'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,1.25,-88.75)'<br>
'd aave(modelolrann,g)'<br>
say result<br>
say '************************************************************************'<br>
say ' lat/lon specified'<br>
say '************************************************************************'<br>
'set x 1 144'<br>
'set y 1 90'<br>
'set z 1'<br>
'set t 1'<br>
'q dims'<br>
say result<br>
'define fieldtemp1 = ave(olr,t=1,t=12)'<br>
'define avemodann = aave(fieldtemp1,g)'<br>
'd avemodann'<br>
say 'global annual averages'<br>
say 'no regrid'<br>
say result<br>
'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,0.,-90)'<br>
'd aave(modelolrann,g)'<br>
say result<br>
'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,1.25,-88.75)'<br>
'd aave(modelolrann,g)'<br>
say result<br>
<br>Linux cjs 2.6.18-194.26.1.el5PAE #1 SMP Fri Oct 29 14:28:58 EDT 2010 i686 i686 i386 GNU/Linux<br>
<br>
Welcome to the OpenGrADS Bundle Distribution<br>
--------------------------------------------<br>
<br>
For additional information enter "grads -h".<br>
<br>
Starting "/usr/local/grads-2.0.a7.oga.3/Contents/Linux/Versions/2.0.a7.oga.3/i686/grads -bclx <a href="http://climatology.gs" target="_blank">climatology.gs</a> " ...<br>
<br>
<br>
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3<br>
Copyright (c) 1988-2008 by Brian Doty and the<br>
Institute for Global Environment and Society (IGES)<br>
GrADS comes with ABSOLUTELY NO WARRANTY<br>
See file COPYRIGHT for more information<br>
<br>
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff<br>
Issue 'q config' command for more information.<br>
Loading User Defined Extensions table </usr/local/grads-2.0.a7.oga.3/Contents/Linux/Versions/2.0.a7.oga.3/i686/gex/udxt> ... ok.<br>
GX Package Initialization: Size = 11 8.5<br>
Running in Batch mode<br>
No hardcopy metafile open<br>
All files closed; all defined objects released;<br>
All GrADS attributes have been reinitialized<br>
<br>
************************************************************************<br>
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff<br>
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3<br>
Copyright (c) 1988-2009 by Brian Doty and the<br>
Institute for Global Environment and Society (IGES)<br>
This program is distributed WITHOUT ANY WARRANTY<br>
See file COPYRIGHT for more information.<br>
<br>
Built Thu Oct 29 17:28:29 EDT 2009 for i686-pc-linux-gnu<br>
<br>
This version of GrADS has been configured with the following options:<br>
o Built on a LITTLE ENDIAN machine<br>
o Command line editing ENABLED<br>
<a href="http://tiswww.case.edu/php/chet/readline/rltop.html" target="_blank">http://tiswww.case.edu/php/chet/readline/rltop.html</a><br>
o printim command for image output ENABLED<br>
<a href="http://www.zlib.net" target="_blank">http://www.zlib.net</a><br>
<a href="http://www.libpng.org/pub/png/libpng.html" target="_blank">http://www.libpng.org/pub/png/libpng.html</a><br>
<a href="http://www.libgd.org/Main_Page" target="_blank">http://www.libgd.org/Main_Page</a><br>
o GRIB2 interface ENABLED<br>
<a href="http://www.ijg.org" target="_blank">http://www.ijg.org</a><br>
<a href="http://www.ece.uvic.ca/~mdadams/jasper" target="_blank">http://www.ece.uvic.ca/~mdadams/jasper</a><br>
<a href="http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2" target="_blank">http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2</a><br>
g2clib-1.0.5<br>
o NetCDF interface ENABLED<br>
<a href="http://www.opendap.org" target="_blank">http://www.opendap.org</a><br>
libnc-dap 4.0.1-beta3-snapshot2009021712 of Feb 21 2009 14:55:15 $<br>
o HDF interface ENABLED<br>
<a href="http://hdfgroup.org" target="_blank">http://hdfgroup.org</a><br>
HDF 4.2r3<br>
HDF5 1.8.2<br>
o Athena Widget GUI ENABLED<br>
o OPeNDAP gridded data interface ENABLED<br>
<a href="http://www.opendap.org" target="_blank">http://www.opendap.org</a><br>
libdap 3.7.10<br>
o OPeNDAP station data interface ENABLED<br>
<a href="http://iges.org/grads/gadoc/supplibs.html" target="_blank">http://iges.org/grads/gadoc/supplibs.html</a><br>
libgadap 2.0.oga.1<br>
o GeoTIFF and KML output ENABLED<br>
<a href="http://www.libtiff.org" target="_blank">http://www.libtiff.org</a><br>
<a href="http://geotiff.osgeo.org" target="_blank">http://geotiff.osgeo.org</a><br>
<br>
For additional information please consult <a href="http://iges.org/grads" target="_blank">http://iges.org/grads</a><br>
<br>
<br>
************************************************************************<br>
<br>
File 1 : c48L48_am3p9<br>
Descriptor: <a href="http://climatology.nc" target="_blank">climatology.nc</a><br>
Binary: <a href="http://climatology.nc" target="_blank">climatology.nc</a><br>
Type = Gridded<br>
Xsize = 144 Ysize = 90 Zsize = 23 Tsize = 12 Esize = 1<br>
Number of Variables = 5<br>
swdn_toa 0 t,y,x SW flux down at TOA<br>
swup_toa 0 t,y,x SW flux up at TOA<br>
olr 0 t,y,x outgoing longwave radiation<br>
olr_clr 0 t,y,x clear skyoutgoing longwave radiation<br>
swup_toa_clr 0 t,y,x clear skySW flux up at TOA<br>
<br>
************************************************************************<br>
lat/lon unspecified<br>
************************************************************************<br>
Default file number is: 1<br>
X is varying Lon = 0 to 360 X = 0.5 to 144.5<br>
Y is varying Lat = -89 to 89 Y = 1 to 90<br>
Z is fixed Lev = 1000 Z = 1<br>
T is fixed Time = 12Z16JAN2000 T = 1<br>
E is fixed Ens = 1 E = 1<br>
<br>
global annual averages<br>
no regrid<br>
Result value = 235.644<br>
<br>
145 360 0.<br>
re(fieldtemp1,145,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)<br>
Result value = 235.651<br>
<br>
145 360 1.25<br>
re(fieldtemp1,145,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)<br>
box_ave error.<br>
Result value = 142.041<br>
<br>
************************************************************************<br>
lat/lon specified<br>
************************************************************************<br>
Default file number is: 1<br>
X is varying Lon = 1.25 to 358.75 X = 1 to 144<br>
Y is varying Lat = -89 to 89 Y = 1 to 90<br>
Z is fixed Lev = 1000 Z = 1<br>
T is fixed Time = 12Z16JAN2000 T = 1<br>
E is fixed Ens = 1 E = 1<br>
<br>
global annual averages<br>
no regrid<br>
Result value = 235.644<br>
<br>
144 358.75 0.<br>
re(fieldtemp1,144,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)<br>
Result value = 235.613<br>
<br>
144 358.75 1.25<br>
re(fieldtemp1,144,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)<br>
Result value = 235.651<br>
<br>
<br>Linux an003 2.6.18-194.26.1.el5 #1 SMP Tue Nov 9 12:54:20 EST 2010 x86_64 x86_64 x86_64 GNU/Linux<br>
<br>
Welcome to the OpenGrADS Bundle Distribution<br>
--------------------------------------------<br>
<br>
For additional information enter "grads -h".<br>
<br>
Starting "/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads -bclx <a href="http://climatology.gs" target="_blank">climatology.gs</a> " ...<br>
<br>
/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads: /opt/ifort/11.1.073/lib/intel64/libimf.so: no version information available (required by /app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads)<br>
<br>
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3<br>
Copyright (c) 1988-2008 by Brian Doty and the<br>
Institute for Global Environment and Society (IGES)<br>
GrADS comes with ABSOLUTELY NO WARRANTY<br>
See file COPYRIGHT for more information<br>
<br>
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff<br>
Issue 'q config' command for more information.<br>
Loading User Defined Extensions table </app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/gex/udxt> ... ok.<br>
GX Package Initialization: Size = 11 8.5<br>
Running in Batch mode<br>
No hardcopy metafile open<br>
All files closed; all defined objects released;<br>
All GrADS attributes have been reinitialized<br>
<br>
************************************************************************<br>
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff<br>
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3<br>
Copyright (c) 1988-2009 by Brian Doty and the<br>
Institute for Global Environment and Society (IGES)<br>
This program is distributed WITHOUT ANY WARRANTY<br>
See file COPYRIGHT for more information.<br>
<br>
Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu<br>
<br>
This version of GrADS has been configured with the following options:<br>
o Built on a LITTLE ENDIAN machine<br>
o Command line editing ENABLED<br>
<a href="http://tiswww.case.edu/php/chet/readline/rltop.html" target="_blank">http://tiswww.case.edu/php/chet/readline/rltop.html</a><br>
o printim command for image output ENABLED<br>
<a href="http://www.zlib.net" target="_blank">http://www.zlib.net</a><br>
<a href="http://www.libpng.org/pub/png/libpng.html" target="_blank">http://www.libpng.org/pub/png/libpng.html</a><br>
<a href="http://www.libgd.org/Main_Page" target="_blank">http://www.libgd.org/Main_Page</a><br>
o GRIB2 interface ENABLED<br>
<a href="http://www.ijg.org" target="_blank">http://www.ijg.org</a><br>
<a href="http://www.ece.uvic.ca/~mdadams/jasper" target="_blank">http://www.ece.uvic.ca/~mdadams/jasper</a><br>
<a href="http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2" target="_blank">http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2</a><br>
g2clib-1.0.5<br>
o NetCDF interface ENABLED<br>
<a href="http://www.opendap.org" target="_blank">http://www.opendap.org</a><br>
libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar 3 2009 14:13:33 $<br>
o HDF interface ENABLED<br>
<a href="http://hdfgroup.org" target="_blank">http://hdfgroup.org</a><br>
HDF 4.2r3<br>
HDF5 1.8.2<br>
o Athena Widget GUI ENABLED<br>
o OPeNDAP gridded data interface ENABLED<br>
<a href="http://www.opendap.org" target="_blank">http://www.opendap.org</a><br>
libdap 3.7.10<br>
o OPeNDAP station data interface ENABLED<br>
<a href="http://iges.org/grads/gadoc/supplibs.html" target="_blank">http://iges.org/grads/gadoc/supplibs.html</a><br>
libgadap 2.0.oga.1<br>
o GeoTIFF and KML output ENABLED<br>
<a href="http://www.libtiff.org" target="_blank">http://www.libtiff.org</a><br>
<a href="http://geotiff.osgeo.org" target="_blank">http://geotiff.osgeo.org</a><br>
<br>
For additional information please consult <a href="http://iges.org/grads" target="_blank">http://iges.org/grads</a><br>
<br>
<br>
************************************************************************<br>
<br>
File 1 : c48L48_am3p9<br>
Descriptor: <a href="http://climatology.nc" target="_blank">climatology.nc</a><br>
Binary: <a href="http://climatology.nc" target="_blank">climatology.nc</a><br>
Type = Gridded<br>
Xsize = 144 Ysize = 90 Zsize = 23 Tsize = 12 Esize = 1<br>
Number of Variables = 5<br>
swdn_toa 0 t,y,x SW flux down at TOA<br>
swup_toa 0 t,y,x SW flux up at TOA<br>
olr 0 t,y,x outgoing longwave radiation<br>
olr_clr 0 t,y,x clear skyoutgoing longwave radiation<br>
swup_toa_clr 0 t,y,x clear skySW flux up at TOA<br>
<br>
************************************************************************<br>
lat/lon unspecified<br>
************************************************************************<br>
Default file number is: 1<br>
X is varying Lon = 0 to 360 X = 0.5 to 144.5<br>
Y is varying Lat = -89 to 89 Y = 1 to 90<br>
Z is fixed Lev = 1000 Z = 1<br>
T is fixed Time = 12Z16JAN2000 T = 1<br>
E is fixed Ens = 1 E = 1<br>
<br>
global annual averages<br>
no regrid<br>
Result value = 235.644<br>
<br>
145 360 0.<br>
re(fieldtemp1,145,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)<br>
Result value = 235.651<br>
<br>
145 360 1.25<br>
re(fieldtemp1,145,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)<br>
Result value = 235.651<br>
<br>
************************************************************************<br>
lat/lon specified<br>
************************************************************************<br>
Default file number is: 1<br>
X is varying Lon = 1.25 to 358.75 X = 1 to 144<br>
Y is varying Lat = -89 to 89 Y = 1 to 90<br>
Z is fixed Lev = 1000 Z = 1<br>
T is fixed Time = 12Z16JAN2000 T = 1<br>
E is fixed Ens = 1 E = 1<br>
<br>
global annual averages<br>
no regrid<br>
Result value = 235.644<br>
<br>
144 358.75 0.<br>
re(fieldtemp1,144,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)<br>
Result value = 235.613<br>
<br>
144 358.75 1.25<br>
re(fieldtemp1,144,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)<br>
Result value = 235.651<br>
<br>
<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><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>