PRO plot_waccm_wam_tuvamp_wn4frq4_latlev_jun, latWACCM, zplevWACCM, tspcWACCM, uspcWACCM, vspcWACCM, latWAM, zplevWAM, tspcWAM, uspcWAM, vspcWAM

SET_PLOT, 'PS'

;psfile = '/local/i/joemci/SW/WAM/tuvamp_June_Mean_WN4Frq4_latlev_WACCM.ps'
psfile = 'tuvamp_June_Mean_WN4Frq4_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 4 and zonal frequency 4
;
tmpaWACCM4 = REFORM(tspcWACCM[75,4,*,0:80])

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

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

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Temperature Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWACCM4[*,*],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, tmpaWACCM4[*,*],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
;
tmpaWACCM4 = REFORM(uspcWACCM[75,5,*,0:80])

tmpaWACCM4 = REVERSE(tmpaWACCM4)

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

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

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Zonal Wind Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWACCM4[*,*],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, tmpaWACCM4[*,*],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 4 and zonal frequency 4
;
tmpaWACCM4 = REFORM(vspcWACCM[75,5,*,0:80])

tmpaWACCM4 = REVERSE(tmpaWACCM4)

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

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

zplevWACCM1 = REVERSE(zplevWACCM)
 
Title = ' WACCM-X Meridional Wind Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWACCM4[*,*],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, tmpaWACCM4[*,*],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/tuvamp_June_Mean_WN4Frq4_latlev_WAM.ps'
psfile = 'tuvamp_June_Mean_WN4Frq4_latlev_WAM.ps'

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

zplevWAM1 = zplevWAM
 
Title = ' WAM Temperature Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWAM4[*,*],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, tmpaWAM4[*,*],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 4 and zonal frequency 4
;
tmpaWAM4 = REFORM(uspcWAM[99,5,*,0:149])

zplevWAM1 = zplevWAM
 
Title = ' WAM Zonal Wind Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWAM4[*,*],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, tmpaWAM4[*,*],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 4 and zonal frequency 4
;
tmpaWAM4 = REFORM(vspcWAM[99,5,*,0:149])

zplevWAM1 = zplevWAM
 
Title = ' WAM Meridional Wind Amplitude Migrating 6 Hour Zonal Wavenumber 4 June Mean '
CONTOUR, tmpaWAM4[*,*],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, tmpaWAM4[*,*],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