print, 'Running plot_zm_u' 

  Id = NCDF_OPEN('/net/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03_MeanUV.nc', /NOWRITE)
  IdUWAM = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WAM/wam_u_mar.nc', /NOWRITE)

;; Coordinates

  ncdf_varget, Id, 'lon', lon         ; longitude
  ncdf_varget, Id, 'lat', lat         ; latitude
  ncdf_varget, Id, 'lev', zplev         ; log pressure, -ln(5.e-7/p) (p in hPa)
  ncdf_varget, Id, 'time', time_min       ; time in minutes

  ncdf_varget, IdUWAM, 'lon', lonWAM         ; longitude
  ncdf_varget, IdUWAM, 'lat', latWAM         ; latitude
  ncdf_varget, IdUWAM, 'lev', zplevWAM       ; log pressure, -ln(5.e-7/p) (p in hPa)
  ncdf_varget, IdUWAM, 'time', time_minWAM   ; time in minutes

  nX = n_elements(lon)
  nY = n_elements(lat)
  nZ = n_elements(zplev)
  nT = n_elements(time_min)

  nxWAM = n_elements(lonWAM)
  nyWAM = n_elements(latWAM)
  nzWAM = n_elements(zplevWAM)
  nTWAM = n_elements(time_minWAM)

; Dimensions lon,lat,lev,time

  ncdf_varget, Id, 'U', uN      ;neutral zonal wind

  NCDF_CLOSE, Id

; Dimensions lon,lat,lev,time

  ncdf_varget, IdUWAM, 'u', unWAM      ;neutral zonal wind

  NCDF_CLOSE, IdUWAM

  uNWZM = 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
   
        uNWZM[iLt,iZ,iT] = MEAN(uN[*,iLt,iZ,iT])
       
;      print, 'iLt, iT, uwzonalmean[iLt,iT] ', iLt, iT, uwzonalmean[iLt,iT]

      ENDFOR

    ENDFOR

  ENDFOR ;Latitude

  uNWZMR = REVERSE(uNWZM,2)
  
  zplevWACCM = REVERSE(zplev)

  uNWZMInt = FLTARR(nY,nZWAM,nT)

  FOR iLt=0,nY-1 DO BEGIN
    FOR iT=0,nT-1 DO BEGIN
      uNWZMInt[iLT,*,iT] = interpol(REFORM(uNWZMR[iLt,*,iT]), zplevWACCM, zplevWAM)
    ENDFOR
  ENDFOR
;
;  Do WAM zonal mean
;
  uNZMWAM = FLTARR(nYWAM,nZWAM,nTWAM)

  FOR iLt=0,nYWAM-1 DO BEGIN
;
;  Find zonal mean for entire longitude range
;
    FOR iT=0,nTWAM-1 DO BEGIN

      FOR iZ=0,nZWAM-1 DO BEGIN
   
        uNZMWAM[iLt,iZ,iT] = MEAN(uNWAM[*,iLt,iZ,iT])

      ENDFOR

    ENDFOR

  ENDFOR ;Latitude

;
;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='u_waccm_wam_UT00_March_ZM.ps'
  loadct, 39

  plot, REFORM(uNWZMInt[47,*,0]), zplevWAM,  /ylog, xrange=[-60,55], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Zonal Wind Zonal Mean March Mean UT00 Equator',$
         xtitle = 'Zonal Wind ', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

  oplot, REFORM(uNZMWAM[46,*,0]), zplevWAM, linestyle=2

  plot, REFORM(uNWZMInt[24,*,0]), zplevWAM,  /ylog, xrange=[-20,40], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM/WAM Zonal Wind Zonal Mean March Mean UT00 45S',$
         xtitle = 'Zonal Wind ', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
         charsize = 1.3 

  oplot, REFORM(uNZMWAM[23,*,0]), zplevWAM, linestyle=2
 
  Title = ' WACCM-X Zonal Wind Zonal Mean March Mean UT00 '
  CONTOUR, REFORM(uNWZM[*,*,0]),lat, zplev, 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(uNWZM[*,*,0]),lat, zplev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[pBot,pTop],xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

  Title = ' WAM Zonal Wind Zonal Mean March Mean UT00 '
  CONTOUR, REFORM(uNZMWAM[*,*,0]),latWAM, zplevWAM, 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(uNZMWAM[*,*,0]),latWAM, zplevWAM, 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'


END