PRO plot_waccm_t_prof, un, zplev, uspc

; 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/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
;
tmpa1 = REFORM(uspc[71,0,46,0:80])
tmpa2 = REFORM(uspc[72,1,46,0:80])
tmpa3 = REFORM(uspc[73,2,46,0:80])
tmpa4 = REFORM(uspc[74,3,46,0:80])

tmpa1 = REVERSE(tmpa1)
tmpa2 = REVERSE(tmpa2)
tmpa3 = REVERSE(tmpa3)
tmpa4 = REVERSE(tmpa4)

zplev1 = REVERSE(zplev)

plot, tmpa1, zplev1,  /ylog, xrange=[0,1000], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Zonal Mean March Mean ',$
	 xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

plot, tmpa2, zplev1,  /ylog, xrange=[0,70], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating Diurnal Tide March Mean ',$
	 xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3

plot, tmpa3, zplev1,  /ylog, xrange=[0,40], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating Semidiurnal Tide March Mean ',$
	 xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

plot, tmpa4, zplev1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating Terdiurnal Tide March Mean ',$
	 xtitle = 'Temperature Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

DEVICE, /CLOSE

SET_PLOT, 'x'

END