Re: [ferret_users] Monthly climatologies repeated along long time axis

[William Kessler <william.s.kessler@xxxxxxxx>]

The only thing I would add would be to overplot the pieces to make sure everything is right. Especially ACL and ACLMOD. The modulo operation can reduce the amplitude of the annual cycle if, for example, ACL values are on the first of the month (as in some data sets, e.g. ERA-Interim) but month_reg is centered on the middle of the month. month_reg is linearly interpolating between values which can reduce peaks. 

Similarly, you might get slightly different values from month_reg vs month_irreg for values on the first of the month. It is possible to define an analog for month_irreg with values on the first day of each month in that case. I could dig up an example.

Whenever I regrid as here, I always make a few test plots of the intermediate steps to check. This is so easy to do in Ferret, and is one reason to do the calculation in explicit steps as in my example.

Billy

> On Jul 10, 2017, at 9:16 AM, Paul Goddard <pgoddard@xxxxxxxxxxxxxxxxx> wrote:
> 
> Perfect, thanks Billy
> 
> Best,
> Paul
> 
> On Sun, Jul 9, 2017 at 6:22 PM, William S. Kessler <william.s.kessler@xxxxxxxx> wrote:
> If A is your 20-year time series and ACL is your 12-month climatology:
> 
> First, put the climatology on a modulo time axis:
> 
> let ACLMOD = ACL[gt=month_reg@mod]
> 
> Extend it to the time axis of A:
> 
> let ACLRG = ACLMOD[gt=A]
> 
> let ADM = A - ACLRG  ! anomalies
> 
> Sent from my phone .... Please excuse my brevity .... and my typing.
> 
> > On Jul 9, 2017, at 4:06 PM, Paul Goddard <pgoddard@xxxxxxxxxxxxxxxxx> wrote:
> >
> > Hi all,
> >
> > It is likely that an answer is already on the ferret users archive, if so, just refer me to that page, but I could not find what I was looking for.
> >
> > I have a monthly climatology, so 12 time steps.  I would like to subtract the climatology off a 20 year monthly data set (of 240 time steps) to end with an anomaly with 240 time steps. How can I regrid the 12 climatologies to a 240 time step axis with the 12 months of the climatology just repeated 20 times.
> >
> > Thanks in advance,
> > Paul
>