[Cdat-discussion] problem reading reanalysis file

[Jeff W. <js...@fa...>]

Here's a simple script to read reanalsis precip:

import cdms
filenameverif = 'prate.sfc.gauss.2002.nc'
fileverif = cdms.open(filenameverif,mode='r')
datreanl = fileverif['prate']
print datreanl.shape
for n in range(datreanl.rank()):
    print datreanl.getAxis(n)
testdata = cdms.MV.filled(datreanl[0,:,:])

(the file prate.sfc.gauss.2002.nc is available at
ftp://ftp.cdc.noaa.gov/Datasets/ncep.reanalysis/surface_gauss)

When I run it with cdat-3.3, python 2.3 I get this error:

/sw/lib/python2.3/xmllib.py:10: DeprecationWarning: The xmllib module is obsolete.  Use xml.sax instead.
  DeprecationWarning)
(1460, 94, 192)
   id: time
   Designated a time axis.
   units:  hours since 1-1-1 00:00:0.0
   Length: 1460
   First:  17540448.0
   Last:   17549202.0
   Other axis attributes:
      avg_period: 0000-00-00 06:00:00
      long_name: Time
      delta_t: 0000-00-00 06:00:00
      actual_range: [ 17540448., 17549202.,]
   Python id:  9ca170

   id: lat
   Designated a latitude axis.
   units:  degrees_north
   Length: 94
   First:  88.5419998169
   Last:   -88.5419998169
   Other axis attributes:
      long_name: Latitude
      actual_range: [ 88.54199982,-88.54199982,]
   Python id:  9c5af8

   id: lon
   Designated a longitude axis.
   units:  degrees_east
   Length: 192
   First:  0.0
   Last:   358.125
   Other axis attributes:
      long_name: Longitude
      actual_range: [   0.   , 358.125,]
   Python id:  9ca328

Traceback (most recent call last):
  File "testcdms.py", line 8, in ?
    testdata = cdms.MV.filled(datreanl[0,:,:])
  File "/sw/lib/python2.3/site-packages/cdms/variable.py", line 78, in __getitem__
    return AbstractVariable.__getitem__(self, key)
  File "/sw/lib/python2.3/site-packages/cdms/avariable.py", line 967, in __getitem__
    return self.getSlice(numericSqueeze=1, raw=0, *speclist)
  File "/sw/lib/python2.3/site-packages/cdms/avariable.py", line 398, in getSlice
    return self.subSlice(*specs, **keys)
  File "/sw/lib/python2.3/site-packages/cdms/avariable.py", line 351, in subSlice
    id = self.id)
  File "/sw/lib/python2.3/site-packages/cdms/tvariable.py", line 41, in __init__
    self.initDomain(axes)
  File "/sw/lib/python2.3/site-packages/cdms/tvariable.py", line 96, in initDomain
    raise CDMSError, "Wrong number of axes to initialize domain."
cdms.error.CDMSError: Wrong number of axes to initialize domain.

I swear this used to work, but the file was re-created recently after a
disk failure.  Is a problem with my script, the file, or cdms?

-Jeff

-- 
Jeffrey S. Whitaker         Phone : (303)497-6313
NOAA/OAR/CDC  R/CDC1        FAX   : (303)497-6449
325 Broadway                Web   : http://www.cdc.noaa.gov/~jsw
Boulder, CO, USA 80305-3328 Office: Skaggs Research Cntr 1D-124