'reinit'
'enable print grads.outfile'
say '** OUTPUT FILE:grads.outfile'
'clear'
'set grads off'

 
year1=1948
while(year1<=2010)
'sdfopen air.sig995.'year1'.nc'
'set lat -106 -106'
'set lon 45 45'
lsleapyear=math_fmod(year1,4)
if lsleapyear=0
 LastDayFeb=29Feb
else
 LastDayFeb=28Feb
endif
*year1=year1+1
*endwhile


*assign l# to the data range t1 thru t2
'define a1=ave(air.1,time=01jan%year1,time=31jan%year%1)'
'define a2=ave(air.1,time=01feb%year1,time=LastDayFeb%year1)'
'define a3=ave(air.1,time=01mar%year1,time=31mar%year1)'
'define a4=ave(air.1,time=01apr%year1,time=30apr%year1)'
'define a5=ave(air.1,time=01may%year1,time=31may%year1)'
'define a6=ave(air.1,time=01jun%year1,time=30jun%year1)'
'define a7=ave(air.1,time=01jul%year1,time=31jul%year1)'
'define a8=ave(air.1,time=01aug%year1,time=31aug%year1)'
'define a9=ave(air.1,time=01sep%year1,time=30sep%year1)'
'define a10=ave(air.1,time=01oct%year1,time=31oct%year1)'
'define a11=ave(air.1,time=01nov%year1,time=30nov%year1)'
'define a12=ave(air.1,time=01dec%year1,time=31dec%year1)'

*Calculate standard deviation for l1-l12
'define sda1=sqrt(ave(pow(ave(air.1,time=01jan%year1,time=31jan%year%1)-air.1,2),time=01jan%year1,time=31jan%year%1))'
'define sda2=sqrt(ave(pow(ave(air.1,time=01feb%year1,time=LastDayFeb%year1)-air.1,2),time=01feb%year1,time=LastDayFeb%year1))'
'define sda3=sqrt(ave(pow(ave(air.1,time=01mar%year1,time=31mar%year1)-air.1,2),time=01mar%year1,t=31mar%year1))'
'define sda4=sqrt(ave(pow(ave(air.1,time=01apr%year1,time=30apr%year1)-air.1,2),time=01apr%year1,time=30apr%year1))'
'define sda5=sqrt(ave(pow(ave(air.1,time=01may%year1,time=31may%year1)-air.1,2),time=01may%year1,time=31may%year1))'
'define sda6=sqrt(ave(pow(ave(air.1,time=01jun%year1,time=30jun%year1)-air.1,2),time=01jun%year1,time=30jun%year1))'
'define sda7=sqrt(ave(pow(ave(air.1,time=01jul%year1,time=31jul%year1)-air.1,2),time=01jul%year1,time=31jul%year1))'
'define sda8=sqrt(ave(pow(ave(air.1,time=01aug%year1,time=31aug%year1)-air.1,2),time=01aug%year1,time=31aug%year1))'
'define sda9=sqrt(ave(pow(ave(air.1,time=01sep%year1,time=30sep%year1)-air.1,2),time=01sep%year1,time=30sep%year1))'
'define sda10=sqrt(ave(pow(ave(air.1,time=01oct%year1,time=31oct%year1)-air.1,2),time=01oct%year1,time=31oct%year1))'
'define sda11=sqrt(ave(pow(ave(air.1,time=01nov%year1,time=30nov%year1)-air.1,2),time=01nov%year1,time=30nov%year1))'
'define sda12=sqrt(ave(pow(ave(air.1,time=01dec%year1,time=31dec%year1)-air.1,2),time=01dec%year1,time=31dec%year1))'
*
year1=year1+1
endwhile