Re: Alternative for rmMonAnnCycTLL?

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Thu Feb 03 2011 - 16:12:38 MST

Hi Soumik,
There are various other functions that do the same thing. Take a look at
this list:
http://www.ncl.ucar.edu/Document/Functions/climo.shtml

There are various forms of rmMonAnnCycTLL (ex. rmAnnCyc1e1D), along with
functions that calculate the climatology (ex. clmMonLLT), and other
functions that take the climatology and remove it from an array (ex.
calcMonAnomLLT).

Now, of course, you could remove the annual cycle manually with a do
loop like this:

; arr is an array dimensioned time x lat x lon, starts in Jan,
; ends in Dec. (UNTESTED)

do gg = 0,11
    arrT = arr(gg::12,:,:)
    mon_avg = dim_avg_n_Wrap(arrT,0)
    printVarSummary(mon_avg)
    mon_avg_c = conform(arrT,mon_avg,(/1,2/))
    printVarSummary(mon_avg_c)
    delete(mon_avg)
    arr(gg::12,:,:) = (/ arrT - mon_avg_c /)
    delete(mon_avg_c)
    delete(arrT)
end do

Good luck,
Adam

On 02/03/2011 03:12 PM, sbasu@iarc.uaf.edu wrote:
> Hi,
>
> Is there any other way to calculate monthly anomalies from monthly SST
> data other than using "rmMonAnnCycTLL" ?
>
> Thanks,
> Soumik
>
>

-- 
__________________________________________________
Adam Phillips 
asphilli@ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
Climate and Global Dynamics Division         fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Feb 3 16:12:44 2011

This archive was generated by hypermail 2.1.8 : Fri Feb 04 2011 - 09:26:53 MST