<div dir="ltr"><div class="gmail_default" style="font-size:small">Did you look at our examples page?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style=""><a href="http://www.ncl.ucar.edu/Applications/">http://www.ncl.ucar.edu/Applications/</a><br></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">There are hundreds of examples showing you how to create different types of graphics.</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">In particular, see example climo_1.ncl at:</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><a href="http://www.ncl.ucar.edu/Applications/climo.shtml">http://www.ncl.ucar.edu/Applications/climo.shtml</a><br></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">--Mary</div><div class="gmail_default" style=""><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 20, 2015 at 8:52 PM, Steven Kong <span dir="ltr"><<a href="mailto:kongksk@gmail.com" target="_blank">kongksk@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>hi all, <br><br></div><div>Current i would like to a plot a t-test output, a significant change of daily mean surface temperature between 01 Janaury and 02 January. Could any one help me?<br><br>;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 01 January;;;;;;;;;;;;;;<br><br>temp1=a->T2(0,:,:)<br>temp1=0<br> do i = 0,23<br> temp1 = temp1 + a->T2(i,:,:)<br> end do<br><br> mean1 = (temp1)/24<br><br>;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 02 January;;;;;;;;;;;;;;<br><br>temp2=b->T2(0,:,:)<br>temp2=0<br> do i = 0,23<br> temp2 = temp2 + b->T2(i,:,:)<br> end do<br><br> mean2 = (temp2)/24<br><br>;-------------T-Test: Area with 95% significant<br>difference----------------------------------------------<br><br>r_ave1 = avg(mean1)<br>r_var1 = variance(mean1)<br><br>r_ave2 = avg(mean2)<br>r_var2 = variance(mean1)<br><br>probt = 100.*(1. -ttest(r_ave1,r_var1,24,r_ave2,r_var2,24, True, False)) ; all areas whose value > 95 are significant<br>aveX = where(probt.lt.95.,r_ave1@_FillValue, r_ave1)<br><br><br> <br></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>