PRO plot_waccm_tuv_wnphs1to6wm_prof, zplevWACCM, tphsWACCM, uphsWACCM, vphsWACCM

; This procedure was created on Tuesday, June 9th, 2009.
; This procedure was created by SOARS Protege McArthur (Mack) Jones Jr.
; This procedure will produce a contour w/ lat. on the x-axis and altitude on the y-axis

;newcolor

SET_PLOT, 'PS'

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

psfile = '/local/i/joemci/SW/WAM/tuvphs_March_Mean_WNPhs1to6WM_Prof_WACCM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WACCM temperature migrating westward wavenumbers 1-12 
;
tmpaWACCM1 = REFORM(tphsWACCM[72,1,47,0:80])
tmpaWACCM2 = REFORM(tphsWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(tphsWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(tphsWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(tphsWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(tphsWACCM[77,6,47,0:80])
;tmpaWACCM7 = REFORM(tphsWACCM[78,7,47,0:80])
;tmpaWACCM8 = REFORM(tphsWACCM[79,8,47,0:80])
;tmpaWACCM9 = REFORM(tphsWACCM[80,9,47,0:80])
;tmpaWACCM10 = REFORM(tphsWACCM[81,10,47,0:80])
;tmpaWACCM11 = REFORM(tphsWACCM[82,11,47,0:80])
;tmpaWACCM12 = REFORM(tphsWACCM[83,12,47,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.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating Westward WNs 1-6 1993 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', 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.35, 0.88, 'WN 1 Freq 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.35, 0.85, 'WN 2 Freq 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.35, 0.82, 'WN 3 Freq 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.35, 0.79, 'WN 4 Freq 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.35, 0.76, 'WN 5 Freq 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.35, 0.73, 'WN 6 Freq 6 (Dashed) ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.30, 'WN 7 Freq 7 (Solid) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.27, 'WN 8 Freq 8 (Dotted) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.24, 'WN 9 Freq 9 (Dashed) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.21, 'WN 10 Freq 10 (Solid) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.18, 'WN 11 Freq 11 (Dotted) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.15, 'WN 12 Freq 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM zonal wind migrating westward wavenumbers 1-12
;
tmpaWACCM1 = REFORM(uphsWACCM[72,1,47,0:80])
tmpaWACCM2 = REFORM(uphsWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(uphsWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(uphsWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(uphsWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(uphsWACCM[77,6,47,0:80])
;tmpaWACCM7 = REFORM(uphsWACCM[78,7,47,0:80])
;tmpaWACCM8 = REFORM(uphsWACCM[79,8,47,0:80])
;tmpaWACCM9 = REFORM(uphsWACCM[80,9,47,0:80])
;tmpaWACCM10 = REFORM(uphsWACCM[81,10,47,0:80])
;tmpaWACCM11 = REFORM(uphsWACCM[82,11,47,0:80])
;tmpaWACCM12 = REFORM(uphsWACCM[83,12,47,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.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Migrating Westward WNs 1-6 1993 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', 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.45, 0.91, 'WN 1 Freq 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.45, 0.88, 'WN 2 Freq 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.45, 0.85, 'WN 3 Freq 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.45, 0.82, 'WN 4 Freq 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.45, 0.79, 'WN 5 Freq 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.45, 0.76, 'WN 6 Freq 6 (Dashed) ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.30, 'WN 7 Freq 7 (Solid) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.27, 'WN 8 Freq 8 (Dotted) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.24, 'WN 9 Freq 9 (Dashed) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.21, 'WN 10 Freq 10 (Solid) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.18, 'WN 11 Freq 11 (Dotted) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.15, 'WN 12 Freq 12 (Dashed) ', color=8, charsize=1.5
;
; WACCM meridional wind migrating westward wavenumbers 1-12
;
tmpaWACCM1 = REFORM(vphsWACCM[72,1,47,0:80])
tmpaWACCM2 = REFORM(vphsWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(vphsWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(vphsWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(vphsWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(vphsWACCM[77,6,47,0:80])
;tmpaWACCM7 = REFORM(vphsWACCM[78,7,47,0:80])
;tmpaWACCM8 = REFORM(vphsWACCM[79,8,47,0:80])
;tmpaWACCM9 = REFORM(vphsWACCM[80,9,47,0:80])
;tmpaWACCM10 = REFORM(vphsWACCM[81,10,47,0:80])
;tmpaWACCM11 = REFORM(vphsWACCM[82,11,47,0:80])
;tmpaWACCM12 = REFORM(vphsWACCM[83,12,47,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.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Migrating Westward WNs 1-6 1993 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', 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.4, 0.84, 'WN 1 Freq 1 (Solid) ', charsize=1.5
xyouts, /NORMAL, 0.4, 0.81, 'WN 2 Freq 2 (Dotted) ', charsize=1.5
xyouts, /NORMAL, 0.4, 0.78, 'WN 3 Freq 3 (Dashed) ', charsize=1.5
xyouts, /NORMAL, 0.4, 0.75, 'WN 4 Freq 4 (Solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.4, 0.72, 'WN 5 Freq 5 (Dotted) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.4, 0.69, 'WN 6 Freq 6 (Dashed) ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.30, 'WN 7 Freq 7 (Solid) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.27, 'WN 8 Freq 8 (Dotted) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.24, 'WN 9 Freq 9 (Dashed) ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.21, 'WN 10 Freq 10 (Solid) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.18, 'WN 11 Freq 11 (Dotted) ', color=8, charsize=1.5
;xyouts, /NORMAL, 0.15, 0.15, 'WN 12 Freq 12 (Dashed) ', color=8, charsize=1.5

DEVICE, /CLOSE

SET_PLOT, 'X'

END