PRO plot_waccm_wam_tuv_wvn1to12_june, tnWACCM, tnWAM, unWACCM, unWAM, vnWACCM, vnWAM, zgm

; This procedure was created on Tuesday, August 31st, 2010.
; This procedure will produce line plots of temperature or wind amplitudes (y-axis) and wavenumber (x-axis)

dims = SIZE(tnWACCM)

nX = dims(1)
nY = dims(2)
nT = dims(4)

tZSubWACCM = FLTARR(nX,nY,nT)
uZSubWACCM = FLTARR(nX,nY,nT)
vZSubWACCM = FLTARR(nX,nY,nT)
;
;  Find altitude index for this time
;    
FOR iLn=0,nX-1 DO BEGIN

  FOR iLt=0,nY-1 DO BEGIN

;    FOR iT = 0, nT-1 DO BEGIN

      zWLnLtT = zgm[iLn,iLt,*,0]
      zInd = WHERE(zWLnLtT GE 275.0 AND zWLnLtT LE 295.0)

      IF zInd[0] EQ -1 THEN BEGIN
  
        MESSAGE, 'Cannot find altitude index in WACCM'
     
      ENDIF ELSE BEGIN
  
        iZ = zInd[0]
     
      ENDELSE

      tZSubWACCM(iLn,iLt,*) = tnWACCM(iLn,iLt,iZ,*)
      uZSubWACCM(iLn,iLt,*) = unWACCM(iLn,iLt,iZ,*)
      vZSubWACCM(iLn,iLt,*) = vnWACCM(iLn,iLt,iZ,*)

;    ENDFOR

  ENDFOR

ENDFOR

SET_PLOT, 'PS'
;
;Temperature amplitudes
;
;psfile = '/local/i/joemci/SW/WAM/twvn1-12_June_Mean_Eq_WACCM_WAM.ps'
psfile = 'twvn1-12_June_Mean_Eq_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmp1 = REFORM(tZSubWACCM[*,47,0])
tmp2 = REFORM(tZSubWACCM[*,47,1])
tmp3 = REFORM(tZSubWACCM[*,47,2])
tmp4 = REFORM(tZSubWACCM[*,47,3])
tmp5 = REFORM(tZSubWACCM[*,47,4])
tmp6 = REFORM(tZSubWACCM[*,47,5])
tmp7 = REFORM(tZSubWACCM[*,47,6])
tmp8 = REFORM(tZSubWACCM[*,47,7])

tmp1 = REVERSE(tmp1)
tmp2 = REVERSE(tmp2)
tmp3 = REVERSE(tmp3)
tmp4 = REVERSE(tmp4)
tmp5 = REVERSE(tmp5)
tmp6 = REVERSE(tmp6)
tmp7 = REVERSE(tmp7)
tmp8 = REVERSE(tmp8)

tmpa1 = ABS(FFT(tmp1,-1))
tmpa2 = ABS(FFT(tmp2,-1))
tmpa3 = ABS(FFT(tmp3,-1))
tmpa4 = ABS(FFT(tmp4,-1))
tmpa5 = ABS(FFT(tmp5,-1))
tmpa6 = ABS(FFT(tmp6,-1))
tmpa7 = ABS(FFT(tmp7,-1))
tmpa8 = ABS(FFT(tmp8,-1))
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmp1 = REFORM(tnWAM[*,46,136,0])
tmp2 = REFORM(tnWAM[*,46,136,3])
tmp3 = REFORM(tnWAM[*,46,136,6])
tmp4 = REFORM(tnWAM[*,46,136,9])
tmp5 = REFORM(tnWAM[*,46,136,12])
tmp6 = REFORM(tnWAM[*,46,136,15])
tmp7 = REFORM(tnWAM[*,46,136,18])
tmp8 = REFORM(tnWAM[*,46,136,21])

tmpa1WAM = ABS(FFT(tmp1,-1))
tmpa2WAM = ABS(FFT(tmp2,-1))
tmpa3WAM = ABS(FFT(tmp3,-1))
tmpa4WAM = ABS(FFT(tmp4,-1))
tmpa5WAM = ABS(FFT(tmp5,-1))
tmpa6WAM = ABS(FFT(tmp6,-1))
tmpa7WAM = ABS(FFT(tmp7,-1))
tmpa8WAM = ABS(FFT(tmp8,-1))

plot, indgen(13),tmpa1[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT00 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa1WAM[0:12],linestyle=1

plot, indgen(13),tmpa2[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT03 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa2WAM[0:12],linestyle=1

plot, indgen(13),tmpa3[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT06 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa3WAM[0:12],linestyle=1

plot, indgen(13),tmpa4[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT09 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa4WAM[0:12],linestyle=1

plot, indgen(13),tmpa5[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT12 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa5WAM[0:12],linestyle=1

plot, indgen(13),tmpa6[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT15 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa6WAM[0:12],linestyle=1

plot, indgen(13),tmpa7[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT18 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa7WAM[0:12],linestyle=1

plot, indgen(13),tmpa8[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'WACCM/WAM Temperature June Mean WNs 1-12 Amplitude 285km UT21 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Temperature Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa8WAM[0:12],linestyle=1

DEVICE, /CLOSE

SET_PLOT, 'x'
;
;Zonal Wind amplitudes
;
SET_PLOT, 'PS'

;psfile = '/local/i/joemci/SW/WAM/uwvn1-12_June_Mean_Eq_WACCM_WAM.ps'
psfile = 'uwvn1-12_June_Mean_Eq_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmpu1 = REFORM(uZSubWACCM[*,47,0])
tmpu2 = REFORM(uZSubWACCM[*,47,1])
tmpu3 = REFORM(uZSubWACCM[*,47,2])
tmpu4 = REFORM(uZSubWACCM[*,47,3])
tmpu5 = REFORM(uZSubWACCM[*,47,4])
tmpu6 = REFORM(uZSubWACCM[*,47,5])
tmpu7 = REFORM(uZSubWACCM[*,47,6])
tmpu8 = REFORM(uZSubWACCM[*,47,7])

tmpu1 = REVERSE(tmpu1)
tmpu2 = REVERSE(tmpu2)
tmpu3 = REVERSE(tmpu3)
tmpu4 = REVERSE(tmpu4)
tmpu5 = REVERSE(tmpu5)
tmpu6 = REVERSE(tmpu6)
tmpu7 = REVERSE(tmpu7)
tmpu8 = REVERSE(tmpu8)

tmpa1 = ABS(FFT(tmpu1,-1))
tmpa2 = ABS(FFT(tmpu2,-1))
tmpa3 = ABS(FFT(tmpu3,-1))
tmpa4 = ABS(FFT(tmpu4,-1))
tmpa5 = ABS(FFT(tmpu5,-1))
tmpa6 = ABS(FFT(tmpu6,-1))
tmpa7 = ABS(FFT(tmpu7,-1))
tmpa8 = ABS(FFT(tmpu8,-1))
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpu1 = REFORM(unWAM[*,46,136,0])
tmpu2 = REFORM(unWAM[*,46,136,3])
tmpu3 = REFORM(unWAM[*,46,136,6])
tmpu4 = REFORM(unWAM[*,46,136,9])
tmpu5 = REFORM(unWAM[*,46,136,12])
tmpu6 = REFORM(unWAM[*,46,136,15])
tmpu7 = REFORM(unWAM[*,46,136,18])
tmpu8 = REFORM(unWAM[*,46,136,21])

tmpa1WAM = ABS(FFT(tmpu1,-1))
tmpa2WAM = ABS(FFT(tmpu2,-1))
tmpa3WAM = ABS(FFT(tmpu3,-1))
tmpa4WAM = ABS(FFT(tmpu4,-1))
tmpa5WAM = ABS(FFT(tmpu5,-1))
tmpa6WAM = ABS(FFT(tmpu6,-1))
tmpa7WAM = ABS(FFT(tmpu7,-1))
tmpa8WAM = ABS(FFT(tmpu8,-1))

plot, indgen(13),tmpa1[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT00 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa1WAM[0:12],linestyle=1

plot, indgen(13),tmpa2[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT03 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa2WAM[0:12],linestyle=1

plot, indgen(13),tmpa3[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT06 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa3WAM[0:12],linestyle=1

plot, indgen(13),tmpa4[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT09 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa4WAM[0:12],linestyle=1

plot, indgen(13),tmpa5[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT12 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa5WAM[0:12],linestyle=1

plot, indgen(13),tmpa6[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT15 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa6WAM[0:12],linestyle=1

plot, indgen(13),tmpa7[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT18 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa7WAM[0:12],linestyle=1

plot, indgen(13),tmpa8[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'WACCM/WAM Zonal Wind June Mean WNs 1-12 Amplitude 285km UT21 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Zonal Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa8WAM[0:12],linestyle=1

DEVICE, /CLOSE

SET_PLOT, 'x'
;
;Meridional Wind amplitudes
;
SET_PLOT, 'PS'

;psfile = '/local/i/joemci/SW/WAM/vwvn1-12_June_Mean_12N_WACCM_WAM.ps'
psfile = 'vwvn1-12_June_Mean_12N_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmpv1 = REFORM(vZSubWACCM[*,41,0])
tmpv2 = REFORM(vZSubWACCM[*,41,1])
tmpv3 = REFORM(vZSubWACCM[*,41,2])
tmpv4 = REFORM(vZSubWACCM[*,41,3])
tmpv5 = REFORM(vZSubWACCM[*,41,4])
tmpv6 = REFORM(vZSubWACCM[*,41,5])
tmpv7 = REFORM(vZSubWACCM[*,41,6])
tmpv8 = REFORM(vZSubWACCM[*,41,7])

tmpv1 = REVERSE(tmpv1)
tmpv2 = REVERSE(tmpv2)
tmpv3 = REVERSE(tmpv3)
tmpv4 = REVERSE(tmpv4)
tmpv5 = REVERSE(tmpv5)
tmpv6 = REVERSE(tmpv6)
tmpv7 = REVERSE(tmpv7)
tmpv8 = REVERSE(tmpv8)

tmpa1 = ABS(FFT(tmpv1,-1))
tmpa2 = ABS(FFT(tmpv2,-1))
tmpa3 = ABS(FFT(tmpv3,-1))
tmpa4 = ABS(FFT(tmpv4,-1))
tmpa5 = ABS(FFT(tmpv5,-1))
tmpa6 = ABS(FFT(tmpv6,-1))
tmpa7 = ABS(FFT(tmpv7,-1))
tmpa8 = ABS(FFT(tmpv8,-1))
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpv1 = REFORM(vnWAM[*,40,136,0])
tmpv2 = REFORM(vnWAM[*,40,136,3])
tmpv3 = REFORM(vnWAM[*,40,136,6])
tmpv4 = REFORM(vnWAM[*,40,136,9])
tmpv5 = REFORM(vnWAM[*,40,136,12])
tmpv6 = REFORM(vnWAM[*,40,136,15])
tmpv7 = REFORM(vnWAM[*,40,136,18])
tmpv8 = REFORM(vnWAM[*,40,136,21])

tmpa1WAM = ABS(FFT(tmpv1,-1))
tmpa2WAM = ABS(FFT(tmpv2,-1))
tmpa3WAM = ABS(FFT(tmpv3,-1))
tmpa4WAM = ABS(FFT(tmpv4,-1))
tmpa5WAM = ABS(FFT(tmpv5,-1))
tmpa6WAM = ABS(FFT(tmpv6,-1))
tmpa7WAM = ABS(FFT(tmpv7,-1))
tmpa8WAM = ABS(FFT(tmpv8,-1))

plot, indgen(13),tmpa1[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT00 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa1WAM[0:12],linestyle=1

plot, indgen(13),tmpa2[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT03 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa2WAM[0:12],linestyle=1

plot, indgen(13),tmpa3[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT06 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa3WAM[0:12],linestyle=1

plot, indgen(13),tmpa4[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT09 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa4WAM[0:12],linestyle=1

plot, indgen(13),tmpa5[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT12 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa5WAM[0:12],linestyle=1

plot, indgen(13),tmpa6[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT15 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa6WAM[0:12],linestyle=1

plot, indgen(13),tmpa7[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT18 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa7WAM[0:12],linestyle=1

plot, indgen(13),tmpa8[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,10.0],ystyle = 1,title = 'WACCM/WAM Meridional Wind June Mean WNs 1-12 Amplitude 285km UT21 12N',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

oplot, indgen(13),tmpa8WAM[0:12],linestyle=1

DEVICE, /CLOSE

SET_PLOT, 'x'

END