ttest warning message and code efficiency

Hi there,

I am using the ttest function on a global data. My arrays are of shape (12,
64, 128, 40) and correspond to month, lat, lon, time. I am currently running
this in a nested do loop which is not very efficient as I have cases where a
warning is written to the screen:

warning:ttest: encountered 1 cases where var1 and/or var2 were less than or
equal to 0. Output set to missing values in these cases

Would someone have a suggestion of how to do this outside of a do loop that
would improve the runtime of the script? My code looks as follows:

do mon=0,11
 do i=0,63
  do j=0,127
  aveX = dim_avg(FSNS_hi_new2(mon,i,j,:))
  aveY = dim_avg(FSNS_lo_new2(mon,i,j,:))
  varX = dim_variance(FSNS_hi_new2(mon,i,j,:))
  varY = dim_variance(FSNS_lo_new2(mon,i,j,:))
  sX = dimsizes(FSNS_hi_new2(mon,i,j,:))
  sY = dimsizes(FSNS_lo_new2(mon,i,j,:))
  iflag = True
  prob = 100.*(1. - ttest(aveX,varX,sX,aveY,varY,sY,iflag,False))
; printVarSummary(prob)
; print(prob)
  FSNSprob(mon,i,j) = prob
  end do
 end do
end do

Any suggestions would be greatly appreciated.
Thanks, Gina.

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 23 2009 - 07:39:06 MDT