dw1_t_cm=total(dw1_t,4)/20.
dw1_t_mc=fltarr(96,81,12)
dw1_t_mc(*,*,0)=total(dw1_t_cm(*,*,0:30),3)/31.
dw1_t_mc(*,*,1)=total(dw1_t_cm(*,*,31:31+27),3)/28.
dw1_t_mc(*,*,2)=total(dw1_t_cm(*,*,59:59+30),3)/31
dw1_t_mc(*,*,3)=total(dw1_t_cm(*,*,90:90+29),3)/30
dw1_t_mc(*,*,4)=total(dw1_t_cm(*,*,120:120+30),3)/31
dw1_t_mc(*,*,5)=total(dw1_t_cm(*,*,151:151+29),3)/30
dw1_t_mc(*,*,6)=total(dw1_t_cm(*,*,181:181+30),3)/31
dw1_t_mc(*,*,7)=total(dw1_t_cm(*,*,212:212+30),3)/31
dw1_t_mc(*,*,8)=total(dw1_t_cm(*,*,243:243+29),3)/30
dw1_t_mc(*,*,9)=total(dw1_t_cm(*,*,273:273+30),3)/31
dw1_t_mc(*,*,10)=total(dw1_t_cm(*,*,304:304+29),3)/30
dw1_t_mc(*,*,11)=total(dw1_t_cm(*,*,334:334+30),3)/31

plot,dw1_t_mc(11,24:30,0),zlp(24:30),xrange=[0,10],yrange=[80,100],$
xtitle='K',ytitle='Approx. height (km)',$
title='DW1 Amplitude at 69S, Jan, WACCMX'
for iyr=0,19 do begin
   for imn=1,30 do begin
      oplot,dw1_t(11,24:30,imn,iyr),zlp(24:30),color=150,thick=0.8
   endfor
endfor
oplot,dw1_t_mc(11,24:30,0),zlp(24:30),thick=3

plot,findgen(12)+1,dw1_t_mc(11,27,*),xst=1,$
xtitle='Month',ytitle='K',title='DW1 climatology at 69S 90km, WACCMX'

plot,lat,dw1_t_mc(*,27,0),xst=1,$
xtitle='Lat (deg)',ytitle='K',title='DW1 climatology at 90km, Jan., WACCMX'

sw2_t_cm=total(sw2_t,4)/20.
sw2_t_mc=fltarr(96,81,12)
sw2_t_mc(*,*,0)=total(sw2_t_cm(*,*,0:30),3)/31.
sw2_t_mc(*,*,1)=total(sw2_t_cm(*,*,31:31+27),3)/28.
sw2_t_mc(*,*,2)=total(sw2_t_cm(*,*,59:59+30),3)/31
sw2_t_mc(*,*,3)=total(sw2_t_cm(*,*,90:90+29),3)/30
sw2_t_mc(*,*,4)=total(sw2_t_cm(*,*,120:120+30),3)/31
sw2_t_mc(*,*,5)=total(sw2_t_cm(*,*,151:151+29),3)/30
sw2_t_mc(*,*,6)=total(sw2_t_cm(*,*,181:181+30),3)/31
sw2_t_mc(*,*,7)=total(sw2_t_cm(*,*,212:212+30),3)/31
sw2_t_mc(*,*,8)=total(sw2_t_cm(*,*,243:243+29),3)/30
sw2_t_mc(*,*,9)=total(sw2_t_cm(*,*,273:273+30),3)/31
sw2_t_mc(*,*,10)=total(sw2_t_cm(*,*,304:304+29),3)/30
sw2_t_mc(*,*,11)=total(sw2_t_cm(*,*,334:334+30),3)/31

plot,sw2_t_mc(11,24:30,0),zlp(24:30),xrange=[0,5],yrange=[80,100],$
xtitle='K',ytitle='Approx. height (km)',$
title='SW2 Amplitude at 69S, Jan, WACCMX'
for iyr=0,19 do begin
   for imn=1,30 do begin
      oplot,sw2_t(11,24:30,imn,iyr),zlp(24:30),color=150,thick=0.8
   endfor
endfor
oplot,sw2_t_mc(11,24:30,0),zlp(24:30),thick=3

plot,findgen(12)+1,sw2_t_mc(11,27,*),xst=1,$
xtitle='Month',ytitle='K',title='SW2 climatology at 69S 90km, WACCMX'

plot,lat,sw2_t_mc(*,27,0),xst=1,$
xtitle='Lat (deg)',ytitle='K',title='SW2 climatology at 90km, Jan., WACCMX'

end