how to deal with an annual mean data

Hi,
Still the data from 'http://ecospriv4.epfl.ch/pub/KK10_2k_1850_30m.nc.gz'
(it's too large to be attached). It's an annual mean data, with variable
'land_use(time,lat,lon)'. I want to calculate the time average of the
variable 'land_use' or calculate the EOFs.
I tried 'dim_avg_n_Wrap' and 'run_ave_n_Wrap', they all failed.

ncl 0> f=addfile("KK10_2k_1850_30m.nc","r")
ncl 1> lu=f->land_use
ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

ncl 3>
ncl 4> lua=runave_n_Wrap(lu,10,0,0)
fatal:runave_n: Unable to allocate memory for output array
fatal:Execute: Error occurred at or near line 6009 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:Execute: Error occurred at or near line 4

ncl 5> luavftime=dim_avg_n_Wrap(lu,0)
fatal:NclMalloc Failed:[errno=12]
fatal:dim_avg_n: Could not coerce input data to double, can't continue
fatal:Execute: Error occurred at or near line 2041 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:Execute: Error occurred at or near line 5

And here is time information:

ncl 6> tim=f->time
ncl 7> printVarSummary(tim)

Variable: tim
Type: integer
Total Size: 7604 bytes
            1901 values
Number of Dimensions: 1
Dimensions and sizes: [time | 1901]
Coordinates:
            time: [2000..100]
Number Of Attributes: 3
  long_name : time
  units : years before 1950
  calendar : standard

I tried a monthly mean data, both functions worked. I don't know why it
failed for the annual mean data. What's wrong with my 'contributed.ncl'?
Thanks a lot!
Best regards,
Mi

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Feb 17 19:36:02 2011