PRO plot_WACCM_tuv_wn4frq1to12_prof_jun, zplevWACCM, tspcWACCM, uspcWACCM, vspcWACCM


;newcolor

SET_PLOT, 'PS'

;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8

;psfile = '/local/i/joemci/SW/WAM/tuvamp_June_Mean_WN4Frq1to12_Prof_WACCM.ps'
psfile = 'tuvamp_June_Mean_WN4Frq1to12_Prof_WACCM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WACCM temperature westward wavenumber 4 and zonal frequencies 1-12 
;
tmpaWACCM1 = REFORM(tspcWACCM[75,1,56,0:80])
tmpaWACCM2 = REFORM(tspcWACCM[75,2,56,0:80])
tmpaWACCM3 = REFORM(tspcWACCM[75,3,56,0:80])
tmpaWACCM4 = REFORM(tspcWACCM[75,4,56,0:80])
tmpaWACCM5 = REFORM(tspcWACCM[75,5,56,0:80])
tmpaWACCM6 = REFORM(tspcWACCM[75,6,56,0:80])
tmpaWACCM7 = REFORM(tspcWACCM[75,7,56,0:80])
tmpaWACCM8 = REFORM(tspcWACCM[75,8,56,0:80])
tmpaWACCM9 = REFORM(tspcWACCM[75,9,56,0:80])
tmpaWACCM10 = REFORM(tspcWACCM[75,10,56,0:80])
tmpaWACCM11 = REFORM(tspcWACCM[75,11,56,0:80])
tmpaWACCM12 = REFORM(tspcWACCM[75,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,14], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Westward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM temperature eastward wavenumber 4 and zonal frequencies 1-12
;
tmpaWACCM1 = REFORM(tspcWACCM[67,1,56,0:80])
tmpaWACCM2 = REFORM(tspcWACCM[67,2,56,0:80])
tmpaWACCM3 = REFORM(tspcWACCM[67,3,56,0:80])
tmpaWACCM4 = REFORM(tspcWACCM[67,4,56,0:80])
tmpaWACCM5 = REFORM(tspcWACCM[67,5,56,0:80])
tmpaWACCM6 = REFORM(tspcWACCM[67,6,56,0:80])
tmpaWACCM7 = REFORM(tspcWACCM[67,7,56,0:80])
tmpaWACCM8 = REFORM(tspcWACCM[67,8,56,0:80])
tmpaWACCM9 = REFORM(tspcWACCM[67,9,56,0:80])
tmpaWACCM10 = REFORM(tspcWACCM[67,10,56,0:80])
tmpaWACCM11 = REFORM(tspcWACCM[67,11,56,0:80])
tmpaWACCM12 = REFORM(tspcWACCM[67,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,14], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Eastward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM zonal wind westward wavenumber 4 and zonal frequencies 1-12
;
tmpaWACCM1 = REFORM(uspcWACCM[75,1,56,0:80])
tmpaWACCM2 = REFORM(uspcWACCM[75,2,56,0:80])
tmpaWACCM3 = REFORM(uspcWACCM[75,3,56,0:80])
tmpaWACCM4 = REFORM(uspcWACCM[75,4,56,0:80])
tmpaWACCM5 = REFORM(uspcWACCM[75,5,56,0:80])
tmpaWACCM6 = REFORM(uspcWACCM[75,6,56,0:80])
tmpaWACCM7 = REFORM(uspcWACCM[75,7,56,0:80])
tmpaWACCM8 = REFORM(uspcWACCM[75,8,56,0:80])
tmpaWACCM9 = REFORM(uspcWACCM[75,9,56,0:80])
tmpaWACCM10 = REFORM(uspcWACCM[75,10,56,0:80])
tmpaWACCM11 = REFORM(uspcWACCM[75,11,56,0:80])
tmpaWACCM12 = REFORM(uspcWACCM[75,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,20], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Westward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM zonal wind eastward wavenumber 4 and zonal frequencies 1-12
;
tmpaWACCM1 = REFORM(uspcWACCM[67,1,56,0:80])
tmpaWACCM2 = REFORM(uspcWACCM[67,2,56,0:80])
tmpaWACCM3 = REFORM(uspcWACCM[67,3,56,0:80])
tmpaWACCM4 = REFORM(uspcWACCM[67,4,56,0:80])
tmpaWACCM5 = REFORM(uspcWACCM[67,5,56,0:80])
tmpaWACCM6 = REFORM(uspcWACCM[67,6,56,0:80])
tmpaWACCM7 = REFORM(uspcWACCM[67,7,56,0:80])
tmpaWACCM8 = REFORM(uspcWACCM[67,8,56,0:80])
tmpaWACCM9 = REFORM(uspcWACCM[67,9,56,0:80])
tmpaWACCM10 = REFORM(uspcWACCM[67,10,56,0:80])
tmpaWACCM11 = REFORM(uspcWACCM[67,11,56,0:80])
tmpaWACCM12 = REFORM(uspcWACCM[67,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,20], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Eastward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM meridional wind westward wavenumber 4 and zonal frequencies 1-12
;
tmpaWACCM1 = REFORM(vspcWACCM[75,1,56,0:80])
tmpaWACCM2 = REFORM(vspcWACCM[75,2,56,0:80])
tmpaWACCM3 = REFORM(vspcWACCM[75,3,56,0:80])
tmpaWACCM4 = REFORM(vspcWACCM[75,4,56,0:80])
tmpaWACCM5 = REFORM(vspcWACCM[75,5,56,0:80])
tmpaWACCM6 = REFORM(vspcWACCM[75,6,56,0:80])
tmpaWACCM7 = REFORM(vspcWACCM[75,7,56,0:80])
tmpaWACCM8 = REFORM(vspcWACCM[75,8,56,0:80])
tmpaWACCM9 = REFORM(vspcWACCM[75,9,56,0:80])
tmpaWACCM10 = REFORM(vspcWACCM[75,10,56,0:80])
tmpaWACCM11 = REFORM(vspcWACCM[75,11,56,0:80])
tmpaWACCM12 = REFORM(vspcWACCM[75,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,14], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Westward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM temperature eastward wavenumber 4 and zonal frequencies 1-12
;
tmpaWACCM1 = REFORM(vspcWACCM[67,1,56,0:80])
tmpaWACCM2 = REFORM(vspcWACCM[67,2,56,0:80])
tmpaWACCM3 = REFORM(vspcWACCM[67,3,56,0:80])
tmpaWACCM4 = REFORM(vspcWACCM[67,4,56,0:80])
tmpaWACCM5 = REFORM(vspcWACCM[67,5,56,0:80])
tmpaWACCM6 = REFORM(vspcWACCM[67,6,56,0:80])
tmpaWACCM7 = REFORM(vspcWACCM[67,7,56,0:80])
tmpaWACCM8 = REFORM(vspcWACCM[67,8,56,0:80])
tmpaWACCM9 = REFORM(vspcWACCM[67,9,56,0:80])
tmpaWACCM10 = REFORM(vspcWACCM[67,10,56,0:80])
tmpaWACCM11 = REFORM(vspcWACCM[67,11,56,0:80])
tmpaWACCM12 = REFORM(vspcWACCM[67,12,56,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)
tmpaWACCM8 = REVERSE(tmpaWACCM8)
tmpaWACCM9 = REVERSE(tmpaWACCM9)
tmpaWACCM10 = REVERSE(tmpaWACCM10)
tmpaWACCM11 = REVERSE(tmpaWACCM11)
tmpaWACCM12 = REVERSE(tmpaWACCM12)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0,14], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Eastward WN 4 Frqs 1-12 1993 June Mean 15N',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1
oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=1
oplot, tmpaWACCM6, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM7, zplevWACCM1, color=4
oplot, tmpaWACCM8, zplevWACCM1, color=4, linestyle=1
oplot, tmpaWACCM9, zplevWACCM1, color=4, linestyle=2
oplot, tmpaWACCM10, zplevWACCM1, color=8
oplot, tmpaWACCM11, zplevWACCM1, color=8, linestyle=1
oplot, tmpaWACCM12, zplevWACCM1, color=8, linestyle=2

xyouts, /NORMAL, 0.65, 0.48, 'Frequency 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.45, 'Frequency 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.42, 'Frequency 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.39, 'Frequency 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.36, 'Frequency 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.33, 'Frequency 6 (Dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.30, 'Frequency 7 (Solid) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Frequency 8 (Dotted) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Frequency 9 (Dashed) ', color=4, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Frequency 10 (Solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Frequency 11 (Dotted) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Frequency 12 (Dashed) ', color=8, charsize=1.5

DEVICE, /CLOSE

SET_PLOT, 'x'

END