[ncl-talk] trend_manken data size limit?

[David Adams]

Hi all,
I am using the trend_manken
function to look at the significance of  precipitation intensity trends
(using TRMM data) for 15 years (every 3 hours data).  When I use 10,000
data points it takes about 3 seconds to run, 20,000 data points about 10
seconds, 25,000 about 15 seconds.  When I use the entire time series 29,737
data points, it never finishes.  I have waited more than 20 minutes and it
keeps spinning.  I switched from double to floating point, but that makes
no difference.

;---------------------------------
; Mann Kendall test
;---------------------------------
pa = trend_manken(precip_ave, False, 0)
print(pa)



Any ideas?

saludos,
Dave
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150404/26a7b0a3/attachment.html