<html>
  <head>

    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    &nbsp;&nbsp;&nbsp; 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.&nbsp; 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.&nbsp; However, I am
    unable to get that to work.&nbsp; 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:&nbsp; Error for variable 'u'</tt><tt><br>
    </tt><tt>&nbsp; Error ocurred at column 1</tt><tt><br>
    </tt><br>
    I don't know what to try next.&nbsp; 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=&gt;u 0 t,z,x,y Ocean Surface Zonal Currents [m/s]</tt><tt><br>
    </tt><tt>v=&gt;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.&nbsp; 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.&nbsp; 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>&nbsp;&nbsp;&nbsp; time = 1 ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; year = 1 ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; depth = 1 ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; latitude = 481 ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; longitude = 1201 ;</tt><tt><br>
    </tt><tt>variables:</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; int time(time) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; time:units = "day since 1992-10-05 00:00:00" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; time:long_name = "Day since 1992-10-05 00:00:00" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; float year(year) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; year:units = "time in years" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; year:long_name = "Time in fractional year" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; float depth(depth) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; depth:units = "meter" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; depth:long_name = "Depth" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double latitude(latitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; latitude:units = "degrees-north" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; latitude:long_name = "Latitude" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double longitude(longitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; longitude:units = "degrees-east" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; longitude:long_name = "Longitude" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double u(time, depth, latitude, longitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; u:units = "meter/sec" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; u:long_name = "Ocean Surface Zonal Currents" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; u:missing_value = 0. ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double v(time, depth, latitude, longitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; v:units = "meter/sec" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; v:long_name = "Ocean Surface Meridional Currents" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; v:missing_value = 0. ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double um(time, depth, latitude, longitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; um:units = "meter/sec" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; um:long_name = "Ocean Surface Zonal Currents
      Maximum Mask" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; um:missing_value = 0. ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; double vm(time, depth, latitude, longitude) ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; vm:units = "meter/sec" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; vm:long_name = "Ocean Surface Meridional Currents
      Maximum Mask" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; vm:missing_value = 0. ;</tt><tt><br>
    </tt><tt><br>
    </tt><tt>// global attributes:</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.VARIABLE = "Ocean Surface Currents" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.DATATYPE = "1/72 YEAR Interval" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.DATASUBTYPE = "unfiltered" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.GEORANGE = "20 to 420 -80 to 80" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.PERIOD = "Jun.07,2013" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.YEAR = "2013" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.DESCRIPTION = "OSCAR Third Degree Sea
      Surface Velocity" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.CREATION_DATE = "00:33 09-Jun-2013" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.VERSION = 2009.f ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :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>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :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>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.COMPANY = "Earth &amp; Space Research,
      Seattle, WA" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :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>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :NC_GLOBAL.NOTE1 = "Maximum Mask velocity is the
      geostrophic component at all points + any concurrent Ekman and
      buoyancy components" ;</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; :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>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Welcome to the OpenGrADS Bundle Distribution</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; --------------------------------------------</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&nbsp;&nbsp; -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
      &lt;/opt/opengrads/Linux/Versions/2.0.1.oga.1/x86_64/gex/udxt&gt;
      ... 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>&nbsp; Descriptor:
      ../../sources/oscar/oscar_vel7550.nc.gz.nc.ctl</tt><tt><br>
    </tt><tt>&nbsp; Binary: ../../sources/oscar/oscar_vel7550.nc.gz.nc</tt><tt><br>
    </tt><tt>&nbsp; Type = Gridded</tt><tt><br>
    </tt><tt>&nbsp; Xsize = 1201&nbsp; Ysize = 481&nbsp; Zsize = 1&nbsp; Tsize = 1&nbsp; Esize =
      1</tt><tt><br>
    </tt><tt>&nbsp; Number of Variables = 2</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp;&nbsp; u&nbsp; 0&nbsp; t,z,x,y&nbsp; Ocean Surface Zonal Currents [m/s]</tt><tt><br>
    </tt><tt>&nbsp;&nbsp;&nbsp;&nbsp; v&nbsp; 0&nbsp; t,z,x,y&nbsp; 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:&nbsp; 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&nbsp; 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:&nbsp; Error for variable 'u'</tt><tt><br>
    </tt><tt>&nbsp; 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>&lt;RC&gt; 1 &lt;/RC&gt;</tt><tt><br>
    </tt><tt>&lt;/IPC&gt;</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>