[gradsusr] OpenGrADS: GrADS 2.0.a7.oga.3 regridding...
Arlindo da Silva
dasilva at alum.mit.edu
Wed Feb 16 20:24:03 EST 2011
On Tue, Feb 15, 2011 at 8:01 PM, Charles Seman <Charles.Seman at noaa.gov>wrote:
> Dear GrADS Users,
>
> 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...
>
> The script climatology.gs (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 climatology.gs, 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,
Could it be a memory problem? The 32-bit binaries have some limitation on
the amount of memory it can use.
> 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?).
>
>
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
/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/libs
into
/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/gex
and this message should disappear.
> Thanks for all your help, and please let me know if you need more
> information and/or have any questions about this issue.
>
>
Let me know if this fixes your problem.
Arlindo
> Thanks,
> Chuck
>
> --
>
> Please note that Charles.Seman at noaa.gov should be considered my NOAA
> email address, not cjs at gfdl.noaa.gov.
>
> ********************************************************************
> Charles Seman Charles.Seman at noaa.gov
> U.S. Department of Commerce / NOAA / OAR
> Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
> 201 Forrestal Road fax: (609) 987-5063
> Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
> ********************************************************************
>
> "The contents of this message are mine personally and do not reflect any
> official or unofficial position of the United States Federal Government,
> the United States Department of Commerce, or NOAA."
>
>
> 'reinit'
> 'q config'
> say
> say
> '************************************************************************'
> say result
> say
> '************************************************************************'
> say
> 'sdfopen climatology.nc'
> 'q file'
> say result
> say
> '************************************************************************'
> say ' lat/lon unspecified'
> say
> '************************************************************************'
> 'q dims'
> say result
> 'set z 1'
> 'set t 1'
> 'define fieldtemp1 = ave(olr,t=1,t=12)'
> 'define avemodann = aave(fieldtemp1,g)'
> 'd avemodann'
> say 'global annual averages'
> say 'no regrid'
> say result
> 'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,0.,-90)'
> 'd aave(modelolrann,g)'
> say result
> 'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,1.25,-88.75)'
> 'd aave(modelolrann,g)'
> say result
> say
> '************************************************************************'
> say ' lat/lon specified'
> say
> '************************************************************************'
> 'set x 1 144'
> 'set y 1 90'
> 'set z 1'
> 'set t 1'
> 'q dims'
> say result
> 'define fieldtemp1 = ave(olr,t=1,t=12)'
> 'define avemodann = aave(fieldtemp1,g)'
> 'd avemodann'
> say 'global annual averages'
> say 'no regrid'
> say result
> 'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,0.,-90)'
> 'd aave(modelolrann,g)'
> say result
> 'define modelolrann = regrid2(fieldtemp1,2.5,2.5,ba_p1,1.25,-88.75)'
> 'd aave(modelolrann,g)'
> say result
>
> Linux cjs 2.6.18-194.26.1.el5PAE #1 SMP Fri Oct 29 14:28:58 EDT 2010 i686
> i686 i386 GNU/Linux
>
> Welcome to the OpenGrADS Bundle Distribution
> --------------------------------------------
>
> For additional information enter "grads -h".
>
> Starting
> "/usr/local/grads-2.0.a7.oga.3/Contents/Linux/Versions/2.0.a7.oga.3/i686/grads
> -bclx climatology.gs " ...
>
>
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2008 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.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds
> hdf5 opendap-grids,stn athena geotiff
> Issue 'q config' command for more information.
> 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.
> GX Package Initialization: Size = 11 8.5
> Running in Batch mode
> No hardcopy metafile open
> All files closed; all defined objects released;
> All GrADS attributes have been reinitialized
>
> ************************************************************************
> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds
> hdf5 opendap-grids,stn athena geotiff
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2009 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> This program is distributed WITHOUT ANY WARRANTY
> See file COPYRIGHT for more information.
>
> Built Thu Oct 29 17:28:29 EDT 2009 for i686-pc-linux-gnu
>
> This version of GrADS has been configured with the following options:
> o Built on a LITTLE ENDIAN machine
> o Command line editing ENABLED
> http://tiswww.case.edu/php/chet/readline/rltop.html
> o printim command for image output ENABLED
> http://www.zlib.net
> http://www.libpng.org/pub/png/libpng.html
> http://www.libgd.org/Main_Page
> o GRIB2 interface ENABLED
> http://www.ijg.org
> http://www.ece.uvic.ca/~mdadams/jasper
> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
> g2clib-1.0.5
> o NetCDF interface ENABLED
> http://www.opendap.org
> libnc-dap 4.0.1-beta3-snapshot2009021712 of Feb 21 2009 14:55:15 $
> o HDF interface ENABLED
> http://hdfgroup.org
> HDF 4.2r3
> HDF5 1.8.2
> o Athena Widget GUI ENABLED
> o OPeNDAP gridded data interface ENABLED
> http://www.opendap.org
> libdap 3.7.10
> o OPeNDAP station data interface ENABLED
> http://iges.org/grads/gadoc/supplibs.html
> libgadap 2.0.oga.1
> o GeoTIFF and KML output ENABLED
> http://www.libtiff.org
> http://geotiff.osgeo.org
>
> For additional information please consult http://iges.org/grads
>
>
> ************************************************************************
>
> File 1 : c48L48_am3p9
> Descriptor: climatology.nc
> Binary: climatology.nc
> Type = Gridded
> Xsize = 144 Ysize = 90 Zsize = 23 Tsize = 12 Esize = 1
> Number of Variables = 5
> swdn_toa 0 t,y,x SW flux down at TOA
> swup_toa 0 t,y,x SW flux up at TOA
> olr 0 t,y,x outgoing longwave radiation
> olr_clr 0 t,y,x clear skyoutgoing longwave radiation
> swup_toa_clr 0 t,y,x clear skySW flux up at TOA
>
> ************************************************************************
> lat/lon unspecified
> ************************************************************************
> Default file number is: 1
> X is varying Lon = 0 to 360 X = 0.5 to 144.5
> Y is varying Lat = -89 to 89 Y = 1 to 90
> Z is fixed Lev = 1000 Z = 1
> T is fixed Time = 12Z16JAN2000 T = 1
> E is fixed Ens = 1 E = 1
>
> global annual averages
> no regrid
> Result value = 235.644
>
> 145 360 0.
> re(fieldtemp1,145,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)
> Result value = 235.651
>
> 145 360 1.25
> re(fieldtemp1,145,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)
> box_ave error.
> Result value = 142.041
>
> ************************************************************************
> lat/lon specified
> ************************************************************************
> Default file number is: 1
> X is varying Lon = 1.25 to 358.75 X = 1 to 144
> Y is varying Lat = -89 to 89 Y = 1 to 90
> Z is fixed Lev = 1000 Z = 1
> T is fixed Time = 12Z16JAN2000 T = 1
> E is fixed Ens = 1 E = 1
>
> global annual averages
> no regrid
> Result value = 235.644
>
> 144 358.75 0.
> re(fieldtemp1,144,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)
> Result value = 235.613
>
> 144 358.75 1.25
> re(fieldtemp1,144,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)
> Result value = 235.651
>
>
> 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
>
> Welcome to the OpenGrADS Bundle Distribution
> --------------------------------------------
>
> For additional information enter "grads -h".
>
> Starting
> "/app/grads-2.0.a7.oga.3/bin/Linux/Versions/2.0.a7.oga.3/x86_64/grads
> -bclx climatology.gs " ...
>
> /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)
>
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2008 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.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds
> hdf5 opendap-grids,stn athena geotiff
> Issue 'q config' command for more information.
> 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.
> GX Package Initialization: Size = 11 8.5
> Running in Batch mode
> No hardcopy metafile open
> All files closed; all defined objects released;
> All GrADS attributes have been reinitialized
>
> ************************************************************************
> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds
> hdf5 opendap-grids,stn athena geotiff
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2009 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> This program is distributed WITHOUT ANY WARRANTY
> See file COPYRIGHT for more information.
>
> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
>
> This version of GrADS has been configured with the following options:
> o Built on a LITTLE ENDIAN machine
> o Command line editing ENABLED
> http://tiswww.case.edu/php/chet/readline/rltop.html
> o printim command for image output ENABLED
> http://www.zlib.net
> http://www.libpng.org/pub/png/libpng.html
> http://www.libgd.org/Main_Page
> o GRIB2 interface ENABLED
> http://www.ijg.org
> http://www.ece.uvic.ca/~mdadams/jasper
> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
> g2clib-1.0.5
> o NetCDF interface ENABLED
> http://www.opendap.org
> libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar 3 2009 14:13:33 $
> o HDF interface ENABLED
> http://hdfgroup.org
> HDF 4.2r3
> HDF5 1.8.2
> o Athena Widget GUI ENABLED
> o OPeNDAP gridded data interface ENABLED
> http://www.opendap.org
> libdap 3.7.10
> o OPeNDAP station data interface ENABLED
> http://iges.org/grads/gadoc/supplibs.html
> libgadap 2.0.oga.1
> o GeoTIFF and KML output ENABLED
> http://www.libtiff.org
> http://geotiff.osgeo.org
>
> For additional information please consult http://iges.org/grads
>
>
> ************************************************************************
>
> File 1 : c48L48_am3p9
> Descriptor: climatology.nc
> Binary: climatology.nc
> Type = Gridded
> Xsize = 144 Ysize = 90 Zsize = 23 Tsize = 12 Esize = 1
> Number of Variables = 5
> swdn_toa 0 t,y,x SW flux down at TOA
> swup_toa 0 t,y,x SW flux up at TOA
> olr 0 t,y,x outgoing longwave radiation
> olr_clr 0 t,y,x clear skyoutgoing longwave radiation
> swup_toa_clr 0 t,y,x clear skySW flux up at TOA
>
> ************************************************************************
> lat/lon unspecified
> ************************************************************************
> Default file number is: 1
> X is varying Lon = 0 to 360 X = 0.5 to 144.5
> Y is varying Lat = -89 to 89 Y = 1 to 90
> Z is fixed Lev = 1000 Z = 1
> T is fixed Time = 12Z16JAN2000 T = 1
> E is fixed Ens = 1 E = 1
>
> global annual averages
> no regrid
> Result value = 235.644
>
> 145 360 0.
> re(fieldtemp1,145,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)
> Result value = 235.651
>
> 145 360 1.25
> re(fieldtemp1,145,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)
> Result value = 235.651
>
> ************************************************************************
> lat/lon specified
> ************************************************************************
> Default file number is: 1
> X is varying Lon = 1.25 to 358.75 X = 1 to 144
> Y is varying Lat = -89 to 89 Y = 1 to 90
> Z is fixed Lev = 1000 Z = 1
> T is fixed Time = 12Z16JAN2000 T = 1
> E is fixed Ens = 1 E = 1
>
> global annual averages
> no regrid
> Result value = 235.644
>
> 144 358.75 0.
> re(fieldtemp1,144,LINEAR,0.,2.5,73,LINEAR,-90,2.5,ba)
> Result value = 235.613
>
> 144 358.75 1.25
> re(fieldtemp1,144,LINEAR,1.25,2.5,72,LINEAR,-88.75,2.5,ba)
> Result value = 235.651
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110216/3b40d3be/attachment-0003.html
More information about the gradsusr
mailing list