The codes I run are below.
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.1981-2010.ltm.v6(1).nc","r")
prc2 = f->precip(:,{-35:37},{344:50})
yyyymm = f->time
printVarSummary(yyyymm)
ntime = dimsizes(yyyymm)
year = yyyymm/100
yrs= year(0)
yrl = year (ntime-1)
yyyymm = yyyymm_time(yrs,yrl,"double")
yrfrac = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)
printVarSummary(prc2)
dimp = dimsizes(prc2)
ntim = dimp(0)
nlat = dimp(1)
mlon = dimp(2)
len = (/3,6,12,24/)
klen = dimsizes(len)
spi = new((/klen,ntim,nlat,mlon/),float)
do k=0,klen-1
spi(k,:,:,:) = dim_spi_n(prc2, len(k), False, 0)
end do
print(yyyymm+sprintf("%8.2f",spi(0,:,:,:))\
+sprintf("%8.2f",spi(1,:,:,:))\
+sprintf("%8.2f",spi(2,:,:,:))\
+sprintf("%8.2f",spi(3,:,:,:)))
;wks = gsn_open_wks ("ps","spi")
;res = True
;res@gsnDraw = False
;res@gsnFrame = False
;res@vpHeightF = 0.4
;es@vpWidthF = 0.8
;res@vpXF = 0.1
;res@trYMinF = -3.0
;res@trYMaxF = 3.0
;res@gsnYRefLine = 0.0
;res@xyMonoDashPattern = True
;res@xyLineThicknessF = 1
;plt = new ( klen, "graphic")
;xyLineColors = (/"black","red","blue","green","brown","magenta"/)
;do n=0,klen-1
; res@xyLineColors = xyLineColors(n)
; res@gsnCenterString = "run="+len(n)
;plt(n) = gsn_csm_xy (wks,yrfrac,spi(n,:,:,:),res)
;end do
;resP = True ; modify the panel plot
;resP@txString = "SPI: Boulder: 1894-2010"
;resP@gsnMaximize = True
;gsn_panel(wks,plt,(/2,2/),resP)
and this are the errors
Variable: yyyymm
Type: double
Total Size: 96 bytes
12 values
Number of Dimensions: 1
Dimensions and sizes: [time | 12]
Coordinates:
time: [-657073..-656739]
Number Of Attributes: 14
long_name : Time
units : days since 1800-1-1 00:00:00
actual_range : ( -657073, -656739 )
delta_t : 0000-01-00 00:00:00
standard_name : time
axis : T
coordinate_defines : start
totalling_period : 0000-01-00 00:00:00
avg_period : 0030-00-00 00:00:00
prev_avg_period : 0000-01-00 00:00:00
climatology : climatology_bounds
climo_period : 1981/01/01 - 2010/12/31
ltm_range : ( 66109, 77035 )
interpreted_actual_range : 0001/01/01 00:00:00 - 0001/12/01 00:00:00
fatal:Argument type mismatch on argument (0) of (yyyymm_time) can not coerce
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 15 in file
spi2se
end
---------- Forwarded message ----------
From: Mary Haley <haley@ucar.edu>
Date: Wed, Jul 31, 2013 at 8:42 PM
Subject: Re: time
To: Setareh Rahimi <setareh.rahimi@gmail.com>
I'm taking this offline.
What exactly does your new code look like? Please send me the whole script
and not just a subset.
What do you mean "type of time is double"? I never saw a "time" array.
What is the exact error message?
--Mary
On Jul 31, 2013, at 12:38 PM, Setareh Rahimi wrote:
> Type of time is double, I repllaced the type of time with double, but got
the same error again!!!
>
>
> --
> S.Rahimi
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- S.Rahimi
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 31 14:02:46 2013
This archive was generated by hypermail 2.1.8 : Thu Aug 01 2013 - 15:55:04 MDT