PRO plot_waccm_wam_tuvphs_wn5frq5_latlev, latWACCM, zplevWACCM, tphsWACCM, uphsWACCM, vphsWACCM, latWAM, zplevWAM, tphsWAM, uphsWAM, vphsWAM

SET_PLOT, 'PS'

;psfile = '/local/i/joemci/SW/WAM/tuvphs_March_Mean_WN5Frq5_latlev_WACCM.ps'
psfile = 'tuvphs_March_Mean_WN5Frq5_latlev_WACCM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39

xTitle = 'Latitude (Degrees)'
yTitle = 'Pressure (hPa)'
!P.CHARSIZE = 1.2

pTop = 1.E-9
;pBot = 1.E-1
pBot = 1000.0
;
; WACCM temperature westward wavenumber 5 and zonal frequency 5
;
tmppWACCM5 = REFORM(tphsWACCM[76,5,*,0:80])

nlat = N_ELEMENTS(latWACCM)
for ilat = 0,nlat-1 do begin

  temp = REFORM(tmppWACCM5[ilat,*])
  temp = REVERSE(temp)
  tmppWACCM5[ilat,*] = temp
 
endfor

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Temperature Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
; WACCM zonal wind westward wavenumber 5 and zonal frequency 5
;
tmppWACCM5 = REFORM(uphsWACCM[76,5,*,0:80])

tmppWACCM5 = REVERSE(tmppWACCM5)

nlat = N_ELEMENTS(latWACCM)
for ilat = 0,nlat-1 do begin

  temp = REFORM(tmppWACCM5[ilat,*])
  temp = REVERSE(temp)
  tmppWACCM5[ilat,*] = temp
 
endfor

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Zonal Wind Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
; WACCM meridional wind westward wavenumber 5 and zonal frequency 5
;
tmppWACCM5 = REFORM(vphsWACCM[76,5,*,0:80])

tmppWACCM5 = REVERSE(tmppWACCM5)

nlat = N_ELEMENTS(latWACCM)
for ilat = 0,nlat-1 do begin

  temp = REFORM(tmppWACCM5[ilat,*])
  temp = REVERSE(temp)
  tmppWACCM5[ilat,*] = temp
 
endfor

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Meridional Wind Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWACCM5[*,*],latWACCM, zplevWACCM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

DEVICE, /CLOSE

SET_PLOT, 'x'

SET_PLOT, 'PS'

;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8

;psfile = '/home/joemci/MackStuff/MyVersion/WAM/tuvphs_March_Mean_WN5Frq5_latlev_WAM.ps'
psfile = 'tuvphs_March_Mean_WN5Frq5_latlev_WAM.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39
;
; WAM temperature westward wavenumber 5 and zonal frequency 5 
;
tmppWAM5 = REFORM(tphsWAM[100,5,*,0:149])

zplevWAM1 = zplevWAM
 
Title = ' WAM Temperature Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
; WAM zonal wind westward wavenumber 5 and zonal frequency 5
;
tmppWAM5 = REFORM(uphsWAM[100,5,*,0:149])

zplevWAM1 = zplevWAM
 
Title = ' WAM Zonal Wind Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
; WAM meridional wind westward wavenumber 5 and zonal frequency 5
;
tmppWAM5 = REFORM(vphsWAM[100,5,*,0:149])

zplevWAM1 = zplevWAM
 
Title = ' WAM Meridional Wind Phase Migrating 4.8 Hour Zonal Wavenumber 5 March Mean '
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10, /fill,c_colors=BYTSCL(indgen(10)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, tmppWAM5[*,*],latWAM, zplevWAM1, nlev=10,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

DEVICE, /CLOSE

SET_PLOT, 'x'

END