[ncl-talk] do loop and rtest problem

[Hyun Marcia]

Hi, I'm trying to significant testing for correlation(rtest)

And I want to use my own degree of freedom.

;====  rtest ====================================
 ;==== dimension : sst_DJF_anom(time,lat,lon), indx_sst(time)
;=============================================
rr_DJF=sst_DJF_anom(lat|:,lon|:,time|:)


do LAG=0,dimsizes(indx_sst)-1
 crr_indx=escar(indx_sst,LAG)
 crr_DJF=escar(rr_DJF,LAG)
end do

c =0
do i=1,dimsizes(indx_sst)
 c= c+crr_indx*crr_DJF
end do

d = 1+ 2*c


df = dimsizes(indx_sst)/d

 r_test=rtest(rtmp, df, 0)

But, there are some errors occured in this script.

I think the do loop has ploblems but i can't find it.

Anyone who tell me what problems it is?

Thank.

Marcia Hyun
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140813/5c432105/attachment.html