<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-&gt;T2(0,:,:)<br>temp1=0<br>  do i = 0,23<br>    temp1 = temp1 + a-&gt;T2(i,:,:)<br>  end do<br><br>  mean1 = (temp1)/24<br><br>;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 02 January;;;;;;;;;;;;;;<br><br>temp2=b-&gt;T2(0,:,:)<br>temp2=0<br>  do i = 0,23<br>    temp2 = temp2 + b-&gt;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 &gt; 95 are significant<br>aveX = where(probt.lt.95.,r_ave1@_FillValue, r_ave1)<br><br><br> <br></div></div>