PRO plot_wam_u_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/uamp_March_Mean_Prof_WAM.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[95,0,31,0:149])
tmpa2 = REFORM(uspc[96,1,31,0:149])
tmpa3 = REFORM(uspc[97,2,20,0:149])
tmpa4 = REFORM(uspc[98,3,31,0:149])

zplev1 = zplev[0:149]

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

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

plot, tmpa3, zplev1,  /ylog, xrange=[0,45], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating Semidiurnal Tide March Mean 50S',$
	 xtitle = 'Zonal Wind 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 = 'WAM Zonal Wind Migrating Terdiurnal Tide March Mean 30S',$
	 xtitle = 'Zonal Wind Amplitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3 

tmpa1 = REFORM(uspc[95,0,46,0:149])
tmpa2 = REFORM(uspc[96,1,46,0:149])
tmpa3 = REFORM(uspc[97,2,46,0:149])
tmpa4 = REFORM(uspc[98,3,46,0:149])

zplev1 = zplev[0:149]

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

plot, tmpa2, zplev1,  /ylog, xrange=[0,75], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = '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

plot, tmpa3, zplev1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = '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 

plot, tmpa4, zplev1,  /ylog, xrange=[0,25], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = '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 

DEVICE, /CLOSE

SET_PLOT, 'x'

END