PRO plot_waccm_wam_v_prof_phase, zplevWACCM, zplevWAM, uphsWACCM, uphsWAM

; 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/vphs_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
;
tmppWACCM1 = REFORM(uphsWACCM[71,0,36,0:80])
tmppWACCM2 = REFORM(uphsWACCM[72,1,36,0:80])
tmppWACCM3 = REFORM(uphsWACCM[73,2,23,0:80])
tmppWACCM4 = REFORM(uphsWACCM[74,3,28,0:80])

tmppWACCMIn1 = REVERSE(tmppWACCM1)
tmppWACCMIn2 = REVERSE(tmppWACCM2)
tmppWACCMIn3 = REVERSE(tmppWACCM3)
tmppWACCMIn4 = REVERSE(tmppWACCM4)

zplevWACCM1 = REVERSE(zplevWACCM)

tmppWACCM1 = interpol(tmppWACCMIn1, zplevWACCM1, zplevWAM)
tmppWACCM2 = interpol(tmppWACCMIn2, zplevWACCM1, zplevWAM)
tmppWACCM3 = interpol(tmppWACCMIn3, zplevWACCM1, zplevWAM)
tmppWACCM4 = interpol(tmppWACCMIn4, zplevWACCM1, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWAM1 = REFORM(uphsWAM[95,0,36,0:149])
tmppWAM2 = REFORM(uphsWAM[96,1,36,0:149])
tmppWAM3 = REFORM(uphsWAM[97,2,23,0:149])
tmppWAM4 = REFORM(uphsWAM[98,3,28,0:149])

zplevWAM1 = zplevWAM[0:149]

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

oplot, tmppWAM2, zplevWAM1, linestyle=2

plot, tmppWACCM2, zplevWAM1,  /ylog, xrange=[0,24], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Diurnal Tide Phase March 20S',$
         xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmppWAM2, zplevWAM1, linestyle=2

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

oplot, tmppWAM3, zplevWAM1, linestyle=2

plot, tmppWACCM3, zplevWAM1,  /ylog, xrange=[5,12], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Semidiurnal Tide Phase March 45S',$
         xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM3, zplevWAM1, linestyle=2

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

oplot, tmppWAM4, zplevWAM1, linestyle=2

plot, tmppWACCM4, zplevWAM1,  /ylog, xrange=[0,8], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Terdiurnal Tide Phase March 35S',$
         xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM4, zplevWAM1, linestyle=2

;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWACCM1 = REFORM(uphsWACCM[71,0,46,0:80])
tmppWACCM2 = REFORM(uphsWACCM[72,1,46,0:80])
tmppWACCM3 = REFORM(uphsWACCM[73,2,46,0:80])
tmppWACCM4 = REFORM(uphsWACCM[74,3,46,0:80])

tmppWACCMIn1 = REVERSE(tmppWACCM1)
tmppWACCMIn2 = REVERSE(tmppWACCM2)
tmppWACCMIn3 = REVERSE(tmppWACCM3)
tmppWACCMIn4 = REVERSE(tmppWACCM4)

zplevWACCM1 = REVERSE(zplevWACCM)

tmppWACCM1 = interpol(tmppWACCMIn1, zplevWACCM1, zplevWAM)
tmppWACCM2 = interpol(tmppWACCMIn2, zplevWACCM1, zplevWAM)
tmppWACCM3 = interpol(tmppWACCMIn3, zplevWACCM1, zplevWAM)
tmppWACCM4 = interpol(tmppWACCMIn4, zplevWACCM1, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWAM1 = REFORM(uphsWAM[95,0,46,0:149])
tmppWAM2 = REFORM(uphsWAM[96,1,46,0:149])
tmppWAM3 = REFORM(uphsWAM[97,2,46,0:149])
tmppWAM4 = REFORM(uphsWAM[98,3,46,0:149])

zplevWAM1 = zplevWAM[0:149]

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

oplot, tmppWAM2, zplevWAM1, linestyle=2

plot, tmppWACCM2, zplevWAM1,  /ylog, xrange=[0,24], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Diurnal Tide Phase March Equator ',$
	 xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3

oplot, tmppWAM2, zplevWAM1, linestyle=2

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

oplot, tmppWAM3, zplevWAM1, linestyle=2

plot, tmppWACCM3, zplevWAM1,  /ylog, xrange=[0,12], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Semidiurnal Tide Phase March Equator ',$
	 xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmppWAM3, zplevWAM1, linestyle=2

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

oplot, tmppWAM4, zplevWAM1, linestyle=2

plot, tmppWACCM4, zplevWAM1,  /ylog, xrange=[0,8], xstyle = 1,$
         yrange = [1e3,1e-2],ystyle = 1,title = 'WACCM/WAM Meridional Wind Migrating Terdiurnal Tide Phase March Equator ',$
	 xtitle = 'Meridional Wind Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmppWAM4, zplevWAM1, linestyle=2

DEVICE, /CLOSE

SET_PLOT, 'x'

END