PRO plot_recon_waccm,tnrecon,lat,lon,lev,time,nx,ny,nz,nt

  print, 'Running plot_recon_waccm' 

;  Id = NCDF_OPEN('/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03_MeanTZ3.nc', /NOWRITE)
;
;; Coordinates
;
;  ncdf_varget, Id, 'lon', lon         ; longitude
;  ncdf_varget, Id, 'lat', lat         ; latitude
;  ncdf_varget, Id, 'lev', lev         ; log pressure, -ln(5.e-7/p) (p in hPa)
;  ncdf_varget, Id, 'time', time       ; time in minutes
;
;  nX = n_elements(lon)
;  nY = n_elements(lat)
;  nZ = n_elements(lev)
;  nT = n_elements(time)
;
;  ncdf_varget, Id, 'T', tN      ;neutral temperature
;
;  NCDF_CLOSE, Id
;
   iLat = 47
;
;    lonOrig = lon
;
;  Convert longitude from 0-360 to -180-180
;
;    lon = lon - 180.0
;
;  print, lon
;
;  tNWZM = FLTARR(nY,nZ,nT)
;
;  FOR iLt=0,nY-1 DO BEGIN
;
;  Find zonal mean for entire longitude range
;
;    FOR iT=0,nT-1 DO BEGIN
;
;      FOR iZ=0,nZ-1 DO BEGIN
;   
;        tNWZM[iLt,iZ,iT] = MEAN(tN[*,iLt,iZ,iT])
;       
;      print, 'iLt, iT, twzonalmean[iLt,iT] ', iLt, iT, twzonalmean[iLt,iT]
;
;      ENDFOR
;
;    ENDFOR
;
;  ENDFOR ;Latitude
;
;  nLonTemp = nXWAM+nXWAM/2
;  tNWAMZMTemp = FLTARR(nLonTemp,nYWAM,nTWAM)
;  tNWAMZMOrig[*,*,*] = tNWAMZM[*,*,*]
;  tNWAMZMTemp[0:nXWAM-1,*,*] = tNWAMZM[0:nXWAM-1,*,*]
;  tNWAMZMTemp[nXWAM:nLonTemp-1,*,*] = tNWAMZM[0:nXWAM/2-1,*,*]
;  tNWAMZM[0:nXWAM-1,*,*] = tNWAMZMTemp[nXWAM/2:nLonTemp-1,*,*]

;
;March 20N High Resolution Temperature Deviation from Zonal Mean
;
;  xTitle = ' Longitude (Degrees)'
;  yTitle = 'Temperature Diurnal Tide Mean (K)'
;  !P.CHARSIZE = 1.5
;  
;  pTop = 1.E-10
;  pBot = 1000.0
;  
;  set_plot, 'ps'
;  device, /color, /land, filename='idl_mtm_waccm_UT00-21_March_Equator_Diurnal.ps'
;  loadct, 39
;  
;  Title = ' WACCM-X Temperature Diurnal Tide March Month Mean UT00 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT03 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT06 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT09 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT12 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT15 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT18 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT21 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  device,/close
;  set_plot, 'x'
;
;
;March 20N High Resolution Temperature Deviation from Zonal Mean
;
;  xTitle = ' Longitude (Degrees)'
;  yTitle = 'Temperature Diurnal Tide Mean (K)'
;  !P.CHARSIZE = 1.5
;  
;  pTop = 1.E-4
;  pBot = 1000.0
;  
;  set_plot, 'ps'
;  device, /color, /land, filename='idl_mtm_waccm_UT00-21_March_Equator_Diurnal_lower.ps'
;  loadct, 39
;  
;  Title = ' WACCM-X Temperature Diurnal Tide March Month Mean UT00 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT03 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT06 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT09 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT12 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT15 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT18 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT21 '
;  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2
;
;  device,/close
;  set_plot, 'x'
;
;March 20N High Resolution Temperature Deviation from Zonal Mean
;
  xTitle = ' Longitude (Degrees)'
  yTitle = 'Temperature Semi-Diurnal Tide Mean (K)'
  !P.CHARSIZE = 1.5
  
  pTop = 1.E-10
  pBot = 1000.0
  
  set_plot, 'ps'
  device, /color, /land, filename='idl_mtm_waccm_UT00-21_March_Equator_Semi-Diurnal.ps'
  loadct, 39
  
  Title = ' WACCM-X Temperature Diurnal Tide March Month Mean UT00 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT03 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT06 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT09 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT12 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT15 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT18 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT21 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  device,/close
  set_plot, 'x'

;
;March 20N High Resolution Temperature Deviation from Zonal Mean
;
  xTitle = ' Longitude (Degrees)'
  yTitle = 'Temperature Diurnal Tide Mean (K)'
  !P.CHARSIZE = 1.5
  
  pTop = 1.E-4
  pBot = 1000.0
  
  set_plot, 'ps'
  device, /color, /land, filename='idl_mtm_waccm_UT00-21_March_Equator_Semi-Diurnal_lower.ps'
  loadct, 39
  
  Title = ' WACCM-X Temperature Diurnal Tide March Month Mean UT00 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,0]),lon, lev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT03 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,1]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT06 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,2]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT09 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,3]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT12 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,4]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT15 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,5]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT18 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,6]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WACCM-X Temperature Diurnal Tide March Mean UT21 '
  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[0.0,360.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNRecon[*,iLat,*,7]),lon, lev, nlev=24,/follow,xrange=[0.0,360.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  device,/close
  set_plot, 'x'

;  Title0 = ' WAM Temperature Deviation from Zonal Mean HR March 20N 285km UT00'
;  plot, lonWAM, tNWAMZM[*,*,0], xrange=[-90.0,90.0],yrange=[1000.0,1.E-09],xstyle=9,ystyle=1,/ylog,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;  axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'
;  save, lonw, lonworig, latw, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_March_20N_HR_FullLon.sav'



;tZSubWAMOrig = tZSubWOrig
;tZSubWAM = tZSubW     
;lonWAMOrig = lonWOrig

;SAVE,tZSubWAM,tZSubWAMOrig,latWAM,lonWAMOrig,presWAM,timeWAM,nXWAM,nYWAM,nZWAM,nTWAM, file = 'wam_mtm_UT0-24_Mar_P7E-08.sav'
  

END