PRO plot_waccm_wam_t_prof_zmdev, zplevWACCM, zplevWAM, tnWACCM, tnWAM, unWACCM, unWAM, vnWACCM, vnWAM

; 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/tn_un_vn_March_Mean_Prof_WACCM_WAM_eqzmdev.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(tnWACCM[0,47,0:80,0])
tmpaWACCM2 = REFORM(unWACCM[0,47,0:80,0])
tmpaWACCM3 = REFORM(vnWACCM[0,47,0:80,0])
;tmpaWACCM4 = REFORM(vnWACCM[74,3,47,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(tnWAM[0,46,0:149,0])
tmpaWAM2 = REFORM(unWAM[0,46,0:149,0])
tmpaWAM3 = REFORM(vnWAM[0,46,0:149,0])
;tmpaWAM4 = REFORM(vnWAM[98,3,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[-75,35], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Temperature Deviation March Mean 00 LST Equator ',$
	 xtitle = 'Temperature (K)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM1, zplevWAM1,  /ylog, xrange=[-15,20], xstyle = 1,$
         yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Temperature Deviation March Mean 00 LST Equator ',$
	 xtitle = 'Temperature (K)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

oplot, tmpaWAM1, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[-50,90], xstyle = 1,$
	 yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Zonal Wind Deviation March Mean 00 LST Equator ',$
	xtitle = 'Zonal Wind (m/s)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM2, zplevWAM1,  /ylog, xrange=[-20,20], xstyle = 1,$
	 yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Zonal Wind Deviation March Mean 00 LST Equator ',$
	xtitle = 'Zonal Wind (m/s)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	charsize = 1.3

oplot, tmpaWAM2, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[-30,35], xstyle = 1,$
	 yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM & WAM Meridional Wind Deviation March Mean 00 LST Equator ',$
	xtitle = 'Meridional Wind (m/s)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

plot, tmpaWACCM3, zplevWAM1,  /ylog, xrange=[-10,15], xstyle = 1,$
	 yrange = [1e3,1e-4],ystyle = 1,title = 'WACCM & WAM Meridional Wind Deviation March Mean 00 LST Equator ',$
	xtitle = 'Meridional Wind (m/s)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	charsize = 1.3 

oplot, tmpaWAM3, zplevWAM1, linestyle=2

DEVICE, /CLOSE

SET_PLOT, 'x'

END