PRO plot_waccm_wam_tuv_wnphs1to6wm_prof, zplevWACCM, tphsWACCM, uphsWACCM, vphsWACCM, zplevWAM, tphsWAM, uphsWAM, vphsWAM

; This procedure will produce profile plots of phase for migrating wavenumbers 1-6

;newcolor

  !P.CHARSIZE = 1.2
  !P.CHARTHICK = 2.5
  !P.THICK = 5.5

SET_PLOT, 'PS'

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

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

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WAM temperature migrating wavenumbers 1-6
;
tmppWACCM1 = REFORM(tphsWACCM[72,1,47,0:80])
tmppWACCM2 = REFORM(tphsWACCM[73,2,47,0:80])
tmppWACCM3 = REFORM(tphsWACCM[74,3,47,0:80])
tmppWACCM4 = REFORM(tphsWACCM[75,4,47,0:80])
tmppWACCM5 = REFORM(tphsWACCM[76,5,47,0:80])
tmppWACCM6 = REFORM(tphsWACCM[77,6,47,0:80])

tmppWACCM1 = REVERSE(tmppWACCM1)
tmppWACCM2 = REVERSE(tmppWACCM2)
tmppWACCM3 = REVERSE(tmppWACCM3)
tmppWACCM4 = REVERSE(tmppWACCM4)
tmppWACCM5 = REVERSE(tmppWACCM5)
tmppWACCM6 = REVERSE(tmppWACCM6)

zplevWACCM1 = REVERSE(zplevWACCM)
;
; WACCM top down so reverse needed
;
plot, tmppWACCM1, zplevWACCM1, /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM2, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.7, 0.88, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.87,0.89],[0.92,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.7, 0.85, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.87,0.86],[0.92,0.86]], /NORMAL, linestyle = 1 

plot, tmppWACCM3, zplevWACCM1, /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM4, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.7, 0.88, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.87,0.89],[0.92,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.7, 0.85, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.87,0.86],[0.92,0.86]], /NORMAL, linestyle = 1 

plot, tmppWACCM5, zplevWACCM1, /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Temperature Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM6, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.7, 0.88, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.87,0.89],[0.92,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.7, 0.85, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.87,0.86],[0.92,0.86]], /NORMAL, linestyle = 1 
;
; WACCM zonal wind migrating wavenumbers 1-6
;
tmppWACCM1 = REFORM(uphsWACCM[72,1,47,0:80])
tmppWACCM2 = REFORM(uphsWACCM[73,2,47,0:80])
tmppWACCM3 = REFORM(uphsWACCM[74,3,47,0:80])
tmppWACCM4 = REFORM(uphsWACCM[75,4,47,0:80])
tmppWACCM5 = REFORM(uphsWACCM[76,5,47,0:80])
tmppWACCM6 = REFORM(uphsWACCM[77,6,47,0:80])

tmppWACCM1 = REVERSE(tmppWACCM1)
tmppWACCM2 = REVERSE(tmppWACCM2)
tmppWACCM3 = REVERSE(tmppWACCM3)
tmppWACCM4 = REVERSE(tmppWACCM4)
tmppWACCM5 = REVERSE(tmppWACCM5)
tmppWACCM6 = REVERSE(tmppWACCM6)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmppWACCM1, zplevWACCM1, /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM2, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.14, 0.88, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.31,0.89],[0.36,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.14, 0.85, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.31,0.86],[0.36,0.86]], /NORMAL, linestyle = 1 

plot, tmppWACCM3, zplevWACCM1, /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM4, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.14, 0.88, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.31,0.89],[0.36,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.14, 0.85, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.31,0.86],[0.36,0.86]], /NORMAL, linestyle = 1 

plot, tmppWACCM5, zplevWACCM1, /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Zonal Wind Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWACCM6, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.14, 0.88, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.31,0.89],[0.36,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.14, 0.85, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.31,0.86],[0.36,0.86]], /NORMAL, linestyle = 1 
;
; WACCM meridional wind migrating wavenumbers 1-12
;
tmpaWACCM1 = REFORM(vphsWACCM[72,1,47,0:80])
tmpaWACCM2 = REFORM(vphsWACCM[73,2,47,0:80])
tmpaWACCM3 = REFORM(vphsWACCM[74,3,47,0:80])
tmpaWACCM4 = REFORM(vphsWACCM[75,4,47,0:80])
tmpaWACCM5 = REFORM(vphsWACCM[76,5,47,0:80])
tmpaWACCM6 = REFORM(vphsWACCM[77,6,47,0:80])

tmpaWACCM1 = REVERSE(tmpaWACCM1)
tmpaWACCM2 = REVERSE(tmpaWACCM2)
tmpaWACCM3 = REVERSE(tmpaWACCM3)
tmpaWACCM4 = REVERSE(tmpaWACCM4)
tmpaWACCM5 = REVERSE(tmpaWACCM5)
tmpaWACCM6 = REVERSE(tmpaWACCM6)

zplevWACCM1 = REVERSE(zplevWACCM)

plot, tmpaWACCM1, zplevWACCM1,  /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM2, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.2, 0.82, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.37,0.83],[0.42,0.83]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.2, 0.79, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.37,0.80],[0.42,0.80]], /NORMAL, linestyle = 1 

plot, tmpaWACCM3, zplevWACCM1,  /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM4, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.2, 0.88, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.37,0.89],[0.42,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.2, 0.85, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.37,0.86],[0.42,0.86]], /NORMAL, linestyle = 1 

plot, tmpaWACCM5, zplevWACCM1,  /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM Meridional Wind Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmpaWACCM6, zplevWACCM1, linestyle=1

xyouts, /NORMAL, 0.2, 0.88, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.37,0.89],[0.42,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.2, 0.85, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.37,0.86],[0.42,0.86]], /NORMAL, linestyle = 1 

DEVICE, /CLOSE

SET_PLOT, 'X'
;
;  Do WAM plots
;
SET_PLOT, 'PS'

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

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

DEVICE, /COLOR, /LAND, FILENAME = psfile
LOADCT, 39, ncolors = 10
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWAM1 = REFORM(tphsWAM[96,1,46,0:149])
tmppWAM2 = REFORM(tphsWAM[97,2,46,0:149])
tmppWAM3 = REFORM(tphsWAM[98,3,46,0:149])
tmppWAM4 = REFORM(tphsWAM[99,4,46,0:149])
tmppWAM5 = REFORM(tphsWAM[100,5,46,0:149])
tmppWAM6 = REFORM(tphsWAM[101,6,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmppWAM1, zplevWAM1, /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM2, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.25, 0.83, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.42,0.84],[0.47,0.84]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.25, 0.80, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.42,0.81],[0.47,0.81]], /NORMAL, linestyle = 1 

plot, tmppWAM3, zplevWAM1, /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM4, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.25, 0.83, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.42,0.84],[0.47,0.84]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.25, 0.80, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.42,0.81],[0.47,0.81]], /NORMAL, linestyle = 1 

plot, tmppWAM5, zplevWAM1, /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Temperature Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Temperature Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM6, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.25, 0.83, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.42,0.84],[0.47,0.84]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.25, 0.80, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.42,0.81],[0.47,0.81]], /NORMAL, linestyle = 1 
;
; WAM zonal wind migrating wavenumbers 1-6
;
tmppWAM1 = REFORM(uphsWAM[96,1,46,0:149])
tmppWAM2 = REFORM(uphsWAM[97,2,46,0:149])
tmppWAM3 = REFORM(uphsWAM[98,3,46,0:149])
tmppWAM4 = REFORM(uphsWAM[99,4,46,0:149])
tmppWAM5 = REFORM(uphsWAM[100,5,46,0:149])
tmppWAM6 = REFORM(uphsWAM[101,6,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmppWAM1, zplevWAM1, /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM2, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.35, 0.88, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.52,0.89],[0.57,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.35, 0.85, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.52,0.86],[0.57,0.86]], /NORMAL, linestyle = 1 


plot, tmppWAM3, zplevWAM1, /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM4, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.25, 0.88, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.42,0.89],[0.47,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.25, 0.85, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.42,0.86],[0.47,0.86]], /NORMAL, linestyle = 1 


plot, tmppWAM5, zplevWAM1, /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Zonal Wind Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Zonal Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM6, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.45, 0.88, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.62,0.89],[0.67,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.45, 0.85, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.62,0.86],[0.67,0.86]], /NORMAL, linestyle = 1 
;
; WAM meridional wind migrating wavenumbers 1-6
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWAM1 = REFORM(uphsWAM[96,1,46,0:149])
tmppWAM2 = REFORM(uphsWAM[97,2,46,0:149])
tmppWAM3 = REFORM(uphsWAM[98,3,46,0:149])
tmppWAM4 = REFORM(uphsWAM[99,4,46,0:149])
tmppWAM5 = REFORM(uphsWAM[100,5,46,0:149])
tmppWAM6 = REFORM(uphsWAM[101,6,46,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmppWAM1, zplevWAM1,  /ylog, xrange=[0.0,360.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Migrating WNs 1-2 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM2, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.35, 0.88, 'WN 1 Freq 1 ', charsize=1.5
plots, [[0.52,0.89],[0.57,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.35, 0.85, 'WN 2 Freq 2 ', charsize=1.5
plots, [[0.52,0.86],[0.57,0.86]], /NORMAL, linestyle = 1 

plot, tmppWAM3, zplevWAM1,  /ylog, xrange=[0.0,120.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Migrating WNs 3-4 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM4, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.25, 0.88, 'WN 3 Freq 3 ', charsize=1.5
plots, [[0.42,0.89],[0.47,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.25, 0.85, 'WN 4 Freq 4 ', charsize=1.5
plots, [[0.42,0.86],[0.47,0.86]], /NORMAL, linestyle = 1 

plot, tmppWAM5, zplevWAM1,  /ylog, xrange=[0.0,72.0], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional Wind Migrating WNs 5-6 March Mean Equator',$
         xtitle = 'Meridional Wind Phase (Degrees)', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM6, zplevWAM1, linestyle=1

xyouts, /NORMAL, 0.45, 0.88, 'WN 5 Freq 5 ', charsize=1.5
plots, [[0.62,0.89],[0.67,0.89]], /NORMAL, linestyle = 0 
xyouts, /NORMAL, 0.45, 0.85, 'WN 6 Freq 6 ', charsize=1.5
plots, [[0.62,0.86],[0.67,0.86]], /NORMAL, linestyle = 1 

DEVICE, /CLOSE

SET_PLOT, 'X'

END