<div dir="ltr"><div class="gmail_default" style="font-size:small">Please see example "climo_1.ncl" at:</div><div class="gmail_default" style="font-size:small"><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="">I admittedly don't know much about ttest, but this example uses it in a way that sounds similar to your request.</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, Nov 3, 2015 at 7:54 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><div>Hi all,<br><br></div>I would like to plot t-test figure, between the average of future and the average of the present periods. Both period consist of 2 continuous days as follow:<br><br></div><div>Future: 01 to 02 January 2100 <br></div><div>Present: 01 to 02 Janaury 2010<br><br></div>Could anyone help me on this? Your help is highly appreciated. <br><div><div><div><br><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;script;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br><br>begin<br> a1 = addfile("./JAN2100/01-05/wrfout_d02_2100-01-01_00:00:00","r")<br> a2 = addfile("./JAN2100/01-05/wrfout_d02_2100-01-02_00:00:00","r")<br><br> b1 = addfile("./JAN2010/01-05/wrfout_d02_2010-01-01_00:00:00","r")<br> b2 = addfile("./JAN2010/01-05/wrfout_d02_2010-01-02_00:00:00","r")<br><br> wks = gsn_open_wks("x11","T-test1-90%")<br> gsn_define_colormap(wks,"BlAqGrYeOrReVi200") ; Create a plot workstation<br><br> opts = True ; Set some Basic Plot options<br> opts@MainTitle = "GEOGRID FIELDS"<br> opts@InitTime = False ; Do not plot time or footers<br> opts@Footer = False<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;Get the Variables;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br>tempa1 = a1->T2<br>tempa2 = a2->T2<br><br>tempb1 = b1->T2<br>tempb2 = b2->T2<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;Calculate average and variance for Future;;;;;;;;;;;;;;;;;;;;<br><br>aveX1 = dim_avg_n_Wrap(tempa1,0) ; average over the 0th dim<br>varX1 = dim_variance_n_Wrap(tempa1,0) ; compute variance<br><br>aveX2 = dim_avg_n_Wrap(tempa2,0) ; average over the 0th dim<br>varX2 = dim_variance_n_Wrap(tempa2,0) ; compute variance<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;Calculate average and variance for Present;;;;;;;;;;;;;;;;;;;;<br><br>aveY1 = dim_avg_n_Wrap(tempb1,0) ; average over the 0th dim<br>varY1 = dim_variance_n_Wrap(tempb1,0) ; compute variance<br><br>aveY2 = dim_avg_n_Wrap(tempb2,0) ; average over the 0th dim<br>varY2 = dim_variance_n_Wrap(tempb2,0) ; compute variance<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;Do T-Test for 1st day and 2nd day;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br>alphat1 = 100.*(1. - ttest(aveX1,varX1,24, aveY1,varY1,24, True, False))<br>aveX1 = where(alphat1.lt.90.,aveX1@_FillValue, aveX1)<br><br>alphat2 = 100.*(1. - ttest(aveX2,varX2,24, aveY2,varY2,24, True, False))<br>aveX2 = where(alphat2.lt.90.,aveX2@_FillValue, aveX2)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;Do T-Test for the average of 1st and 2nd days;;;;;;;;;;;;;;;;;;<br><br>result = (alphat1 + alphat2)/2<br>;result = (aveX1 + aveX2)/2<br><br></div></div></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>