;$Id: ppuyuzmzm.pro,v 1.2 2006/06/19 15:03:03 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  plots only the mean uy velocity
;  uses data from save file
;  run first puymzm on the compute machine
;
;  mv idl.ps ~/tex/hydro/shearlayer/fig/ppuymzm_256_512a_Om02.ps
;
default,nu,2e-3
@param
nlev=20
bar='!20!s!u$!n!r!6'
!p.charsize=1.8
;!x.margin=[10.6,1.2]
;
loadct,5
!x.title='!8z!6/!8d!6'
!y.title='stress'
!p.title='!6'
restore,'uyuzmz.sav'
stress=fm
;
;  integrate
;
gm=integr(fm,x=zzz)
gml=integr(fml,x=zzz)
gmu=integr(fmu,x=zzz)
!y.range=0
;
restore,'uymz.sav'
flow=fm
;
!p.multi=[0,1,3]
plot,zzz,flow
plot,fm,gm,xr=[-1,1]*.15,yr=[-1,1]*1e-3,ps=2
p=linfit(fm,gm)
print,'p/nu=',p/nu
xx=[-1,1]*.13 & oplot,xx,xx*p(1)+p(0),col=122 & oplot,xx,xx*p(1),col=122
;
plot,deriv(zzz,flow),stress,ps=1,xr=[-1,1]*.4,yr=[-1,1]*1e-3
p=linfit(deriv(zzz,flow),stress)
xx=[-1,1]*.22
oplot,xx,xx*p(1)+p(0),col=122 & oplot,xx,xx*p(1),col=122
print,'p/nu=',p/nu
;
ibot=109 & itop=235
print,'uybot,uytop=',flow(ibot),flow(itop)
;
END