PRO plot_waccm_wam_recon_lon_15n_jun,zgm,tn,tnWAM,lon,lonWAM,twv1ZAllT,twv1ZAllTWAM

; 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

  LTimes = ['18','20','22','0','2','4','6']

  twvtotal1to3 = REFORM(twv1ZAllT[*,*,*,1])
  twvtotal1to6 = REFORM(twv1ZAllT[*,*,*,1])
  twvtotal1to9 = REFORM(twv1ZAllT[*,*,*,1])

  FOR iwvn=2,3 DO twvtotal1to3[*,*,*] = twvtotal1to3[*,*,*] + twv1ZAllT[*,*,*,iwvn]
  FOR iwvn=2,6 DO twvtotal1to6[*,*,*] = twvtotal1to6[*,*,*] + twv1ZAllT[*,*,*,iwvn]
  FOR iwvn=2,9 DO twvtotal1to9[*,*,*] = twvtotal1to9[*,*,*] + twv1ZAllT[*,*,*,iwvn]

  twvtotal1to3_shift = SHIFT(twvtotal1to3,72,0,0)
  twvtotal1to6_shift = SHIFT(twvtotal1to6,72,0,0)
  twvtotal1to9_shift = SHIFT(twvtotal1to9,72,0,0)

  twvtotal1to3WAM = REFORM(twv1ZAllTWAM[*,*,*,1])
  twvtotal1to6WAM = REFORM(twv1ZAllTWAM[*,*,*,1])
  twvtotal1to9WAM = REFORM(twv1ZAllTWAM[*,*,*,1])

  FOR iwvn=2,3 DO twvtotal1to3WAM[*,*,*] = twvtotal1to3WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn]
  FOR iwvn=2,6 DO twvtotal1to6WAM[*,*,*] = twvtotal1to6WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn]
  FOR iwvn=2,9 DO twvtotal1to9WAM[*,*,*] = twvtotal1to9WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn]

  twvtotal1to3WAM_shift = SHIFT(twvtotal1to3WAM,96,0,0)
  twvtotal1to6WAM_shift = SHIFT(twvtotal1to6WAM,96,0,0)
  twvtotal1to9WAM_shift = SHIFT(twvtotal1to9WAM,96,0,0)
  
  tn_shift = SHIFT(tn,72,0,0,0)

  tnWAM_shift = SHIFT(tnWAM,96,0,0,0)    
  
  lon_shift = lon - 180.0

  lonWAM_shift = lonWAM - 180.0  
																				    ;  nTAllW = dimsize[3]
  iLt = 56 ;15 North
  iLtWAM = 55 ;15 North WAM
  
  dimszgm = size(zgm)
  nx = dimszgm(1)
  ny = dimszgm(2)
  nT = dimszgm(4)

  tNZSub = FLTARR(nX,nT)
  tNZSub_shift = FLTARR(nX,nT)
  twvtotal1to3Sub = FLTARR(nX,nT)
  twvtotal1to3Sub_shift = FLTARR(nX,nT)
  twvtotal1to6Sub = FLTARR(nX,nT)
  twvtotal1to6Sub_shift = FLTARR(nX,nT)
  twvtotal1to9Sub = FLTARR(nX,nT)
  twvtotal1to9Sub_shift = FLTARR(nX,nT)
;
;  Find altitude index for this time
;    
  FOR iLn=0,nX-1 DO BEGIN

    FOR iT = 0, nT-1 DO BEGIN

      zWLnLtT = zgm[iLn,iLt,*,iT]
      zInd = WHERE(zWLnLtT GE 275.0 AND zWLnLtT LE 295.0)

      IF zInd[0] EQ -1 THEN BEGIN
  
        MESSAGE, 'Cannot find altitude index in WACCM'
     
      ENDIF ELSE BEGIN
  
        iZ = zInd[0]
     
      ENDELSE

      tNZSub[iLn,iT] = tN[iLn,iLt,iZ,iT]
      tNZSub_shift[iLn,iT] = tN_shift[iLn,iLt,iZ,iT]
      twvtotal1to3Sub[iLn,iT] = twvtotal1to3[iLn,iLt,iT]
      twvtotal1to3Sub_shift[iLn,iT] = twvtotal1to3_shift[iLn,iLt,iT]
      twvtotal1to6Sub[iLn,iT] = twvtotal1to6[iLn,iLt,iT]
      twvtotal1to6Sub_shift[iLn,iT] = twvtotal1to6_shift[iLn,iLt,iT]
      twvtotal1to9Sub[iLn,iT] = twvtotal1to9[iLn,iLt,iT]
      twvtotal1to9Sub_shift[iLn,iT] = twvtotal1to9_shift[iLn,iLt,iT]

    ENDFOR

  ENDFOR
;																				    ;
;June Equator Mean Temperature 															
;																		
  SET_PLOT, 'PS'																
																		
  device, /color, /land, bits_per_pixel=8,filename='mtm_UT00-UT23_JuneMean_15N_Tcomp_lon.ps'						
  LOADCT, 39, ncolors = 10															
  
  xTitle = 'Longitude (Degrees)'
  yTitle = 'Temperature (K)'															
; 																		
  Title0 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT00'									
  plot, lon_shift, tNZSub_shift[*,0], xrange=[-90.0,90.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,0])+555., color=8, thick=4									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,0])+555., color=6, thick=4									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,0])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title0 = ' WAM Temperature Reconstruction June Mean Equator 285km UT00' 											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,0],xrange=[-90.0,90.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,0])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,0])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,0])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												  
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												  
																				  
  Title1 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT01'									
  plot, lon_shift, tNZSub_shift[*,1], xrange=[-105.0,75.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title1,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,1])+555., color=8, thick=4									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,1])+555., color=6, thick=4									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,1])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title1 = ' WAM Temperature Reconstruction June Mean Equator 285km UT01' 											
  plot, lonWAM_shift,tnWAM_shift[*,iLtWAM,142,1],xrange=[-105.0,75.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title1,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,1])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,1])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,1])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												  
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												  
																				  
  Title2 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT02'									
  plot, lon_shift, tNZSub_shift[*,2], xrange=[-120.0,60.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title2,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,2])+555., color=8, thick=4									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,2])+555., color=6, thick=4									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,2])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												  
																				  
  Title2 = ' WAM Temperature Reconstruction June Mean Equator 285km UT02' 											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,2],xrange=[-120.0,60.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title2,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,2])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,2])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,2])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title3 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT03'										
  plot, lon_shift, tNZSub_shift[*,3], xrange=[-135.0,45.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,3])+555., color=8, thick=4									       
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,3])+555., color=6, thick=4										 
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,3])+555., color=2, thick=4										 
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title3 = ' WAM Temperature Reconstruction June Mean Equator 285km UT03'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,3], xrange=[-135.0,45.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,3])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,3])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,3])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												   
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												   
																				   
  Title4 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT04'											
  plot, lon_shift, tNZSub_shift[*,4], xrange=[-150.0,30.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title4,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		   
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,4])+555., color=8, thick=4									      
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,4])+555., color=6, thick=4									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,4])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title4 = ' WAM Temperature Reconstruction June Mean Equator 285km UT04'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,4], xrange=[-150.0,30.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title4,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,4])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,4])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,4])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title5 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT05'									
  plot, lon_shift, tNZSub_shift[*,5], xrange=[-165.0,15.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title5,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,5])+555., color=8, thick=4									      
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,5])+555., color=6, thick=4									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,5])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title5 = ' WAM Temperature Reconstruction June Mean Equator 285km UT05'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,5], xrange=[-165.0,15.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title5,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,5])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,5])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,5])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												    
																				    
  Title6 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT06'											
  plot, lon_shift, tNZSub_shift[*,6],xrange=[-180.0,0.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,6])+555., color=8, thick=4  									     
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,6])+555., color=6, thick=4  									     
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,6])+555., color=2, thick=4  									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												    
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
																				    
  Title6 = ' WAM Temperature Reconstruction June Mean Equator 285km UT06'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,6], xrange=[-180.0,0.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,6])+935., color=8, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,6])+935., color=6, thick=4 									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,6])+935., color=2, thick=4 									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												    
																				    
  Title7 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT07'											
  plot, lon, tNZSub[*,7],xrange=[165.0,345.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title7,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]   
  oplot, lon, reform(twvtotal1to3Sub[*,7])+555., color=8, thick=4									     
  oplot, lon, reform(twvtotal1to6Sub[*,7])+555., color=6, thick=4									     
  oplot, lon, reform(twvtotal1to9Sub[*,7])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												    
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
																				    
  Title7 = ' WAM Temperature Reconstruction June Mean Equator 285km UT07'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,7], xrange=[165.0,345.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title7,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,7])+935., color=8, thick=4 								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,7])+935., color=6, thick=4 								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,7])+935., color=2, thick=4 								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title8 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT08'								
  plot, lon, tNZSub[*,8],xrange=[150.0,330.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title8,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,8])+555., color=8, thick=4									   
  oplot, lon, reform(twvtotal1to6Sub[*,8])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,8])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
																	
  Title8 = ' WAM Temperature Reconstruction June Mean Equator 285km UT08'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,8], xrange=[150.0,330.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title8,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,8])+935., color=8, thick=4									
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,8])+935., color=6, thick=4									
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,8])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title9 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT09'									
  plot, lon, tNZSub[*,9], xrange=[135.0,315.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lon, reform(twvtotal1to3Sub[*,9])+555., color=8, thick=4									       
  oplot, lon, reform(twvtotal1to6Sub[*,9])+555., color=6, thick=4										
  oplot, lon, reform(twvtotal1to9Sub[*,9])+555., color=2, thick=4										
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title9 = ' WAM Temperature Reconstruction June Mean Equator 285km UT09'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,9], xrange=[135.0,315.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,9])+935., color=8, thick=4									
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,9])+935., color=6, thick=4									
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,9])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
																		
  Title10 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT10'									
  plot, lon, tNZSub[*,10], xrange=[120.0,300.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title10,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lon, reform(twvtotal1to3Sub[*,10])+555., color=8, thick=4									       
  oplot, lon, reform(twvtotal1to6Sub[*,10])+555., color=6, thick=4										
  oplot, lon, reform(twvtotal1to9Sub[*,10])+555., color=2, thick=4										
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title10 = ' WAM Temperature Reconstruction June Mean Equator 285km UT10'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,10], xrange=[120.0,300.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title10,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,10])+935., color=8, thick=4									
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,10])+935., color=6, thick=4									
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,10])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
																		
  Title11 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT11'									
  plot, lon, tNZSub[*,11], xrange=[105.0,285.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title11,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lon, reform(twvtotal1to3Sub[*,11])+555., color=8, thick=4									     
  oplot, lon, reform(twvtotal1to6Sub[*,11])+555., color=6, thick=4										
  oplot, lon, reform(twvtotal1to9Sub[*,11])+555., color=2, thick=4										
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
																		
  Title11 = ' WAM Temperature Reconstruction June Mean Equator 285km UT11'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,11], xrange=[105.0,285.0],yrange=[840,990],xstyle=9,ystyle=1,subtitle=Title11,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,11])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,11])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,11])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
																	
  Title12 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT12'								
  plot, lon, tNZSub[*,12], xrange=[90.0,270.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,12])+555., color=8, thick=4									    
  oplot, lon, reform(twvtotal1to6Sub[*,12])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,12])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title12 = ' WAM Temperature Reconstruction June Mean Equator 285km UT12'								
  plot, lonWAM, tnWAM[*,iLtWAM,142,12], xrange=[90.0,270.0],yrange=[850,990],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,12])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,12])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,12])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title13 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT13'								
  plot, lon, tNZSub[*,13], xrange=[75.0,255.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title13,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,13])+555., color=8, thick=4									    
  oplot, lon, reform(twvtotal1to6Sub[*,13])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,13])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title13 = ' WAM Temperature Reconstruction June Mean Equator 285km UT13'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,13], xrange=[75.0,255.0],yrange=[850,990],xstyle=9,ystyle=1,subtitle=Title13,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,13])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,13])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,13])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									

  Title14 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT14'								
  plot, lon, tNZSub[*,14], xrange=[60.0,240.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title14,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,14])+555., color=8, thick=4									    
  oplot, lon, reform(twvtotal1to6Sub[*,14])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,14])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title14 = ' WAM Temperature Reconstruction June Mean Equator 285km UT14'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,14], xrange=[60.0,240.0],yrange=[850,990],xstyle=9,ystyle=1,subtitle=Title14,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,14])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,14])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,14])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title15 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT15'								
  plot, lon, tNZSub[*,15], xrange=[45.0,225.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,15])+555., color=8, thick=4									    
  oplot, lon, reform(twvtotal1to6Sub[*,15])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,15])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title15 = ' WAM Temperature Reconstruction June Mean Equator 285km UT15'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,15], xrange=[45.0,225.0],yrange=[860,1010],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,15])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,15])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,15])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
																	
  Title16 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT16'								
  plot, lon, tNZSub[*,16], xrange=[30.0,210.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title16,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,16])+555., color=8, thick=4									     
  oplot, lon, reform(twvtotal1to6Sub[*,16])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,16])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title16 = ' WAM Temperature Reconstruction June Mean Equator 285km UT16'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,16], xrange=[30.0,210.0],yrange=[860,1010],xstyle=9,ystyle=1,subtitle=Title16,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,16])+935., color=8, thick=4								
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,16])+935., color=6, thick=4								
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,16])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
																	
  Title17 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT17'								
  plot, lon, tNZSub[*,17], xrange=[15.0,195.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title17,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon, reform(twvtotal1to3Sub[*,17])+555., color=8, thick=4									    
  oplot, lon, reform(twvtotal1to6Sub[*,17])+555., color=6, thick=4									
  oplot, lon, reform(twvtotal1to9Sub[*,17])+555., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'									
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5								
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5									
																	
  Title17 = ' WAM Temperature Reconstruction June Mean Equator 285km UT17'											
  plot, lonWAM, tnWAM[*,iLtWAM,142,17], xrange=[15.0,195.0],yrange=[860,1010],xstyle=9,ystyle=1,subtitle=Title17,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lonWAM, reform(twvtotal1to3WAM[*,iLtWAM,17])+935., color=8, thick=4									
  oplot, lonWAM, reform(twvtotal1to6WAM[*,iLtWAM,17])+935., color=6, thick=4									
  oplot, lonWAM, reform(twvtotal1to9WAM[*,iLtWAM,17])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
																		
  Title18 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT18'									
  plot, lon_shift, tNZSub_shift[*,18], xrange=[0.0,180.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,18])+555., color=8, thick=4  								
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,18])+555., color=6, thick=4  								
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,18])+555., color=2, thick=4  								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'										
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5									
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5									
																		
  Title18 = ' WAM Temperature Reconstruction June Mean Equator 285km UT18'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,18], xrange=[0.0,180.0],yrange=[850,1000],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,18])+935., color=8, thick=4								
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,18])+935., color=6, thick=4								
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,18])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title19 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT19'										
  plot, lon_shift, tNZSub_shift[*,19], xrange=[-15.0,165.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title19,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,19])+555., color=8, thick=4  									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,19])+555., color=6, thick=4  									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,19])+555., color=2, thick=4  									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
																			
  Title19 = ' WAM Temperature Reconstruction June Mean Equator 285km UT19'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,19], xrange=[-15.0,165.0],yrange=[850,1000],xstyle=9,ystyle=1,subtitle=Title19,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,19])+935., color=8, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,19])+935., color=6, thick=4								
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,19])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title20 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT20'										
  plot,  lon_shift, tNZSub_shift[*,20], xrange=[-30.0,150.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title20,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,20])+555., color=8, thick=4  									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,20])+555., color=6, thick=4  									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,20])+555., color=2, thick=4  									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
																			
  Title20 = ' WAM Temperature Reconstruction June Mean Equator 285km UT20'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,20], xrange=[-30.0,150.0],yrange=[850,1000],xstyle=9,ystyle=1,subtitle=Title20,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,20])+935., color=8, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,20])+935., color=6, thick=4								
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,20])+935., color=2, thick=4								
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title21 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT21'										
  plot, lon_shift, tNZSub_shift[*,21], xrange=[-45.0,135.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,21])+555., color=8, thick=4  									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,21])+555., color=6, thick=4  									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,21])+555., color=2, thick=4  									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title21 = ' WAM Temperature Reconstruction June Mean Equator 285km UT21'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,21], xrange=[-45.0,135.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,21])+935., color=8, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,21])+935., color=6, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,21])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
																			
  Title22 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT22'										
  plot, lon_shift, tNZSub_shift[*,22], xrange=[-60.0,120.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title22,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,22])+555., color=8, thick=4  									
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,22])+555., color=6, thick=4  									
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,22])+555., color=2, thick=4  									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'											
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5										
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5											
																			
  Title22 = ' WAM Temperature Reconstruction June Mean Equator 285km UT22'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,22], xrange=[-60.0,120.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title22,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,22])+935., color=8, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,22])+935., color=6, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,22])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												
																				
  Title23 = ' WACCM-X Temperature Reconstruction June Mean Equator 285km UT23'											
  plot, lon_shift, tNZSub_shift[*,23], xrange=[-75.0,105.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title23,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]	
  oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,23])+555., color=8, thick=4  										
  oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,23])+555., color=6, thick=4  										
  oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,23])+555., color=2, thick=4  										
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												
																				
  Title23 = ' WAM Temperature Reconstruction June Mean Equator 285km UT23'											
  plot, lonWAM_shift, tnWAM_shift[*,iLtWAM,142,23], xrange=[-75.0,105.0],yrange=[860,1000],xstyle=9,ystyle=1,subtitle=Title23,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,iLtWAM,23])+935., color=8, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,iLtWAM,23])+935., color=6, thick=4									
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,iLtWAM,23])+935., color=2, thick=4									
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												  
  xyouts, /NORMAL, 0.65, 0.9, 'Zonal Wavenumbers 1-3 ', color=8, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.86, 'Zonal Wavenumbers 1-6 ', color=6, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.82, 'Zonal Wavenumbers 1-9 ', color=2, charsize=1.5											
  xyouts, /NORMAL, 0.65, 0.78, 'All Zonal Wavenumbers ', charsize=1.5												  

  device,/close 																		  
  set_plot, 'x' 																		 

return
end