print, 'Running plot_zm' 

  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

;    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 = ' Latitude (Degrees)'
  yTitle = 'Temperature Deviation from Zonal Mean (hPa)'
  !P.CHARSIZE = 1.5
  
  pTop = 1.E-10
  pBot = 1000.0
  
  set_plot, 'ps'
  device, /color, /land, filename='idl_mtm_waccm_UT00-21_March_ZM.ps'
  loadct, 39
  
  Title = ' WACCM-X Temperature Zonal Mean March Mean UT00 '
  CONTOUR, REFORM(tNWZM[*,*,0]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,0]),lat, 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 Zonal Mean March Mean UT03 '
  CONTOUR, REFORM(tNWZM[*,*,1]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,1]),lat, 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 Zonal Mean March Mean UT06 '
  CONTOUR, REFORM(tNWZM[*,*,2]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,2]),lat, 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 Zonal Mean March Mean UT09 '
  CONTOUR, REFORM(tNWZM[*,*,3]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,3]),lat, 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 Zonal Mean March Mean UT12 '
  CONTOUR, REFORM(tNWZM[*,*,4]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,4]),lat, 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 Zonal Mean March Mean UT15 '
  CONTOUR, REFORM(tNWZM[*,*,5]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,5]),lat, 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 Zonal Mean March Mean UT18 '
  CONTOUR, REFORM(tNWZM[*,*,6]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,6]),lat, 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 Zonal Mean March Mean UT21 '
  CONTOUR, REFORM(tNWZM[*,*,7]),lat, lev, nlev=24, /fill,c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0],yrange=[pBot,pTop], xstyle=1,ystyle=1,/ylog,Title=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
  CONTOUR, REFORM(tNWZM[*,*,7]),lat, lev, nlev=24,/follow,xrange=[-90.0,90.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