[ncl-talk] trend_manken data size limit?

[Mary Haley]

Dave,

Could you write the precip_ave variable to a NetCDF file and send us the
file, so we can try trend_manken on it here?

You can do something like:

fout = addfile("test_trend_manken.nc","c")
fout->pcp = precip_ave

and then if the file is too large (say over 25 megabytes), then you can use
our ftp to upload it. See

http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP

Thanks,

--Mary


On Sat, Apr 4, 2015 at 12:56 PM, David Adams <dave.k.adams at gmail.com> wrote:

> 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
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150405/b57f40c3/attachment.html