<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hello GrADS users,<br>
I am trying to use GrADS with OSCAR data downloaded from here (using
wget):<br>
<br>
<tt>wget
<a class="moz-txt-link-freetext" href="http://podaac-opendap.jpl.nasa.gov/opendap/allData/oscar/preview/L4/oscar_third_deg/oscar_vel7550.nc.gz.nc">http://podaac-opendap.jpl.nasa.gov/opendap/allData/oscar/preview/L4/oscar_third_deg/oscar_vel7550.nc.gz.nc</a></tt><tt><br>
</tt><br>
I am trying to extract lat/lon boxes from this data, so I am
invoking GrADS using the "lats4d.sh" wrapper script. When I try to
open the binary NetCDF file above, I get this error:<br>
<br>
<tt>gadsdf: SDF file has no discernable X coordinate.</tt><tt><br>
</tt><br>
Therefor I am using a "Data Descriptor" (.ctl) file. However, I am
unable to get that to work. Here is the error I am getting with my
.ctl file:<br>
<br>
<tt>gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds
dimension bound</tt><tt><br>
</tt><tt>Data Request Error: Error for variable 'u'</tt><tt><br>
</tt><tt> Error ocurred at column 1</tt><tt><br>
</tt><br>
I don't know what to try next. Here is my .ctl file for this data
-- the only thing that is a little weird is that the source says
"Longitude extends from 20 E to 420 E to avoid a break in major
ocean basins", so I have the XDEF line starting with the value
"20.0" instead of "0.0"... but I have also tried starting XDEF with
"0.0" and it doesn't change the error message:<br>
<br>
<tt>DSET ^oscar_vel7550.nc.gz.nc</tt><tt><br>
</tt><tt>DTYPE netcdf</tt><tt><br>
</tt><tt>XDEF 1201 LINEAR 20.0 0.3</tt><tt><br>
</tt><tt>YDEF 481 LINEAR -80.0 0.3333333333333333</tt><tt><br>
</tt><tt>ZDEF 1 LINEAR 0 1</tt><tt><br>
</tt><tt>tdef 1 linear 00Z7jun2013 1dy</tt><tt><br>
</tt><tt>UNDEF 0.0</tt><tt><br>
</tt><tt>VARS 2</tt><tt><br>
</tt><tt>u=>u 0 t,z,x,y Ocean Surface Zonal Currents [m/s]</tt><tt><br>
</tt><tt>v=>v 0 t,z,x,y Ocean Surface Meridional Currents [m/s]</tt><tt><br>
</tt><tt>ENDVARS</tt><tt><br>
</tt><br>
And below is the output of "ncdump -h", which shows the dimensions,
along with the full output from GrADS. Also, this source is
available via OPeNDAP, so you can find out additional details about
this source by visiting<br>
<br>
<meta http-equiv="content-type" content="text/html;
charset=ISO-8859-1">
<tt><a
href="http://podaac-opendap.jpl.nasa.gov/opendap/allData/oscar/preview/L4/oscar_third_deg/oscar_vel7550.nc.gz.html">http://podaac-opendap.jpl.nasa.gov/opendap/allData/oscar/preview/L4/oscar_third_deg/oscar_vel7550.nc.gz.html</a></tt><br>
<br>
Any advice or help is greatly appreciated. Here is the "ncdump -h"
output and the "lats4d.sh" GrADS output (which also shows I am using
GrADS version v2.0.1.oga.1):<br>
<br>
<tt>netcdf oscar_vel7550.nc.gz {</tt><tt><br>
</tt><tt>dimensions:</tt><tt><br>
</tt><tt> time = 1 ;</tt><tt><br>
</tt><tt> year = 1 ;</tt><tt><br>
</tt><tt> depth = 1 ;</tt><tt><br>
</tt><tt> latitude = 481 ;</tt><tt><br>
</tt><tt> longitude = 1201 ;</tt><tt><br>
</tt><tt>variables:</tt><tt><br>
</tt><tt> int time(time) ;</tt><tt><br>
</tt><tt> time:units = "day since 1992-10-05 00:00:00" ;</tt><tt><br>
</tt><tt> time:long_name = "Day since 1992-10-05 00:00:00" ;</tt><tt><br>
</tt><tt> float year(year) ;</tt><tt><br>
</tt><tt> year:units = "time in years" ;</tt><tt><br>
</tt><tt> year:long_name = "Time in fractional year" ;</tt><tt><br>
</tt><tt> float depth(depth) ;</tt><tt><br>
</tt><tt> depth:units = "meter" ;</tt><tt><br>
</tt><tt> depth:long_name = "Depth" ;</tt><tt><br>
</tt><tt> double latitude(latitude) ;</tt><tt><br>
</tt><tt> latitude:units = "degrees-north" ;</tt><tt><br>
</tt><tt> latitude:long_name = "Latitude" ;</tt><tt><br>
</tt><tt> double longitude(longitude) ;</tt><tt><br>
</tt><tt> longitude:units = "degrees-east" ;</tt><tt><br>
</tt><tt> longitude:long_name = "Longitude" ;</tt><tt><br>
</tt><tt> double u(time, depth, latitude, longitude) ;</tt><tt><br>
</tt><tt> u:units = "meter/sec" ;</tt><tt><br>
</tt><tt> u:long_name = "Ocean Surface Zonal Currents" ;</tt><tt><br>
</tt><tt> u:missing_value = 0. ;</tt><tt><br>
</tt><tt> double v(time, depth, latitude, longitude) ;</tt><tt><br>
</tt><tt> v:units = "meter/sec" ;</tt><tt><br>
</tt><tt> v:long_name = "Ocean Surface Meridional Currents" ;</tt><tt><br>
</tt><tt> v:missing_value = 0. ;</tt><tt><br>
</tt><tt> double um(time, depth, latitude, longitude) ;</tt><tt><br>
</tt><tt> um:units = "meter/sec" ;</tt><tt><br>
</tt><tt> um:long_name = "Ocean Surface Zonal Currents
Maximum Mask" ;</tt><tt><br>
</tt><tt> um:missing_value = 0. ;</tt><tt><br>
</tt><tt> double vm(time, depth, latitude, longitude) ;</tt><tt><br>
</tt><tt> vm:units = "meter/sec" ;</tt><tt><br>
</tt><tt> vm:long_name = "Ocean Surface Meridional Currents
Maximum Mask" ;</tt><tt><br>
</tt><tt> vm:missing_value = 0. ;</tt><tt><br>
</tt><tt><br>
</tt><tt>// global attributes:</tt><tt><br>
</tt><tt> :NC_GLOBAL.VARIABLE = "Ocean Surface Currents" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.DATATYPE = "1/72 YEAR Interval" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.DATASUBTYPE = "unfiltered" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.GEORANGE = "20 to 420 -80 to 80" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.PERIOD = "Jun.07,2013" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.YEAR = "2013" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.DESCRIPTION = "OSCAR Third Degree Sea
Surface Velocity" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.CREATION_DATE = "00:33 09-Jun-2013" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.VERSION = 2009.f ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.SOURCE = "Gary Lagerloef, ESR
(<a class="moz-txt-link-abbreviated" href="mailto:lager@esr.org">lager@esr.org</a>) and Kathleen Dohan, ESR (<a class="moz-txt-link-abbreviated" href="mailto:kdohan@esr.org">kdohan@esr.org</a>)" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.CONTACT = "Kathleen Dohan
(<a class="moz-txt-link-abbreviated" href="mailto:kdohan@esr.org">kdohan@esr.org</a>) or John T. Gunn (<a class="moz-txt-link-abbreviated" href="mailto:gunn@esr.org">gunn@esr.org</a>)" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.COMPANY = "Earth & Space Research,
Seattle, WA" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.REFERENCE = "Bonjean F. and G.S.E.
Lagerloef, 2002 ,\\\"Diagnostic model and analysis of the surface
currents in the tropical Pacific ocean\\\", J. Phys. Oceanogr.,
32, 2,938-2,954" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.NOTE1 = "Maximum Mask velocity is the
geostrophic component at all points + any concurrent Ekman and
buoyancy components" ;</tt><tt><br>
</tt><tt> :NC_GLOBAL.NOTE2 = "Longitude extends from 20 E to
420 E to avoid a break in major ocean basins. Data repeats in
overlap region." ;</tt><tt><br>
</tt><tt>}</tt><tt><br>
</tt><br>
<br>
Here is the complete GrADS output from lats4d.sh:<br>
<br>
<br>
<tt>/opt/opengrads/grads -blc 'run lats4d.gs -q -v -i
../../sources/oscar/oscar_vel7550.nc.gz.nc.ctl -lat 16.8 31.8 -lon
276.0 291.0 -format grib -table
../../sources/oscar/grads.lats.table -ntimes 1 -o
/tmp/tmpCHWp85.grb '</tt><tt><br>
</tt><tt><br>
</tt><tt> Welcome to the OpenGrADS Bundle Distribution</tt><tt><br>
</tt><tt> --------------------------------------------</tt><tt><br>
</tt><tt><br>
</tt><tt>For additional information enter "grads -h".</tt><tt><br>
</tt><tt><br>
</tt><tt>Starting
"/opt/opengrads/Linux/Versions/2.0.1.oga.1/x86_64/grads -blc
'run lats4d.gs -q -v -i
../../sources/oscar/oscar_vel7550.nc.gz.nc.ctl -lat 16.8 31.8 -lon
276.0 291.0 -format grib -table
../../sources/oscar/grads.lats.table -ntimes 1 -o
/tmp/tmpCHWp85.grb ' " ...</tt><tt><br>
</tt><tt><br>
</tt><tt><br>
</tt><tt>Grid Analysis and Display System (GrADS) Version
2.0.1.oga.1</tt><tt><br>
</tt><tt>Copyright (c) 1988-2011 by Brian Doty and the</tt><tt><br>
</tt><tt>Institute for Global Environment and Society (IGES)</tt><tt><br>
</tt><tt>GrADS comes with ABSOLUTELY NO WARRANTY</tt><tt><br>
</tt><tt>See file COPYRIGHT for more information</tt><tt><br>
</tt><tt><br>
</tt><tt>Config: v2.0.1.oga.1 little-endian readline printim grib2
netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile</tt><tt><br>
</tt><tt>Issue 'q config' command for more detailed configuration
information</tt><tt><br>
</tt><tt>Loading User Defined Extensions table
</opt/opengrads/Linux/Versions/2.0.1.oga.1/x86_64/gex/udxt>
... ok.</tt><tt><br>
</tt><tt>GX Package Initialization: Size = 11 8.5 </tt><tt><br>
</tt><tt>Running in Batch mode</tt><tt><br>
</tt><tt>lats4d: Version 2.1.5 of 28 June 2010</tt><tt><br>
</tt><tt>No hardcopy metafile open</tt><tt><br>
</tt><tt>All files closed; all defined objects released;</tt><tt><br>
</tt><tt>All GrADS attributes have been reinitialized</tt><tt><br>
</tt><tt>lats4d: Data file is </tt><tt><br>
</tt><tt>File 1 : </tt><tt><br>
</tt><tt> Descriptor:
../../sources/oscar/oscar_vel7550.nc.gz.nc.ctl</tt><tt><br>
</tt><tt> Binary: ../../sources/oscar/oscar_vel7550.nc.gz.nc</tt><tt><br>
</tt><tt> Type = Gridded</tt><tt><br>
</tt><tt> Xsize = 1201 Ysize = 481 Zsize = 1 Tsize = 1 Esize =
1</tt><tt><br>
</tt><tt> Number of Variables = 2</tt><tt><br>
</tt><tt> u 0 t,z,x,y Ocean Surface Zonal Currents [m/s]</tt><tt><br>
</tt><tt> v 0 t,z,x,y Ocean Surface Meridional Currents [m/s]</tt><tt><br>
</tt><tt><br>
</tt><tt>lats4d: Dimension environment file same as data file</tt><tt><br>
</tt><tt>lats4d: time range: 00Z07JUN2013 00Z07JUN2013 by 1, delta
t: 24 hourly</tt><tt><br>
</tt><tt>lats4d: vertical levels: 0</tt><tt><br>
</tt><tt>lats4d: invalid plev 0 for GRIB output</tt><tt><br>
</tt><tt>lats4d: using HYBRID level number for vertical coordinate</tt><tt><br>
</tt><tt>lats4d: latitudinal range: 16.8 31.8</tt><tt><br>
</tt><tt>lats4d: longitudinal range: 276.0 291.0</tt><tt><br>
</tt><tt>gancgrid error: nc_get_vara_double failed; NetCDF: Index
exceeds dimension bound</tt><tt><br>
</tt><tt>Data Request Error: Error for variable 'u'</tt><tt><br>
</tt><tt> Error ocurred at column 1</tt><tt><br>
</tt><tt>lats4d: error creating grib_only file
/tmp/tmpCHWp85.grb.grb</tt><tt><br>
</tt><tt><RC> 1 </RC></tt><tt><br>
</tt><tt></IPC></tt><tt><br>
</tt><tt>lats4d: exiting from GrADS...</tt><tt><br>
</tt><tt>No hardcopy metafile open</tt><tt><br>
</tt><tt>GX package terminated </tt><tt><br>
</tt><br>
<br>
<br>
<br>
Thank You,<br>
Derek Simkowiak<br>
<br>
</body>
</html>