SPI calculation

Dear all,

I run the SPI codes with another data set,but still get errors!!!, below
are the detailes of the data set , the code I have run and the error I got.
Please help me to sort it out.

Kind regards,
______________________________________________________________________________________________________________________________________
[rahimi@localhost seti]$ ncdump -h precip.mon.total.nc
netcdf precip.mon.total {
dimensions:
    lat = 360 ;
    lon = 720 ;
    nbnds = 2 ;
    time = UNLIMITED ; // (1320 currently)
variables:
    float lat(lat) ;
        lat:long_name = "Latitude" ;
        lat:units = "degrees_north" ;
        lat:standard_name = "latitude" ;
        lat:axis = "Y" ;
        lat:coordinate_defines = "center" ;
        lat:actual_range = 89.75f, -89.75f ;
    float lon(lon) ;
        lon:long_name = "Longitude" ;
        lon:units = "degrees_east" ;
        lon:standard_name = "longitude" ;
        lon:actual_range = 0.25f, 359.75f ;
        lon:axis = "X" ;
        lon:coordinate_defines = "center" ;
    double time(time) ;
        time:long_name = "Time" ;
        time:units = "days since 1800-1-1 00:00:00" ;
        time:actual_range = 36889., 77035. ;
        time:delta_t = "0000-01-00 00:00:00" ;
        time:standard_name = "time" ;
        time:axis = "T" ;
        time:coordinate_defines = "start" ;
        time:totalling_period = "0000-01-00 00:00:00" ;
    float precip(time, lat, lon) ;
        precip:add_offset = 0.f ;
        precip:scale_factor = 1.f ;
        precip:missing_value = -9.96921e+36f ;
        precip:units = "mm" ;
        precip:actual_range = 0.f, 4412.31f ;
        precip:dataset = "GPCC Precipitation Full V6 (0.5x0.5)" ;
        precip:var_desc = "Precipitation" ;
        precip:level_desc = "Surface" ;
        precip:statistic = "Total" ;
        precip:parent_stat = "Observations" ;
        precip:long_name = "GPCC Monthly total of precipitation" ;
        precip:valid_range = 0.f, 8000.f ;

// global attributes:
        :history = "Created 07/2012 based on V5 data obtained via ftp" ;
        :title = "GPCC Full Data Reanalysis Version 6 0.5x0.5 Monthly
Means" ;
        :convention = "CF 1.0" ;
        :Documentation = "
http://www.esrl.noaa.gov/psd/data/gridded/data.gpcc.html" ;
        :Original_Source = "
http://www.dwd.de/en/FundE/Klima/KLIS/int/GPCC/GPCC.htm\n",
            " is the webpage and the data is at
ftp://ftp.dwd.de/pub/data/gpcc/download.html" ;
        :Reference = "Users of the data sets are kindly requested to give
feed back and to refer to GPCC publications on this webpage:
http://www.dwd.de/bvbw/appmanager/bvbw/dwdwwwDesktop/?_nfpb=true&_pageLabel=_dwdwww_klima_umwelt_datenzentren_wzn&T12404518261141645246564gsbDocumentPath=Content%2FOeffentlichkeit%2FKU%2FKU4%2FKU42%2Fteaser__product__access.html&_state=maximized&_windowLabel=T12404518261141645246564&lastPageLabel=_dwdwww_klima_umwelt_datenzentren_wzn"
;
        :original_source = "
ftp://ftp-anon.dwd.de/pub/data/gpcc/html/fulldata_download.html" ;
______________________________________________________________________________________________________________________________________
load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

   f = addfile("precip.mon.total.nc","r")
   prc = f->precip
   dimp = dimsizes(prc)
   ntim = dimp(0)
   nlat = dimp(1)
   mlon = dimp(2)

   len = (/3, 6, 12, 24, 36, 48 /)
   klen = dimsizes(len)
   spi = new((/klen,ntim,nlat,mlon/) ,float,prc@_FillValue)

   do k=0,klen-1
      spi(k,:,:,:) = dim_spi_n(prc, len(k), False, 0)
   end do

   print(sprintf("%8.2f", prc) \
               +sprintf("%8.2f", spi(0,:,:,:)) \
               +sprintf("%8.2f", spi(1,:,:,:)) \
               +sprintf("%8.2f", spi(2,:,:,:)) \
               +sprintf("%8.2f", spi(3,:,:,:)) \
               +sprintf("%8.2f", spi(4,:,:,:)) \
               +sprintf("%8.2f", spi(5,:,:,:)) )
end
__________________________________________________________________________________________________________
fatal:NclMalloc Failed:[errno=12]
fatal:New: could not create new array:[errno=12]
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 16 in file
spi2.ncl

[rahimi@localhost seti]$
---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

-- 
S.Rahimi

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jul 29 16:14:28 2013