xClm(12,lev,lat,lon), x(time,lev,lat,lon)
I have no idea what units tour time variable is so I'll assume yyyymm[*]
mm = yyyymm - (yyyymm/100)*100 ; mm[*] assuming 1 to 12
do nmo=1,12
imo = ind(mm.eq.nmo)
x(imo,:,:,:) = x(imo,:,:,:) - xClm(nmo,:,:,:)
delete(imo) ; may change on next iteration
end do
Good luck
Song Feng wrote:
> Suppose I have only one year daily data (e.g. [time, lev, lat, lon],
> and Ntime=365 or 366), is there a simple way to remove the monthly
> average? For example, I want to compute the daily temperature
> anomalies in January (Feb.....) by remove the mean temperature in
> January (Feb...). NCL did have some functions such as "calcDayAnomTLL"
> to compute daily anomalies but they are not applicable in our case.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jan 09 2009 - 10:13:11 MST
This archive was generated by hypermail 2.2.0 : Mon Jan 12 2009 - 16:21:33 MST