PRO plot_wam_tuv_prof_sep, 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

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

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WAM migrating wavenumbers 1-5
;
tmpaWAM1 = REFORM(tspcWAM[96,1,46,0:149])
tmpaWAM2 = REFORM(tspcWAM[97,2,46,0:149])
tmpaWAM3 = REFORM(tspcWAM[98,3,46,0:149])
tmpaWAM4 = REFORM(tspcWAM[99,4,46,0:149])
tmpaWAM5 = REFORM(tspcWAM[100,5,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,40], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Migrating WNs 2-5 September Mean Equator',$
         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 migrating wavenumbers 1-5
;
tmpaWAM1 = REFORM(uspcWAM[96,1,46,0:149])
tmpaWAM2 = REFORM(uspcWAM[97,2,46,0:149])
tmpaWAM3 = REFORM(uspcWAM[98,3,46,0:149])
tmpaWAM4 = REFORM(uspcWAM[99,4,46,0:149])
tmpaWAM5 = REFORM(uspcWAM[100,5,46,0:149])

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating WNs 2-5 September Mean Equator',$
         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 migrating wavenumbers 1-5
;
tmpaWAM1 = REFORM(vspcWAM[96,1,46,0:149])
tmpaWAM2 = REFORM(vspcWAM[97,2,46,0:149])
tmpaWAM3 = REFORM(vspcWAM[98,3,46,0:149])
tmpaWAM4 = REFORM(vspcWAM[99,4,46,0:149])
tmpaWAM5 = REFORM(vspcWAM[100,5,46,0:149])

plot, tmpaWAM2, zplevWAM1,  /ylog, xrange=[0,20], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Migrating WNs 2-5 September Mean Equator',$
         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
;
; WAM non-migrating wavenumbers 1-5
;
;tmpaWAM1 = REFORM(tspcWAM[94,1,46,0:149])
;tmpaWAM2 = REFORM(tspcWAM[93,2,46,0:149])
;tmpaWAM3 = REFORM(tspcWAM[92,3,46,0:149])
;tmpaWAM4 = REFORM(tspcWAM[91,4,46,0:149])
;tmpaWAM5 = REFORM(tspcWAM[90,5,46,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 March Mean Equator',$
;	  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,46,0:149])
;tmpaWAM2 = REFORM(uspcWAM[93,2,46,0:149])
;tmpaWAM3 = REFORM(uspcWAM[92,3,46,0:149])
;tmpaWAM4 = REFORM(uspcWAM[91,4,46,0:149])
;tmpaWAM5 = REFORM(uspcWAM[90,5,46,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 March Mean Equator',$
;	  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,46,0:149])
;tmpaWAM2 = REFORM(vspcWAM[93,2,46,0:149])
;tmpaWAM3 = REFORM(vspcWAM[92,3,46,0:149])
;tmpaWAM4 = REFORM(vspcWAM[91,4,46,0:149])
;tmpaWAM5 = REFORM(vspcWAM[90,5,46,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 March Mean Equator',$
;	  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