daily climatology

Calculating climatology !! (ref, ref2)

climate data has various time steps, such as hourly, daily, weekly, monthly, years ..

This script is useful for handling the NetCDF file with various time steps.


<Example for daily>

define axis/t=1-jan-1961:31-dec-1990:1/unit=days tlong
let myvar = RANDU(T[gt=tlong])       ! a random dummy data variable (if variable is SST then, SST[gt=tlong])
def axis/t=0:365.2425/edges/npoints=365/t0=1-JAN-0001/units=days/modulo daily1
let myvar_clim=myvar[gt=daily1@mod]
let myvar_anom=myvar-myvar_clim[gt=tlong]


<For drawing and checking data>

set region/t=1-jan-1980:31-dec-1980
set view upper
plot myvar, myvar_clim
set view lower
plot myvar_anom


Not hard huh?

So, Simple 🙂

(if you know ferret little bit)

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.