PRO plot_waccm_wam_u_prof, zplevWACCM, zplevWAM, uspcWACCM, uspcWAM

; 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/tcomp_March_21.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_21_MSDTT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MSTDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MSTDT1.ps'

psfile = '/home/joemci/MackStuff/MyVersion/WAM/uamp_March_Mean_Prof_WACCM_WAM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/WAM/tamp_March_Mean_Prof_WACCM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/WAM/tamp_March_Mean_ZM_WAM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_Mean_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MDT1.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpaWACCM1 = REFORM(uspcWACCM[71,0,36,0:80])
tmpaWACCM2 = REFORM(uspcWACCM[72,1,36,0:80])
tmpaWACCM3 = REFORM(uspcWACCM[73,2,23,0:80])
tmpaWACCM4 = REFORM(uspcWACCM[74,3,28,0:80])

tmpaWACCMIn1 = REVERSE(tmpaWACCM1)
tmpaWACCMIn2 = REVERSE(tmpaWACCM2)
tmpaWACCMIn3 = REVERSE(tmpaWACCM3)
tmpaWACCMIn4 = REVERSE(tmpaWACCM4)

zplevWACCM1 = REVERSE(zplevWACCM)

tmpaWACCM1 = interpol(tmpaWACCMIn1, zplevWACCM1, zplevWAM)
tmpaWACCM2 = interpol(tmpaWACCMIn2, zplevWACCM1, zplevWAM)
tmpaWACCM3 = interpol(tmpaWACCMIn3, zplevWACCM1, zplevWAM)
tmpaWACCM4 = interpol(tmpaWACCMIn4, zplevWACCM1, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpaWAM1 = REFORM(uspcWAM[95,0,36,0:149])
tmpaWAM2 = REFORM(uspcWAM[96,1,36,0:149])
tmpaWAM3 = REFORM(uspcWAM[97,2,23,0:149])
tmpaWAM4 = REFORM(uspcWAM[98,3,28,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[0,70], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Zonal Mean March Mean 20S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[0,70], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Zonal Mean March Mean 20S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[0,75], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Diurnal Tide March Mean 20S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[0,30], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Diurnal Tide March Mean 20S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[0,45], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Semidiurnal Tide March Mean 45S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[0,10], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Semidiurnal Tide March Mean 45S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

plot, tmpaWACCM4, zplevWAM1,  /ylog, xrange=[0,30], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Terdiurnal Tide March Mean 35S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM4, zplevWAM1, linestyle=2

plot, tmpaWACCM4, zplevWAM1,  /ylog, xrange=[0,3], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Migrating Terdiurnal Tide March Mean 35S',$
         xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWAM4, zplevWAM1, linestyle=2
;
;  Do plots at the equator
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpaWACCM1 = REFORM(uspcWACCM[71,0,46,0:80])
tmpaWACCM2 = REFORM(uspcWACCM[72,1,46,0:80])
tmpaWACCM3 = REFORM(uspcWACCM[73,2,46,0:80])
tmpaWACCM4 = REFORM(uspcWACCM[74,3,46,0:80])

tmpaWACCMIn1 = REVERSE(tmpaWACCM1)
tmpaWACCMIn2 = REVERSE(tmpaWACCM2)
tmpaWACCMIn3 = REVERSE(tmpaWACCM3)
tmpaWACCMIn4 = REVERSE(tmpaWACCM4)

zplevWACCM1 = REVERSE(zplevWACCM)

tmpaWACCM1 = interpol(tmpaWACCMIn1, zplevWACCM1, zplevWAM)
tmpaWACCM2 = interpol(tmpaWACCMIn2, zplevWACCM1, zplevWAM)
tmpaWACCM3 = interpol(tmpaWACCMIn3, zplevWACCM1, zplevWAM)
tmpaWACCM4 = interpol(tmpaWACCMIn4, zplevWACCM1, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpaWAM1 = REFORM(uspcWAM[95,0,46,0:149])
tmpaWAM2 = REFORM(uspcWAM[96,1,46,0:149])
tmpaWAM3 = REFORM(uspcWAM[97,2,46,0:149])
tmpaWAM4 = REFORM(uspcWAM[98,3,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[0,65], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Zonal Mean March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[0,65], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Zonal Mean March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[0,75], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Diurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[0,10], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Diurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Semidiurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[0,8], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Semidiurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

plot, tmpaWACCM4, zplevWAM1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Terdiurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM4, zplevWAM1, linestyle=2

plot, tmpaWACCM4, zplevWAM1,  /ylog, xrange=[0,5], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM and WAM Zonal Wind Migrating Terdiurnal Tide March Mean Equator ',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM4, zplevWAM1, linestyle=2

DEVICE, /CLOSE

SET_PLOT, 'x'

END