aura he5 data exceeds dimension

[Jennifer Adams]

Hi, John --
There is an extra dimension for those variables that the sdfopen code  
doesn't find because it's not 'coards compliant'. Check the ncdump  
output and you'll find:

dimensions:
         nWavelDiagnostic = 5 ;
         lat = 720 ;
         lon = 1440 ;
data:
  nWavelDiagnostic = 0, 1, 2, 3, 4 ;

So you have to tell GrADS how to map that "nWavelDiagnostic"  
dimension. I'd use EDEF for that. Plus, the variable names are very  
long, so you have to alias them to something more meaningful to you.  
Once you open the file, check the 'query attr' output for the titles  
of these variables -- they are also very long. Anyway, here's a bare- 
bones descriptor file that works (note the dset entry shouldn't wrap):

dset http://acdisc.sci.gsfc.nasa.gov/opendap/Aura_OMI_Level3/OMAEROe.003/2006/OMI-Aura_L3-OMAEROe_2006m0101_v003-2008m0815t092235.he5
dtype netcdf
undef -88888 missing_value
unpack scale_factor add_offset
xdef 1440 linear -179.875 0.25
ydef  720 linear -89.875 0.25
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN0001 1mn
edef 5 names  0 1 2 3 4
vars 13
AbsorbingAerosolOpticalThicknessMW=>var1          0  e,y,x  var1
AerosolModelMW=>var2                              0  y,x    var2
AerosolOpticalThicknessMW=>var3                   0  e,y,x  var3
AerosolOpticalThicknessPassedThresholdMean=>var4  0  y,x    var4
AerosolOpticalThicknessPassedThresholdStd=>var5   0  y,x    var5
SingleScatteringAlbedoMW=>var6                    0  e,y,x  var6
SingleScatteringAlbedoPassedThresholdMean=>var7   0  y,x    var7
SingleScatteringAlbedoPassedThresholdStd=>var8    0  y,x    var8
SolarZenithAngle=>var9      0  y,x  solarzenithangle
TerrainReflectivity=>var10  0  y,x  TerrainReflectivity
UVAerosolIndex=>var11       0  y,x  UVAerosolIndex
VISAerosolIndex=>var12      0  y,x  VISAerosolIndex
ViewingZenithAngle=>var13   0  y,x  ViewingZenithAngle
endvars


Jennifer

On Jan 8, 2009, at 10:20 AM, John Huddleston wrote:

> Hi
>
> I am accessing http://acdisc.sci.gsfc.nasa.gov/opendap/Aura_OMI_Level3/OMAEROe.003/2006/OMI-Aura_L3-OMAEROe_2006m0101_v003-2008m0815t092235.he5 
>  with the latest v2.0a4 gradsdap.
>
> It can open the file and the ‘q file’ works fine; however, when I  
> display some of the aerosol optical variables, there is an error:
>
> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds  
> dimension bound
>
> Data Request Error:  Error for variable 'aerosolopticalt'
>   Error ocurred at column 1
> DISPLAY error:  Invalid expression
>   Expression = aerosolopticalt
>
> Is this a problem with the data or bounds limitation in the gradsdap  
> application?
>
> Here is the information on the gradsdap startup
>
> Starting X server under c:\GrADS20/Xming
> Starting gradsdap under c:\GrADS20\win32 ...
>
> Grid Analysis and Display System (GrADS) Version 2.0.a4.oga.1
> 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.a3.oga.1 little-endian readline printim grib2 netcdf  
> hdf4-sds opendap-grids,stn Athena
>
> John Huddleston
>
>
>
>
>
>
>

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/45c7145b/attachment.html