<div dir="ltr"><div>'prob' is the probability<br><br></div>Well, what is 'prob1' ?<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 26, 2015 at 2:58 PM, Sunmin Park <span dir="ltr"><<a href="mailto:mireiyue@gmail.com" target="_blank">mireiyue@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 style="word-wrap:break-word">Dear NCL users<div><br></div><div>I am trying to get rtest and ttest values from cross correlation lag 0. ESCORC was used to calculate cross correlation and I did rtest first then ttest using the equation from here <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml</a>. However, the ttest value are not correct when I printout Min Max values (min=0 max=805.891).</div><div><br></div><div>How I can fix this? </div><div><br></div><div>Here is script I used to calculate.</div><div><br></div><div>massvsext1 = escorc(CMF3,EXT) ; calc cross correlations</div><div>massvsext = mask(massvsext1,landmask,0)</div><div>printVarSummary(massvsext)</div><div><br></div><div>Nr = dimsizes(CMF3&time)</div><div><br></div><div>prob = rtest(massvsext,Nr,0)</div><div>probt = prob1*sqrt((Nr-2)/(1-prob1^2))</div><div><br></div><div>probt!0 = "lat"</div><div>probt!1 = "lon"</div><div>lat = fspan(-90,90,180)</div><div>lon = fspan(-180,180,360)</div><div>probt&lat = lat</div><div>probt&lon = lon</div><div>lat@longname = "latitude"</div><div>probt&lat@units = "degrees_north"</div><div>lon@longname ="longitude"</div><div>probt&lon@units = "degrees_east"</div><div><br></div><div>printVarSummary(probt)</div><div>printMinMax(probt,0)</div><div><br></div><div>Thank you for helps and comments.</div><div>Sun-</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>