Re: climatology question

[kessler@xxxxxxxxxxxxx (Billy Kessler)]

Use climatological_axes to make a monthly climatology
(average annual cycle). Then regrid back to daily.
After all, unless you have REALLY a lot of years, a
daily climatology does not make much sense. For example, 
the average of 43 November 22nds (1957-99) does not give 
a good estimate of the climatological November 22nd.
The chances are it will differ significantly from the
average Nov 21st over those 43 years, due to not having
a sufficiently long record to average over the randomness 
of weather. A better estimate is formed by linearly 
interpolating the monthly averages to daily.

Do this by:

yes? let myvar_cl = myvar[gt=month_reg@mod]
yes? let myvar_cl_daily = myvar_cl[gt=myvar]

Now myvar_cl_daily has exactly the same structure as
the original myvar, which is useful in that you can 
subtract to get daily anomalies from the average
annual cycle:

yes? let myvar_anom = myvar - myvar_cl_daily

Billy K