use "
test.nc"
define axis/t/units=months/edges/t="15-jan-2000":"15-dec-2010":1/calendar=360_days tmonthly
let var_anm=var_data[gt=tmonthly@asn] ! variable anomaly is the new variable name
let month=mod(l[gt=var_anm],12)
let mask1 = if month ge 3 and month le 5 then 1
let mask2 = if month ge 6 and month le 9 then 1
!<<<<<<<<<<<<<<<<<<<<<<<<<masking for mam & jjas<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<!
let var_mam = var_anm*mask1
let var_jjas = var_anm*mask2
define axis/t/units=years/edges/t="15-jan-2000":"15-dec-2010":1/calendar=360_days tyearly
let var_mam_yearly = var_mam[gt=tyearly@ave]
let var_jjas_yearly = var_jjas[gt=tyearly@ave]
list var_mam_yearly
list var_jjas_yearly