<div dir="ltr"><div><div><div>10,000 => 3 sec; 20,000 => 10 sec; 25,000 => 15 sec<br><br></div></div>To my knowledge, there is no 'internal' NCL reason for 29,737 => "infinity" sec<br><br>===<br></div><div>To my knowledge, all TRMM data are float ... not double.<br></div><div>Hence, your comment "I switched from double to floating point" is a surprise.<br></div><div><br></div><div>FYI: <br>*All* NCL computational functions use double precision. Float values are <font>temporarily </font>promoted to double for computations. <br></div><div>The results are then 'demoted' to float and returned. Punch line: users need not create double versions of the data.<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 4, 2015 at 12:56 PM, David Adams <span dir="ltr"><<a href="mailto:dave.k.adams@gmail.com" target="_blank">dave.k.adams@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi all, <br></div>I am using the trend_manken<br></div>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.<br><br>;---------------------------------<br>; Mann Kendall test<br>;---------------------------------<br>pa = trend_manken(precip_ave, False, 0)<br>print(pa)<br><br><br><br></div>Any ideas?<br><br></div>saludos,<br></div>Dave<br></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>