[ferret_users] Seasonal Averages

["nadia.rehman" <nadia.rehman@xxxxxxxxxxxx>]

Hi Users, 

I have got an idea through some e-mails which calculate the seasonal average in a very short way i.e. 

use filename.nc 
use climatological_axes 
cancel data climatological_axes 
let ta_clm=ta[d=1,gt=month_reg@mod] 
fill ta_clm[l=6:9@ave] 

And i also calculated this by a crude way for summer average as: 

use "SRF-ECHAM-1960-70.nc" 
let mod1=ta[l=6:9@ave] 
let mod2=ta[l=18:21@ave] 
let mod3=ta[l=30:33@ave] 
let mod4=ta[l=42:45@ave] 
let mod5=ta[l=54:57@ave] 
let mod6=ta[l=66:69@ave] 
let mod7=ta[l=78:81@ave] 
let mod8=ta[l=90:93@ave] 
let mod9=ta[l=102:105@ave] 
let mod10=ta[l=114:117@ave] 
let mod11=mod1+mod2+mod3+mod4+mod5+mod6+mod7+mod8+mod9+mod10 
let modave=mod11/10 
fill/x=60:95/y=5:40/nolabel modave-273.16 

But the result of both are different.I am using monthly data set. 

Thanks, 

Cheers 


Nadia Rehman Yousafzai
Research Fellow (Climate Section)
Global Change Impact Studies Centre
First Floor Saudi Pak tower 
Blue Area Islamabad Pakistan
Ph :   92-51-9219785
Fax : 92-51-9219787  




________________________________________________________________
Sent via the WebMail system at gcisc.org.pk