PRO plot_WAM_tuv_wn5frq1to12_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

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

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WAM temperature westward wavenumber 5 and zonal frequencies 1-12 
;
tmpaWAM1 = REFORM(tspcWAM[100,1,55,0:149])
tmpaWAM2 = REFORM(tspcWAM[100,2,55,0:149])
tmpaWAM3 = REFORM(tspcWAM[100,3,55,0:149])
tmpaWAM4 = REFORM(tspcWAM[100,4,55,0:149])
tmpaWAM5 = REFORM(tspcWAM[100,5,55,0:149])
tmpaWAM6 = REFORM(tspcWAM[100,6,55,0:149])
tmpaWAM7 = REFORM(tspcWAM[100,7,55,0:149])
tmpaWAM8 = REFORM(tspcWAM[100,8,55,0:149])
tmpaWAM9 = REFORM(tspcWAM[100,9,55,0:149])
tmpaWAM10 = REFORM(tspcWAM[100,10,55,0:149])
tmpaWAM11 = REFORM(tspcWAM[100,11,55,0:149])
tmpaWAM12 = REFORM(tspcWAM[100,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,12], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Westward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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
;
; WAM temperature eastward wavenumber 5 and zonal frequencies 1-12
;
tmpaWAM1 = REFORM(tspcWAM[90,1,55,0:149])
tmpaWAM2 = REFORM(tspcWAM[90,2,55,0:149])
tmpaWAM3 = REFORM(tspcWAM[90,3,55,0:149])
tmpaWAM4 = REFORM(tspcWAM[90,4,55,0:149])
tmpaWAM5 = REFORM(tspcWAM[90,5,55,0:149])
tmpaWAM6 = REFORM(tspcWAM[90,6,55,0:149])
tmpaWAM7 = REFORM(tspcWAM[90,7,55,0:149])
tmpaWAM8 = REFORM(tspcWAM[90,8,55,0:149])
tmpaWAM9 = REFORM(tspcWAM[90,9,55,0:149])
tmpaWAM10 = REFORM(tspcWAM[90,10,55,0:149])
tmpaWAM11 = REFORM(tspcWAM[90,11,55,0:149])
tmpaWAM12 = REFORM(tspcWAM[90,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,12], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Eastward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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
;
; WAM zonal wind westward wavenumber 5 and zonal frequencies 1-12
;
tmpaWAM1 = REFORM(uspcWAM[100,1,55,0:149])
tmpaWAM2 = REFORM(uspcWAM[100,2,55,0:149])
tmpaWAM3 = REFORM(uspcWAM[100,3,55,0:149])
tmpaWAM4 = REFORM(uspcWAM[100,4,55,0:149])
tmpaWAM5 = REFORM(uspcWAM[100,5,55,0:149])
tmpaWAM6 = REFORM(uspcWAM[100,6,55,0:149])
tmpaWAM7 = REFORM(uspcWAM[100,7,55,0:149])
tmpaWAM8 = REFORM(uspcWAM[100,8,55,0:149])
tmpaWAM9 = REFORM(uspcWAM[100,9,55,0:149])
tmpaWAM10 = REFORM(uspcWAM[100,10,55,0:149])
tmpaWAM11 = REFORM(uspcWAM[100,11,55,0:149])
tmpaWAM12 = REFORM(uspcWAM[100,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,24], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Westward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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
;
; WAM zonal wind eastward wavenumber 5 and zonal frequencies 1-12
;
tmpaWAM1 = REFORM(uspcWAM[90,1,55,0:149])
tmpaWAM2 = REFORM(uspcWAM[90,2,55,0:149])
tmpaWAM3 = REFORM(uspcWAM[90,3,55,0:149])
tmpaWAM4 = REFORM(uspcWAM[90,4,55,0:149])
tmpaWAM5 = REFORM(uspcWAM[90,5,55,0:149])
tmpaWAM6 = REFORM(uspcWAM[90,6,55,0:149])
tmpaWAM7 = REFORM(uspcWAM[90,7,55,0:149])
tmpaWAM8 = REFORM(uspcWAM[90,8,55,0:149])
tmpaWAM9 = REFORM(uspcWAM[90,9,55,0:149])
tmpaWAM10 = REFORM(uspcWAM[90,10,55,0:149])
tmpaWAM11 = REFORM(uspcWAM[90,11,55,0:149])
tmpaWAM12 = REFORM(uspcWAM[90,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,24], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Eastward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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
;
; WAM meridional wind westward wavenumber 5 and zonal frequencies 1-12
;
tmpaWAM1 = REFORM(vspcWAM[100,1,55,0:149])
tmpaWAM2 = REFORM(vspcWAM[100,2,55,0:149])
tmpaWAM3 = REFORM(vspcWAM[100,3,55,0:149])
tmpaWAM4 = REFORM(vspcWAM[100,4,55,0:149])
tmpaWAM5 = REFORM(vspcWAM[100,5,55,0:149])
tmpaWAM6 = REFORM(vspcWAM[100,6,55,0:149])
tmpaWAM7 = REFORM(vspcWAM[100,7,55,0:149])
tmpaWAM8 = REFORM(vspcWAM[100,8,55,0:149])
tmpaWAM9 = REFORM(vspcWAM[100,9,55,0:149])
tmpaWAM10 = REFORM(vspcWAM[100,10,55,0:149])
tmpaWAM11 = REFORM(vspcWAM[100,11,55,0:149])
tmpaWAM12 = REFORM(vspcWAM[100,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,12], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Westward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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
;
; WAM temperature eastward wavenumber 5 and zonal frequencies 1-12
;
tmpaWAM1 = REFORM(vspcWAM[90,1,55,0:149])
tmpaWAM2 = REFORM(vspcWAM[90,2,55,0:149])
tmpaWAM3 = REFORM(vspcWAM[90,3,55,0:149])
tmpaWAM4 = REFORM(vspcWAM[90,4,55,0:149])
tmpaWAM5 = REFORM(vspcWAM[90,5,55,0:149])
tmpaWAM6 = REFORM(vspcWAM[90,6,55,0:149])
tmpaWAM7 = REFORM(vspcWAM[90,7,55,0:149])
tmpaWAM8 = REFORM(vspcWAM[90,8,55,0:149])
tmpaWAM9 = REFORM(vspcWAM[90,9,55,0:149])
tmpaWAM10 = REFORM(vspcWAM[90,10,55,0:149])
tmpaWAM11 = REFORM(vspcWAM[90,11,55,0:149])
tmpaWAM12 = REFORM(vspcWAM[90,12,55,0:149])

zplevWAM1 = zplevWAM

plot, tmpaWAM1, zplevWAM1,  /ylog, xrange=[0,12], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Eastward WN 5 Freqs 1-12 1993 June Mean 15N',$
         xtitle = 'Meridional Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM2, zplevWAM1, linestyle=1
oplot, tmpaWAM3, zplevWAM1, linestyle=2
oplot, tmpaWAM4, zplevWAM1, color=3
oplot, tmpaWAM5, zplevWAM1, color=3, linestyle=1
oplot, tmpaWAM6, zplevWAM1, color=3, linestyle=2
oplot, tmpaWAM7, zplevWAM1, color=4
oplot, tmpaWAM8, zplevWAM1, color=4, linestyle=1
oplot, tmpaWAM9, zplevWAM1, color=4, linestyle=2
oplot, tmpaWAM10, zplevWAM1, color=8
oplot, tmpaWAM11, zplevWAM1, color=8, linestyle=1
oplot, tmpaWAM12, zplevWAM1, 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