[Thread Prev][Thread Next][Index]
Re: climatology question
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
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement