PRO plot_WACCM_tuv_prof, zplevWACCM, tspcWACCM, uspcWACCM, vspcWACCM

; 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

!P.CHARSIZE = 2.5
!P.CHARTHICK = 2.5
!P.THICK = 5.5

;psfile = '/local/i/joemci/SW/WAM/tuvamp_March_Mean_Prof_WACCM_WNs2-7.ps'
psfile = 'tuvamp_March_Mean_Prof_WACCM_WNs2-7.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WACCM migrating wavenumbers 2-7
;
tmpaWACCM2 = REFORM(tspcWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(tspcWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(tspcWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(tspcWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(tspcWACCM[77,6,47,0:80])
tmpaWACCM7 = REFORM(tspcWACCM[78,7,47,0:80])

tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,12], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating Tides WNs 2-7 March Mean Equator',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM6, zplevWACCM1, color=8
oplot, tmpaWACCM7, zplevWACCM1, color=8, linestyle=2
xyouts, /NORMAL, 0.15, 0.90, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.87, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.84, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.81, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.78, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.15, 0.75, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WACCM migrating wavenumbers 2-7
;
tmpaWACCM2 = REFORM(uspcWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(uspcWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(uspcWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(uspcWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(tspcWACCM[77,6,47,0:80])
tmpaWACCM7 = REFORM(tspcWACCM[78,7,47,0:80])

tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)

plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,22], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Migrating Tides WNs 2-7 March Mean Equator',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM6, zplevWACCM1, color=8
oplot, tmpaWACCM7, zplevWACCM1, color=8, linestyle=2
xyouts, /NORMAL, 0.15, 0.90, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.87, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.84, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.81, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.78, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.15, 0.75, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WACCM migrating wavenumbers 2-7
;
tmpaWACCM2 = REFORM(vspcWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(vspcWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(vspcWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(vspcWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(tspcWACCM[77,6,47,0:80])
tmpaWACCM7 = REFORM(tspcWACCM[78,7,47,0:80])

tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)
tmpaWACCM7 = REVERSE(tmpaWACCM7)

plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,8], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Migrating Tides WNs 2-7 March Mean Equator',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM3, zplevWACCM1, linestyle=2
oplot, tmpaWACCM4, zplevWACCM1, color=3
oplot, tmpaWACCM5, zplevWACCM1, color=3, linestyle=2
oplot, tmpaWACCM6, zplevWACCM1, color=8
oplot, tmpaWACCM7, zplevWACCM1, color=8, linestyle=2
xyouts, /NORMAL, 0.15, 0.90, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.87, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.15, 0.84, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.81, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.15, 0.78, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.15, 0.75, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WACCM non-migrating wavenumbers 1-5
;
;tmpaWACCM1 = REFORM(tspcWACCM[70,1,46,0:80])
;tmpaWACCM2 = REFORM(tspcWACCM[69,2,46,0:80])
;tmpaWACCM3 = REFORM(tspcWACCM[68,3,46,0:80])
;tmpaWACCM4 = REFORM(tspcWACCM[67,4,46,0:80])
;tmpaWACCM5 = REFORM(tspcWACCM[66,5,46,0:80])
;
;tmpaWACCM1 = REVERSE(tmpaWACCM1)
;tmpaWACCM2 = REVERSE(tmpaWACCM2)
;tmpaWACCM3 = REVERSE(tmpaWACCM3)
;tmpaWACCM4 = REVERSE(tmpaWACCM4)
;tmpaWACCM5 = REVERSE(tmpaWACCM5)
;
;plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,6], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Non-Migrating WNs 2-5 1993 March Mean Equator',$
;         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;         charsize = 1.3 
;
;oplot, tmpaWACCM3, zplevWACCM1, color=3
;oplot, tmpaWACCM4, zplevWACCM1, color=4
;oplot, tmpaWACCM5, zplevWACCM1, color=8
;xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 2 ', charsize=1.5
;xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 3 ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 4 ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 5 ', color=8, charsize=1.5
;
; WACCM non-migrating wavenumbers 1-5
;
;tmpaWACCM1 = REFORM(uspcWACCM[70,1,46,0:80])
;tmpaWACCM2 = REFORM(uspcWACCM[69,2,46,0:80])
;tmpaWACCM3 = REFORM(uspcWACCM[68,3,46,0:80])
;tmpaWACCM4 = REFORM(uspcWACCM[67,4,46,0:80])
;tmpaWACCM5 = REFORM(uspcWACCM[66,5,46,0:80])
;
;tmpaWACCM1 = REVERSE(tmpaWACCM1)
;tmpaWACCM2 = REVERSE(tmpaWACCM2)
;tmpaWACCM3 = REVERSE(tmpaWACCM3)
;tmpaWACCM4 = REVERSE(tmpaWACCM4)
;tmpaWACCM5 = REVERSE(tmpaWACCM5)
;
;plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,8], xstyle = 1,$
;	  yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Non-Migrating WNs 2-5 1993 March Mean Equator',$
;	  xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	  charsize = 1.3 
;
;oplot, tmpaWACCM3, zplevWACCM1, color=3
;oplot, tmpaWACCM4, zplevWACCM1, color=4
;oplot, tmpaWACCM5, zplevWACCM1, color=8
;xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 2 ', charsize=1.5
;xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 3 ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 4 ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 5 ', color=8, charsize=1.5
;;
;; WACCM non-migrating wavenumbers 1-5
;;
;tmpaWACCM1 = REFORM(vspcWACCM[70,1,46,0:80])
;tmpaWACCM2 = REFORM(vspcWACCM[69,2,46,0:80])
;tmpaWACCM3 = REFORM(vspcWACCM[68,3,46,0:80])
;tmpaWACCM4 = REFORM(vspcWACCM[67,4,46,0:80])
;tmpaWACCM5 = REFORM(vspcWACCM[66,5,46,0:80])
;
;tmpaWACCM1 = REVERSE(tmpaWACCM1)
;tmpaWACCM2 = REVERSE(tmpaWACCM2)
;tmpaWACCM3 = REVERSE(tmpaWACCM3)
;tmpaWACCM4 = REVERSE(tmpaWACCM4)
;tmpaWACCM5 = REVERSE(tmpaWACCM5)
;
;plot, tmpaWACCM2, zplevWACCM1,  /xlog,/ylog, xrange=[1E-04,7], xstyle = 1,$
;	  yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Non-Migrating WNs 2-5 1993 March Mean Equator',$
;	  xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	  charsize = 1.3 
;
;oplot, tmpaWACCM3, zplevWACCM1, color=3
;oplot, tmpaWACCM4, zplevWACCM1, color=4
;oplot, tmpaWACCM5, zplevWACCM1, color=8
;xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 2 ', charsize=1.5
;xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 3 ', color=3, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 4 ', color=4, charsize=1.5
;xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 5 ', color=8, charsize=1.5

DEVICE, /CLOSE

SET_PLOT, 'x'

END