PRO plot_wam_tuv_prof_jun, zplevWAM, tspcWAM, uspcWAM, vspcWAM

; 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_June_Mean_Prof_WAM_WNs2-7.ps'
psfile = 'tuvamp_June_Mean_Prof_WAM_WNs2-7.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WAM migrating wavenumbers 2-7
;
tmpaWAM2 = REFORM(tspcWAM[97,2,55,0:149])
tmpaWAM3 = REFORM(tspcWAM[98,3,55,0:149])
tmpaWAM4 = REFORM(tspcWAM[99,4,55,0:149])
tmpaWAM5 = REFORM(tspcWAM[100,5,55,0:149])
tmpaWAM6 = REFORM(tspcWAM[101,6,55,0:149])
tmpaWAM7 = REFORM(tspcWAM[102,7,55,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,32], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Migrating Tides WNs 2-7 June Mean 16N',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM6, zplevWAM1, color=8
oplot, tmpaWAM7, zplevWAM1, color=8, linestyle=2
xyouts, /NORMAL, 0.65, 0.30, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WAM migrating wavenumbers 2-7
;
tmpaWAM2 = REFORM(uspcWAM[97,2,55,0:149])
tmpaWAM3 = REFORM(uspcWAM[98,3,55,0:149])
tmpaWAM4 = REFORM(uspcWAM[99,4,55,0:149])
tmpaWAM5 = REFORM(uspcWAM[100,5,55,0:149])
tmpaWAM6 = REFORM(tspcWAM[101,6,55,0:149])
tmpaWAM7 = REFORM(tspcWAM[102,7,55,0:149])

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,26], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating Tides WNs 2-7 June Mean 16N',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM6, zplevWAM1, color=8
oplot, tmpaWAM7, zplevWAM1, color=8, linestyle=2
xyouts, /NORMAL, 0.65, 0.30, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WAM migrating wavenumbers 2-7
;
tmpaWAM2 = REFORM(vspcWAM[97,2,55,0:149])
tmpaWAM3 = REFORM(vspcWAM[98,3,55,0:149])
tmpaWAM4 = REFORM(vspcWAM[99,4,55,0:149])
tmpaWAM6 = REFORM(tspcWAM[101,6,55,0:149])
tmpaWAM7 = REFORM(tspcWAM[102,7,55,0:149])
tmpaWAM5 = REFORM(vspcWAM[100,5,55,0:149])

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,26], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Migrating Tides WNs 2-7 June Mean 16N',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM6, zplevWAM1, color=8
oplot, tmpaWAM7, zplevWAM1, color=8, linestyle=2
xyouts, /NORMAL, 0.65, 0.30, 'Wavenumber 2 (solid) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.27, 'Wavenumber 3 (dashed) ', charsize=1.5
xyouts, /NORMAL, 0.65, 0.24, 'Wavenumber 4 (solid) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.21, 'Wavenumber 5 (dashed) ', color=3, charsize=1.5
xyouts, /NORMAL, 0.65, 0.18, 'Wavenumber 6 (solid) ', color=8, charsize=1.5
xyouts, /NORMAL, 0.65, 0.15, 'Wavenumber 7 (dashed) ', color=8, charsize=1.5
;
; WAM non-migrating wavenumbers 1-5
;
;tmpaWAM1 = REFORM(tspcWAM[94,1,55,0:149])
;tmpaWAM2 = REFORM(tspcWAM[93,2,55,0:149])
;tmpaWAM3 = REFORM(tspcWAM[92,3,55,0:149])
;tmpaWAM4 = REFORM(tspcWAM[91,4,55,0:149])
;tmpaWAM5 = REFORM(tspcWAM[90,5,55,0:149])
;
;zplevWAM1 = zplevWAM[0:149]
;
;plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,10], xstyle = 1,$
;	  yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Non-Migrating WNs 2-5 June Mean 16N',$
;	  xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	  charsize = 1.3 
;
;oplot, tmpaWAM3, zplevWAM1, color=3
;oplot, tmpaWAM4, zplevWAM1, color=4
;oplot, tmpaWAM5, zplevWAM1, 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
;;
;; WAM non-migrating wavenumbers 1-5
;;
;tmpaWAM1 = REFORM(uspcWAM[94,1,55,0:149])
;tmpaWAM2 = REFORM(uspcWAM[93,2,55,0:149])
;tmpaWAM3 = REFORM(uspcWAM[92,3,55,0:149])
;tmpaWAM4 = REFORM(uspcWAM[91,4,55,0:149])
;tmpaWAM5 = REFORM(uspcWAM[90,5,55,0:149])
;
;plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,5], xstyle = 1,$
;	  yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Non-Migrating WNs 2-5 June Mean 16N',$
;	  xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	  charsize = 1.3 
;
;oplot, tmpaWAM3, zplevWAM1, color=3
;oplot, tmpaWAM4, zplevWAM1, color=4
;oplot, tmpaWAM5, zplevWAM1, 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
;;
;; WAM non-migrating wavenumbers 1-5
;;
;tmpaWAM1 = REFORM(vspcWAM[94,1,55,0:149])
;tmpaWAM2 = REFORM(vspcWAM[93,2,55,0:149])
;tmpaWAM3 = REFORM(vspcWAM[92,3,55,0:149])
;tmpaWAM4 = REFORM(vspcWAM[91,4,55,0:149])
;tmpaWAM5 = REFORM(vspcWAM[90,5,55,0:149])
;
;plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,10], xstyle = 1,$
;	  yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Non-Migrating WNs 2-5 June Mean 16N',$
;	  xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	  charsize = 1.3 
;
;oplot, tmpaWAM3, zplevWAM1, color=3
;oplot, tmpaWAM4, zplevWAM1, color=4
;oplot, tmpaWAM5, zplevWAM1, 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