PRO plot_waccm_wam_recon_lon,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     
;																				    ;
;March Equator Mean Temperature 																    ;March 20N Gravity Wave parameterization High Resolution Temperature
;																				    ;
  SET_PLOT, 'PS'																		    ;  iLt = 58 ;20 North
																				    ;
  device, /color, /land, bits_per_pixel=8,filename='mtm_UT00-UT21_MarchMean_Equator_Tcomp_lon.ps'							    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  LOADCT, 39, ncolors = 10																	    ;  dimsize = size(twzonalmeandev0)
																				    ;  nTAllW = dimsize[3]
  iLt = 47 ;Equator
  
  xTitle = 'Longitude (Degrees)'
  yTitle = 'Temperature (K)'																		    ;
; 																				    ;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHR_20N_zmd_tcomp.ps'
  Title0 = ' Temperature Reconstruction March Mean Equator 285km UT00'										    ;  LOADCT, 39, ncolors = 10
  plot, lon_shift, tn_shift[*,iLt,3,0], xrange=[-90.0,90.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  Title0 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT00'
  oplot, lon_shift, reform(twvtotal1to3_shift[*,47,0])+550., color=8, thick=4 												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon_shift, reform(twvtotal1to6_shift[*,47,0])+550., color=6, thick=4 											   ;  nTAllW = dimsize[3]
  oplot, lon_shift, reform(twvtotal1to9_shift[*,47,0])+550., color=2, thick=4 											   ;  nTAllW = dimsize[3]
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				    ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title0 = ' Temperature Reconstruction March Mean Equator 285km UT00' 													    ;  
  plot, lonWAM_shift, tnWAM_shift[*,46,142,0],xrange=[-90.0,90.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,0])+835., color=8, thick=4											   ;  nTAllW = dimsize[3]
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,0])+835., color=6, thick=4												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,0])+835., color=2, thick=4												   ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title3 = ' Temperature Reconstruction March Mean Equator 285km UT03'													     ;  
  plot, lon_shift, tn_shift[*,iLt,3,1], xrange=[-135.0,45.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon_shift, reform(twvtotal1to3_shift[*,47,1])+550., color=8, thick=4  										 ;  nTAllW = dimsize[3]
  oplot, lon_shift, reform(twvtotal1to6_shift[*,47,1])+550., color=6, thick=4  											 ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon_shift, reform(twvtotal1to9_shift[*,47,1])+550., color=2, thick=4  											 ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title3 = ' Temperature Reconstruction March Mean Equator 285km UT03'													     ;  
  plot, lonWAM_shift, tnWAM_shift[*,iLt,142,3], xrange=[-135.0,45.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,3])+835., color=8, thick=4 											   ;  nTAllW = dimsize[3]
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,3])+835., color=6, thick=4 												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,3])+835., color=2, thick=4 												   ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title6 = ' Temperature Reconstruction March Mean Equator 285km UT06'													     ;  
  plot, lon_shift, tn_shift[*,iLt,3,2],xrange=[-180.0,0.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon_shift, reform(twvtotal1to3_shift[*,47,2])+550., color=8, thick=4  										 ;  nTAllW = dimsize[3]
  oplot, lon_shift, reform(twvtotal1to6_shift[*,47,2])+550., color=6, thick=4  											 ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon_shift, reform(twvtotal1to9_shift[*,47,2])+550., color=2, thick=4  											 ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title6 = ' Temperature Reconstruction March Mean Equator 285km UT06'													     ;  
  plot, lonWAM_shift, tnWAM_shift[*,iLt,142,6], xrange=[-180.0,0.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,6])+835., color=8, thick=4 											   ;  nTAllW = dimsize[3]
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,6])+835., color=6, thick=4 												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,6])+835., color=2, thick=4 												   ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title9 = ' Temperature Reconstruction March Mean Equator 285km UT09'													     ;  
  plot, lon, tn[*,iLt,3,3], xrange=[135.0,315.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon, reform(twvtotal1to3[*,47,3])+550., color=8, thick=4									     ;  nTAllW = dimsize[3]
  oplot, lon, reform(twvtotal1to6[*,47,3])+550., color=6, thick=4										     ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon, reform(twvtotal1to9[*,47,3])+550., color=2, thick=4										     ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title9 = ' Temperature Reconstruction March Mean Equator 285km UT09'													     ;  
  plot, lonWAM, tnWAM[*,iLt,142,9], xrange=[135.0,315.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM, reform(twvtotal1to3WAM[*,46,9])+835., color=8, thick=4 										       ;  nTAllW = dimsize[3]
  oplot, lonWAM, reform(twvtotal1to6WAM[*,46,9])+835., color=6, thick=4 											       ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM, reform(twvtotal1to9WAM[*,46,9])+835., color=2, thick=4 											       ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title12 = ' Temperature Reconstruction March Mean Equator 285km UT12'													     ;  
  plot, lon, tn[*,iLt,3,4], xrange=[90.0,270.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon, reform(twvtotal1to3[*,47,4])+550., color=8, thick=4									     ;  nTAllW = dimsize[3]
  oplot, lon, reform(twvtotal1to6[*,47,4])+550., color=6, thick=4										     ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon, reform(twvtotal1to9[*,47,4])+550., color=2, thick=4										     ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title12 = ' Temperature Reconstruction March Mean Equator 285km UT12'													     ;  
  plot, lonWAM, tnWAM[*,iLt,142,12], xrange=[90.0,270.0],yrange=[750,890],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM, reform(twvtotal1to3WAM[*,46,12])+835., color=8, thick=4										       ;  nTAllW = dimsize[3]
  oplot, lonWAM, reform(twvtotal1to6WAM[*,46,12])+835., color=6, thick=4											       ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM, reform(twvtotal1to9WAM[*,46,12])+835., color=2, thick=4											       ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title15 = ' Temperature Reconstruction March Mean Equator 285km UT15'													     ;  
  plot, lon, tn[*,iLt,3,5], xrange=[45.0,225.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon, reform(twvtotal1to3[*,47,5])+550., color=8, thick=4									     ;  nTAllW = dimsize[3]
  oplot, lon, reform(twvtotal1to6[*,47,5])+550., color=6, thick=4										     ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon, reform(twvtotal1to9[*,47,5])+550., color=2, thick=4										     ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title15 = ' Temperature Reconstruction March Mean Equator 285km UT15'													     ;  
  plot, lonWAM, tnWAM[*,iLt,142,15], xrange=[45.0,225.0],yrange=[760,910],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM, reform(twvtotal1to3WAM[*,46,15])+835., color=8, thick=4										       ;  nTAllW = dimsize[3]
  oplot, lonWAM, reform(twvtotal1to6WAM[*,46,15])+835., color=6, thick=4											       ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM, reform(twvtotal1to9WAM[*,46,15])+835., color=2, thick=4											       ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title18 = ' Temperature Reconstruction March Mean Equator 285km UT18'													     ;  
  plot, lon_shift, tn_shift[*,iLt,3,6], xrange=[0.0,180.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon_shift, reform(twvtotal1to3_shift[*,47,6])+550., color=8, thick=4  										 ;  nTAllW = dimsize[3]
  oplot, lon_shift, reform(twvtotal1to6_shift[*,47,6])+550., color=6, thick=4  											 ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon_shift, reform(twvtotal1to9_shift[*,47,6])+550., color=2, thick=4  											 ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title18 = ' Temperature Reconstruction March Mean Equator 285km UT18'													     ;  
  plot, lonWAM_shift, tnWAM_shift[*,iLt,142,18], xrange=[0.0,180.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,18])+835., color=8, thick=4 											   ;  nTAllW = dimsize[3]
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,18])+835., color=6, thick=4 												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,18])+835., color=2, thick=4 												   ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title21 = ' Temperature Reconstruction March Mean Equator 285km UT21'													     ;  
  plot, lon_shift, tn_shift[*,iLt,3,7], xrange=[-45.0,135.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lon_shift, reform(twvtotal1to3_shift[*,47,7])+550., color=8, thick=4  										 ;  nTAllW = dimsize[3]
  oplot, lon_shift, reform(twvtotal1to6_shift[*,47,7])+550., color=6, thick=4  											 ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lon_shift, reform(twvtotal1to9_shift[*,47,7])+550., color=2, thick=4  											 ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
																				     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
  Title21 = ' Temperature Reconstruction March Mean Equator 285km UT21'													     ;  
  plot, lonWAM_shift, tnWAM_shift[*,iLt,142,21], xrange=[-45.0,135.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
  oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,21])+835., color=8, thick=4 											   ;  nTAllW = dimsize[3]
  oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,21])+835., color=6, thick=4 												   ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
  oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,21])+835., color=2, thick=4 												   ;
  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
  xyouts, /NORMAL, 0.65, 0.9, 'Wavenumbers 1-3 ', color=8, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.86, 'Wavenumbers 1-6 ', color=6, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
  xyouts, /NORMAL, 0.65, 0.82, 'Wavenumbers 1-9 ', color=2, charsize=1.5												    ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4

;  twvtotal1to3 = REFORM(twv[*,*,*,2])
;  twvtotal1to6 = REFORM(twv[*,*,*,2])
;  twvtotal1to9 = REFORM(twv[*,*,*,2])
;
;  FOR iwvn=1,3 DO twvtotal1to3[*,*,*] = twvtotal1to3[*,*,*] + twv[*,*,*,iwvn]
;  FOR iwvn=1,6 DO twvtotal1to6[*,*,*] = twvtotal1to6[*,*,*] + twv[*,*,*,iwvn]
;  FOR iwvn=1,9 DO twvtotal1to9[*,*,*] = twvtotal1to9[*,*,*] + twv[*,*,*,iwvn]
;
;  twvtotal1to3WAM = REFORM(twvWAM[*,*,*,5])
;  twvtotal1to6WAM = REFORM(twvWAM[*,*,*,5])
;  twvtotal1to9WAM = REFORM(twvWAM[*,*,*,5])
;
;  FOR iwvn=1,3 DO twvtotal1to3WAM[*,*,*] = twvtotal1to3WAM[*,*,*] + twvWAM[*,*,*,iwvn]
;  FOR iwvn=1,6 DO twvtotal1to6WAM[*,*,*] = twvtotal1to6WAM[*,*,*] + twvWAM[*,*,*,iwvn]
;  FOR iwvn=1,9 DO twvtotal1to9WAM[*,*,*] = twvtotal1to9WAM[*,*,*] + twvWAM[*,*,*,iwvn]
;
;  twvtotal1to3_shift = SHIFT(twvtotal1to3,72,0,0)
;  twvtotal1to6_shift = SHIFT(twvtotal1to6,72,0,0)
;  twvtotal1to9_shift = SHIFT(twvtotal1to9,72,0,0)
;
;  twvtotal1to3WAM_shift = SHIFT(twvtotal1to3WAM,91,0,0)
;  twvtotal1to6WAM_shift = SHIFT(twvtotal1to6WAM,91,0,0)
;  twvtotal1to9WAM_shift = SHIFT(twvtotal1to9WAM,91,0,0)
;;																				     ;  plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  Title6 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT06'										     ;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt]
;  plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1]), color=8, thick=4
;  oplot, lonw, reform(_shift[*,47,7,2])+550, color=8, thick=4  										     ;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,2])+550., color=6, thick=4 											     ;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav'
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;;  
;;																				     ;
;  xTitleOrig = ' Longitude (Degrees East)'															     ;;  device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT09'										     ;  Title6 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT06'
;  plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]  	     ;  plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,3])+550, color=8, thick=4											     ;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,47,7,3])+550., color=6, thick=4												     ;  oplot, lonw, reform(uwvtotal_shift[*,58,7,2]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;																				     ;;  
;  xTitleOrig = ' Longitude (Degrees East)'															     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon.sav'
;  Title12 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT12'										     ;  dimsize = size(twzonalmeandev0)
;  plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]		     ;  nTAllW = dimsize[3]
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,4])+550, color=8, thick=4											     ;  xTitleOrig = ' Longitude (Degrees East)'
;  oplot, lonworig, reform(uwvtotal[*,47,7,4])+550., color=6, thick=4												     ;;  device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N_T.ps'
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  Title9 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT09'
;;																				     ;  oplot, lonworig, reform(uwvtotal[*,58,7,3]), color=6, thick=4
;  xTitleOrig = ' Longitude (Degrees East)'															     ;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;  Title15 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT15'										     ;;  twzonalmeandev9 = twzonalmeandev
;  plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]		     ;;  twzonalmeandevorig9 = twzonalmeandevorig
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,5])+550, color=8, thick=4											     ;;  
;  oplot, lonworig, reform(uwvtotal[*,47,7,5])+550., color=6, thick=4												     ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon.sav'
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												     ;  dimsize = size(twzonalmeandev0)
;;																				     ;;  device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT18'										     ;  Title12 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT12'
;  plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		     ;  plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,6])+550, color=8, thick=4											     ;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,6])+550., color=6, thick=4 											     ;  oplot, lonworig, reform(uwvtotal[*,58,7,4]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  																				    ;;  
;  Title21 = ' Temperature Deviation from Zonal Mean March Mean Equator 285km UT21'										    ;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon.sav'
;  plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]		    ;  dimsize = size(twzonalmeandev0)
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,7])+550, color=8, thick=4											    ;
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,7])+550., color=6, thick=4												    ;  xTitleOrig = ' Longitude (Degrees East)'
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'												    ;;  device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N_T.ps'
;																				    ;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5]), color=8, thick=4
  device,/close 																		    ;  oplot, lonworig, reform(uwvtotal[*,58,7,5]), color=6, thick=4
  set_plot, 'x' 																		    ;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

;																				    ;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav'
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT18'
;  plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,6]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav'
;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT21'
;  plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,7]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;March Equator Gravity Wave parameterization High Resolution Temperature
;
;  iLt = 47 ;20 North
;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHR_Equator_zmd_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  Title0 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT00'
;  plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,0]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,0]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_Equator_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT03'
;  plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,1]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,1]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_Equator_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT06'
;  plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,2]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,2]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_Equator_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT09'
;  plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,3]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,3]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_Equator_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT12'
;  plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,4]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,4]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_Equator_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March Equator 285km UT15'
;  plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,5]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,5]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_Equator_T.sav'
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_Equator_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT18'
;  plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,6]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,6]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_Equator_T.sav'
;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_31N_FullLon.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_Equator_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT21'
;  plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,7]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,7]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_Equator_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;March 20N 3.1.9 Mean Temperature
;
;  iLt = 58 ;20 North
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_March_20N_zmd_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  Title0 = ' Temperature Deviation from Zonal Mean March 20N 285km UT00'
;  plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_31N_FullLon.sav'
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT03_March_20N_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March 20N 285km UT03'
;  plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,1]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_31N_FullLon_20N_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT06_March_20N_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March 20N 285km UT06'
;  plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,2]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_31N_FullLon_20N_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_March_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March 20N 285km UT09'
;  plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,3]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_31N_FullLon_20N_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_March_20N_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March 20N 285km UT12'
;  plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,4]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_31N_FullLon_20N_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_March_20N_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15'
;  plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,5]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_31N_FullLon_20N_T.sav'
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT18_March_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March 20N 285km UT18'
;  plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,6]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_31N_FullLon_20N_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT21_March_20N_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March 20N 285km UT21'
;  plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,7]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;March Equator 3.1.9 Mean Temperature
;
;  iLt = 47 ;Equator
;;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  SET_PLOT, 'PS'
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_March_Equator_zmd_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  Title0 = ' Temperature Deviation from Zonal Mean March Equator 285km UT00'
;  plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,0]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,0]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_31N_FullLon.sav'
;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT03_March_20N_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March Equator 285km UT03'
;  plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,1]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,1]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT06_March_20N_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March Equator 285km UT06'
;  plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,2]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,2]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_March_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March Equator 285km UT09'
;  plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,3]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,47,7,3]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_March_20N_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March Equator 285km UT12'
;  plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,4]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,47,7,4]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_March_20N_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March Equator 285km UT15'
;  plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,47,7,5]), color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,47,7,5]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_31N_FullLon_Equator_T.sav'
;
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT18_March_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March Equator 285km UT18'
;  plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,6]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,6]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_31N_FullLon_Equator_T.sav'
;;  
;  restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N.sav'
;  dimsize = size(twzonalmeandev0)
;  nTAllW = dimsize[3]
;
;;  device, /color, /land, filename='idl_tmt_UT21_March_20N_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March Equator 285km UT21'
;  plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,7]), color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,47,7,7]), color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;  SET_PLOT, 'PS'
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHRMean_20N_T_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  
;  Title0 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT00'
;  plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
;
;;  device,/close
;;  set_plot, 'x'
;;  
;;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT03'
;  plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,1])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT06'
;  plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,2])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT09'
;  plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,3])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT12'
;  plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,4])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT15'
;  plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,5])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav'
;;
;;  device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT18'
;  plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,6])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT21'
;  plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,7])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;March Equator GWHR Mean Temperature
;
;  iLt = 47 ;Equator
;
;  SET_PLOT, 'PS'
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHRMean_Equator_T_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  
;  Title0 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT00'
;  plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,0])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,0])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_Equator_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT03'
;  plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,1])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,1])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_Equator_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT06'
;  plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,2])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,2])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_Equator_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT09'
;  plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,3])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,3])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_Equator_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT12'
;  plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,4])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,4])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_Equator_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT15'
;  plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,5])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,iLt,7,5])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav'
;;
;;  device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_Equator_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT18'
;  plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,6])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,6])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_Equator_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT21'
;  plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,7])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,iLt,7,7])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;  twzonalmeandevorig21 = twzonalmeandevorig
;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_Equator_T.sav'
;
;  device,/close
;  set_plot, 'x'
;
;  SET_PLOT, 'PS'
;
;  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchMean_20N_T_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  
;  Title0 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT00'
;  plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,0])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav'
;
;;  device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT03'
;  plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1])+530., color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,1])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT06'
;  plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,2])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev6 = twzonalmeandev
;;  twzonalmeandevorig6 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT09'
;  plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,3])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev9 = twzonalmeandev
;;  twzonalmeandevorig9 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT12'
;  plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,4])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev12 = twzonalmeandev
;;  twzonalmeandevorig12 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;  xTitleOrig = ' Longitude (Degrees East)'
;;  device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N_T.ps'
;  Title15 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT15'
;  plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt]
;  oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5])+530, color=8, thick=4
;  oplot, lonworig, reform(uwvtotal[*,58,7,5])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev15 = twzonalmeandev
;;  twzonalmeandevorig15 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav'
;;
;;  device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps'
;  Title18 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT18'
;  plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,6])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev18 = twzonalmeandev
;;  twzonalmeandevorig18 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav'
;;  
;;  device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps'
;  Title21 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT21'
;  plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;;  for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt]
;  oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7])+530, color=8, thick=4
;  oplot, lonw, reform(uwvtotal_shift[*,58,7,7])+530., color=6, thick=4
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;;  twzonalmeandev21 = twzonalmeandev
;;  twzonalmeandevorig21 = twzonalmeandevorig
;;  save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav'
;
;  device,/close
;  set_plot, 'x'
; 

return
end