<div dir="ltr"><div>I was about to delete this from the queue but thought I'd add something.<br><br>FYI only: <br><br></div><div>If the residuals of the line derived via simple linear regression (SLR) are normally distributed, <br></div><div>there is virtually no difference between the regression and the M-K (Thiel-Sen) estimates.<br><br></div><div>See attached which shows the SLR and M-K (Thiel-Sen) fits for the 1st half; 2nd half and whole series.<br><br></div><div>Again: the NCL Mann-Kendall (Thiel-Sen) function will be available in 6.3.0<br></div><div><br></div>time series of global anomalies: <a href="http://hadobs.metoffice.com/hadcrut3/diagnostics/global/nh+sh/">http://hadobs.metoffice.com/hadcrut3/diagnostics/global/nh+sh/</a><br><br>The data are annual means and reside in column 2 of the input file. The file structure is described at:<br><a href="http://hadobs.metoffice.com/hadcrut3/diagnostics/time-series.html">http://hadobs.metoffice.com/hadcrut3/diagnostics/time-series.html</a><br><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 5, 2014 at 2:41 PM, Xi Chang <span dir="ltr"><<a href="mailto:xi.chang01@gmail.com" target="_blank">xi.chang01@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">sorry i missed the fortran code.<div>i found the significant test in the original <span></span>code.</div><div><br></div><div>Thanks</div><div class="HOEnZb"><div class="h5"><div><br></div><div><br><br>On Friday, September 5, 2014, Xi Chang <<a href="mailto:xi.chang01@gmail.com" target="_blank">xi.chang01@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">and how to test whether the trend is significant at 95% or not?<div><br></div><div>Thanks</div><div><br></div><div><span></span><br><br>On Friday, September 5, 2014, Ruping Mo <<a>rupingmo@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thank you very much, Mary. I will give this a try. I also want to thank Dennis for his effort.<div><br></div><div>Ruping<br><br>On Friday, September 5, 2014, Mary Haley <<a>haley@ucar.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Ruping,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">On occasion, we will create documentation for upcoming functions, even though they are not available yet. We always try to put the version number on new functions, so you know if your version of NCL has that function or not.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">We have not actually implemented the wrapper for this function yet, but Dennis had provided the Fortran code, a sample NCL script, and a sample dataset.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you want to give this a try, then download the three attached files, then type:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><font face="courier new, monospace">WRAPIT kens_trend_dp.ncl.f<br></font></div><div class="gmail_default"><font face="courier new, monospace">ncl test.ken_trend.ncl</font><br></div><div class="gmail_default"><br></div><div class="gmail_default">You should see output like this:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="courier new, monospace">(0)<span style="white-space:pre-wrap"> </span>================ TA ==================================</font></div><div class="gmail_default"><font face="courier new, monospace"> n=165 s= 8059 nc= 13530 nt=165 var=503597.0 vartie= 19.7 trend= 0.00436</font></div><div class="gmail_default"><font face="courier new, monospace">(0)<span style="white-space:pre-wrap"> </span> stau=8059</font></div><div class="gmail_default"><font face="courier new, monospace">(0)<span style="white-space:pre-wrap"> </span> z=11.35496228080661</font></div><div class="gmail_default"><font face="courier new, monospace">(0)<span style="white-space:pre-wrap"> </span>prob =1</font></div><div class="gmail_default"><font face="courier new, monospace">(0)<span style="white-space:pre-wrap"> </span>trend=0.004361904761904762</font></div><div class="gmail_default"><br></div><div class="gmail_default">If you have questions about this, please post back to ncl-talk.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Sep 4, 2014 at 2:06 PM, Ruping Mo <span dir="ltr"><<a>rupingmo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear NCL,<div><br></div><div>I found this new function in the function list. Is it available now, or in version 6.3 only?</div><div><br></div>Regards,<div><br></div><div>Ruping</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>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</div></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>