Re: compute daily anomalies by remove the monthly average

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri, 09 Jan 2009 10:13:11 -0700

     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:
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