hourly climatology (6-hourly)

For making climatology in 6 hourly interval, new axis should be defined. (ref)

Let’s start!




set memory/size=200
use slp.79-97.nc

def axis/t=0:`365*24`/calendar=noleap/edges/npoints=`365*4`/t0=1-JAN-0001/units=hours/modulo new_axis
let clim=var[gt=new_axis@mod]

save/clo/file=slp_79-97_clim.nc clim



at t=0:`365*24` : unit is hours, so 1 year axis should have 365 (day) * 24 (hour) points.

at calendar=noleap: noleap means no leap year. most model use noleap calendar type.

(If the data considers leap year, then use “julian”.ref)

at t0=1-Jan-0001: define first time as 1 year, January, 1 day.

at @mod: creating climatology (ref)


It’s very powerful function and,

So, Simple 🙂

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.