PRO plot_waccm_wam_tuv_wvn, tnWACCM, tnWAM, unWACCM, unWAM, vnWACCM, vnWAM

; 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)

SET_PLOT, 'PS'
;
;Temperature amplitudes
;
;psfile = '/local/i/joemci/SW/WAM/twvn_March_Mean_Eq_WACCM_WAM.ps'
psfile = 'twvn_March_Mean_Eq_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmp1 = REFORM(tnWACCM[*,47,4,0])
tmp2 = REFORM(tnWACCM[*,47,4,1])
tmp3 = REFORM(tnWACCM[*,47,4,2])
tmp4 = REFORM(tnWACCM[*,47,4,3])
tmp5 = REFORM(tnWACCM[*,47,4,4])
tmp6 = REFORM(tnWACCM[*,47,4,5])
tmp7 = REFORM(tnWACCM[*,47,4,6])
tmp8 = REFORM(tnWACCM[*,47,4,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,140,0])
tmp2 = REFORM(tnWAM[*,46,140,3])
tmp3 = REFORM(tnWAM[*,46,140,6])
tmp4 = REFORM(tnWAM[*,46,140,9])
tmp5 = REFORM(tnWAM[*,46,140,12])
tmp6 = REFORM(tnWAM[*,46,140,15])
tmp7 = REFORM(tnWAM[*,46,140,18])
tmp8 = REFORM(tnWAM[*,46,140,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 = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5


plot, indgen(13),tmpa2[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa3[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa4[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa5[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa6[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa7[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa8[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,30.0],ystyle = 1,title = 'Temperature March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

DEVICE, /CLOSE

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

;psfile = '/home/joemci/MackStuff/MyVersion/WAM/uwvn_March_Mean_Eq_WACCM_WAM.ps'
psfile = 'uwvn_March_Mean_Eq_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmpu1 = REFORM(unWACCM[*,47,4,0])
tmpu2 = REFORM(unWACCM[*,47,4,1])
tmpu3 = REFORM(unWACCM[*,47,4,2])
tmpu4 = REFORM(unWACCM[*,47,4,3])
tmpu5 = REFORM(unWACCM[*,47,4,4])
tmpu6 = REFORM(unWACCM[*,47,4,5])
tmpu7 = REFORM(unWACCM[*,47,4,6])
tmpu8 = REFORM(unWACCM[*,47,4,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,140,0])
tmpu2 = REFORM(unWAM[*,46,140,3])
tmpu3 = REFORM(unWAM[*,46,140,6])
tmpu4 = REFORM(unWAM[*,46,140,9])
tmpu5 = REFORM(unWAM[*,46,140,12])
tmpu6 = REFORM(unWAM[*,46,140,15])
tmpu7 = REFORM(unWAM[*,46,140,18])
tmpu8 = REFORM(unWAM[*,46,140,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 = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa2[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa3[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa4[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa5[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa6[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa7[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

plot, indgen(13),tmpa8[0:12],xrange=[1,12],xstyle=1,$
         yrange = [0.0,40.0],title = 'Zonal Wind March 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 = '/home/joemci/MackStuff/MyVersion/WAM/vwvn_March_Mean_Eq_WACCM_WAM.ps'
psfile = 'vwvn_March_Mean_Eq_WACCM_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmpv1 = REFORM(vnWACCM[*,47,4,0])
tmpv2 = REFORM(vnWACCM[*,47,4,1])
tmpv3 = REFORM(vnWACCM[*,47,4,2])
tmpv4 = REFORM(vnWACCM[*,47,4,3])
tmpv5 = REFORM(vnWACCM[*,47,4,4])
tmpv6 = REFORM(vnWACCM[*,47,4,5])
tmpv7 = REFORM(vnWACCM[*,47,4,6])
tmpv8 = REFORM(vnWACCM[*,47,4,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[*,46,140,0])
tmpv2 = REFORM(vnWAM[*,46,140,3])
tmpv3 = REFORM(vnWAM[*,46,140,6])
tmpv4 = REFORM(vnWAM[*,46,140,9])
tmpv5 = REFORM(vnWAM[*,46,140,12])
tmpv6 = REFORM(vnWAM[*,46,140,15])
tmpv7 = REFORM(vnWAM[*,46,140,18])
tmpv8 = REFORM(vnWAM[*,46,140,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 = 'Meridional Wind March Mean WNs 1-12 Amplitude 285km UT00 Equator',$
	 xtitle = 'Wavenumber', ytitle = 'Meridional Wind Amplitude',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

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

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

xyouts, /NORMAL, 0.65, 0.88, 'WACCM (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.85, 'WAM (Dotted) ', charsize=1.5

DEVICE, /CLOSE

SET_PLOT, 'x'

END