PRO plotting_recon_wam, tn, lat, lon, uspc, compd_sd_td_qd, uwvtotal

; 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 = ['9','11','13','15','17','19','21','23','1','3','5','7','9']
  LTimes = ['18','20','22','0','2','4','6']

  dims = size(tn)
  nx = dims[1]
  ny = dims[2]
  nt = dims[4]

  compd_sd_td_qdtemp = FLTARR(nx,ny,nt)
  uwvtotaltemp = FLTARR(nx,ny,nt)
  tntemp = FLTARR(nx,ny,nt)

  compd_sd_td_qdtemp[*,*,*] = REFORM(compd_sd_td_qd[*,*,136,*])
  uwvtotaltemp[*,*,*] = REFORM(uwvtotal[*,*,136,*])
  tntemp[*,*,*] = REFORM(tn[*,*,136,*])

  compd_sd_td_qd_shift = SHIFT(compd_sd_td_qdtemp,91,0,0)
  uwvtotal_shift = SHIFT(uwvtotaltemp,91,0,0)
  tn_shift = SHIFT(tntemp,91,0,0)
  
  lonworig = lon
  lonw = lon - 180.0 

;  SET_PLOT, 'PS'

;device, /color, /land, filename='idl_mtm_UT00_MarchGWHR_20N_zmd_tcomp.ps'
;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8
;DEVICE, /COLOR, /LAND, FILENAME = psfile
;LOADCT, 39, ncolors = 10

;psfile = '/home/joemci/MackStuff/MyVersion/tcomp_March_21.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_21_MSDTT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MSTDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MSTDT1.ps'

;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_MeanGWHR_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_Mean_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MDT1.ps'

;
;March 20N Gravity Wave parameterization High Resolution Temperature
;
;  iLt = 58 ;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_20N_zmd_tcomp.ps'
;  LOADCT, 39, ncolors = 10
;  Title0 = ' Temperature Deviation from Zonal Mean March GWHR 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_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_20N_T.ps'
;  Title3 = ' Temperature Deviation from Zonal Mean March GWHR 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_MarchGWHR_31N_FullLon_20N_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_20N_T.ps'
;  Title6 = ' Temperature Deviation from Zonal Mean March GWHR 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_MarchGWHR_31N_FullLon_20N_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_20N_T.ps'
;  Title9 = ' Temperature Deviation from Zonal Mean March GWHR 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_MarchGWHR_31N_FullLon_20N_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_20N_T.ps'
;  Title12 = ' Temperature Deviation from Zonal Mean March GWHR 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_MarchGWHR_31N_FullLon_20N_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_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_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 = 46 ;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'
; 
;
;March Equator Mean Temperature
;
  SET_PLOT, 'PS'

  device, /color, /land, bits_per_pixel=8, filename='idl_mtm_wam_UT00-UT21_MarchMean_Equator_T_tcomp.ps'
  LOADCT, 39, ncolors = 10

  iLt = 46 ;Equator
; 
  Title0 = ' WAM Temperature March Mean Equator ~285km UT00'
  plot, lonw, tn_shift[*,iLt,0], xrange=[-90.0,90.0],yrange=[700,900],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[*,46,0])+825., color=8, thick=4
  oplot, lonw, reform(uwvtotal_shift[*,46,0])+825., 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'

  Title3 = ' WAM Temperature March Mean Equator ~285km UT03'
  plot, lonw, tn_shift[*,iLt,3], xrange=[-135.0,45.0],yrange=[700,900],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[*,46,3])+825., color=8, thick=4
  oplot, lonw, reform(uwvtotal_shift[*,46,3])+825., 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'
;  
  Title6 = ' WAM Temperature March Mean Equator ~285km UT06'
  plot, lonw, tn_shift[*,iLt,6], xrange=[-180.0,0.0],yrange=[700,900],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[*,46,6])+825, color=8, thick=4
  oplot, lonw, reform(uwvtotal_shift[*,46,6])+825., 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'
;  
  xTitleOrig = ' Longitude (Degrees East)'
  Title9 = ' WAM Temperature March Mean Equator ~285km UT09'
  plot, lonworig,tn[*,iLt,136,9],xrange=[135.0,315.0],yrange=[700,900],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[*,46,136,9])+825, color=8, thick=4
  oplot, lonworig, reform(uwvtotal[*,46,136,9])+825., 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'
;  
  xTitleOrig = ' Longitude (Degrees East)'
  Title12 = ' WAM Temperature March Mean Equator ~285km UT12'
  plot, lonworig,tn[*,iLt,136,12],xrange=[90.0,270.0],yrange=[700,900],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[*,46,136,12])+825, color=8, thick=4
  oplot, lonworig, reform(uwvtotal[*,46,136,12])+825., 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'
;  
  xTitleOrig = ' Longitude (Degrees East)'
  Title15 = ' WAM Temperature March Mean Equator ~285km UT15'
  plot,lonworig,tn[*,iLt,136,15],xrange=[45.0,225.0],yrange=[700,900],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[*,46,136,15])+825, color=8, thick=4
  oplot, lonworig, reform(uwvtotal[*,46,136,15])+825., 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'
;
  Title18 = ' WAM Temperature March Mean Equator ~285km UT18'
  plot, lonw, tn_shift[*,iLt,18], xrange=[0.0,180.0],yrange=[700,900],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[*,46,18])+825, color=8, thick=4
  oplot, lonw, reform(uwvtotal_shift[*,46,18])+825., 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'
;  
  Title21 = ' WAM Temperature March Mean Equator ~285km UT21'
  plot, lonw, tn_shift[*,iLt,21], xrange=[-45.0,135.0],yrange=[700,900],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[*,46,21])+825, color=8, thick=4
  oplot, lonw, reform(uwvtotal_shift[*,46,21])+825., 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'

return
end