PRO plot_waccm_wam_t_prof_wvnphase, zplevWACCM, zplevWAM, uphsWACCM, uphsWAM

; 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/twvnphs_March_Mean_Prof_WACCM_WAM.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
;
tmppWACCM1 = REFORM(uphsWACCM[1,47,0:80])
tmppWACCM2 = REFORM(uphsWACCM[2,47,0:80])
tmppWACCM3 = REFORM(uphsWACCM[3,47,0:80])
tmppWACCM4 = REFORM(uphsWACCM[4,47,0:80])
tmppWACCM5 = REFORM(uphsWACCM[5,47,0:80])
tmppWACCM6 = REFORM(uphsWACCM[6,47,0:80])
tmppWACCM7 = REFORM(uphsWACCM[7,47,0:80])
tmppWACCM8 = REFORM(uphsWACCM[8,47,0:80])
tmppWACCM9 = REFORM(uphsWACCM[9,47,0:80])
tmppWACCM10 = REFORM(uphsWACCM[10,47,0:80])
;tmppWACCM2 = REFORM(uphsWACCM[1,37,0:80])
;tmppWACCM3 = REFORM(uphsWACCM[1,24,0:80])
;tmppWACCM4 = REFORM(uphsWACCM[1,29,0:80])

tmppWACCMIn1 = REVERSE(tmppWACCM1)
tmppWACCMIn2 = REVERSE(tmppWACCM2)
tmppWACCMIn3 = REVERSE(tmppWACCM3)
tmppWACCMIn4 = REVERSE(tmppWACCM4)
tmppWACCMIn5 = REVERSE(tmppWACCM5)
tmppWACCMIn6 = REVERSE(tmppWACCM6)
tmppWACCMIn7 = REVERSE(tmppWACCM7)
tmppWACCMIn8 = REVERSE(tmppWACCM8)
tmppWACCMIn9 = REVERSE(tmppWACCM9)
tmppWACCMIn10 = REVERSE(tmppWACCM10)

zplevWACCM1 = REVERSE(zplevWACCM)

tmppWACCM1 = interpol(tmppWACCMIn1, zplevWACCM1, zplevWAM)
tmppWACCM2 = interpol(tmppWACCMIn2, zplevWACCM1, zplevWAM)
tmppWACCM3 = interpol(tmppWACCMIn3, zplevWACCM1, zplevWAM)
tmppWACCM4 = interpol(tmppWACCMIn4, zplevWACCM1, zplevWAM)
tmppWACCM5 = interpol(tmppWACCMIn5, zplevWACCM1, zplevWAM)
tmppWACCM6 = interpol(tmppWACCMIn6, zplevWACCM1, zplevWAM)
tmppWACCM7 = interpol(tmppWACCMIn7, zplevWACCM1, zplevWAM)
tmppWACCM8 = interpol(tmppWACCMIn8, zplevWACCM1, zplevWAM)
tmppWACCM9 = interpol(tmppWACCMIn9, zplevWACCM1, zplevWAM)
tmppWACCM10 = interpol(tmppWACCMIn10, zplevWACCM1, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmppWAM1 = REFORM(uphsWAM[1,46,0:149])
tmppWAM2 = REFORM(uphsWAM[2,46,0:149])
tmppWAM3 = REFORM(uphsWAM[3,46,0:149])
tmppWAM4 = REFORM(uphsWAM[4,46,0:149])
tmppWAM5 = REFORM(uphsWAM[5,46,0:149])
tmppWAM6 = REFORM(uphsWAM[6,46,0:149])
tmppWAM7 = REFORM(uphsWAM[7,46,0:149])
tmppWAM8 = REFORM(uphsWAM[8,46,0:149])
tmppWAM9 = REFORM(uphsWAM[8,46,0:149])
tmppWAM10 = REFORM(uphsWAM[10,46,0:149])
;tmppWAM2 = REFORM(uphsWAM[1,36,0:149])
;tmppWAM3 = REFORM(uphsWAM[1,23,0:149])
;tmppWAM4 = REFORM(uphsWAM[1,28,0:149])

zplevWAM1 = zplevWAM[0:149]

plot, tmppWACCM1, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 1 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmppWAM1, zplevWAM1, linestyle=2

plot, tmppWACCM2, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 2 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmppWAM2, zplevWAM1, linestyle=2

n1 = tmppWAM1 * 0. + 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

plot, tmppWACCM3, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 3 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmppWAM3, zplevWAM1, linestyle=2

n1 = tmppWAM1 * 0. + 60.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 60.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 120.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 120.0
  
oplot, n1, zplevWAM1, linestyle = 3 
  
plot, tmppWACCM4, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 4 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3

oplot, tmppWAM4, zplevWAM1, linestyle=2

n1 = tmppWAM1 * 0. + 45.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 45.0
 
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0.+ 135.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0.- 135.0

oplot, n1, zplevWAM1, linestyle = 3 

plot, tmppWACCM5, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 5 Tide Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM5, zplevWAM1, linestyle=2
 
n1 = tmppWAM1 * 0. + 36.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 36.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 72.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 72.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 108.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 108.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 144.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 144.0
  
oplot, n1, zplevWAM1, linestyle = 3 
  
plot, tmppWACCM6, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 6 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM6, zplevWAM1, linestyle=2
 
n1 = tmppWAM1 * 0. + 30.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 30.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 60.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 60.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 90.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 120.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 120.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 150.0
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 150.0
  
oplot, n1, zplevWAM1, linestyle = 3 
  
plot, tmppWACCM7, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 7 Phase March Equator',$
         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

oplot, tmppWAM7, zplevWAM1, linestyle=2
 
n1 = tmppWAM1 * 0. + 25.7
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 25.7
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 51.4
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 51.4
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 77.1
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 77.1
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 102.8
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 102.8
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 128.5
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 128.5
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. + 154.2
  
oplot, n1, zplevWAM1, linestyle = 3 

n1 = tmppWAM1 * 0. - 154.2
  
oplot, n1, zplevWAM1, linestyle = 3 
  
;plot, tmppWACCM8, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 8 Phase March Equator',$
;         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;         charsize = 1.3 
;
;oplot, tmppWAM8, zplevWAM1, linestyle=2
;
;plot, tmppWACCM9, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 9 Phase March Equator',$
;         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;         charsize = 1.3 
;
;oplot, tmppWAM9, zplevWAM1, linestyle=2
;
;plot, tmppWACCM10, zplevWAM1,  /ylog, xrange=[-180,180], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Temperature Wavenumber 10 Phase March Equator',$
;         xtitle = 'Temperature Phase', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;         charsize = 1.3 
;
;oplot, tmppWAM10, zplevWAM1, linestyle=2

DEVICE, /CLOSE

SET_PLOT, 'x'

END