;$Id: ppuymzm.pro,v 1.8 2006/12/21 06:49:06 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
;  mv idl.ps ~/tex/procs/prague06b/fig/ppuymzm_256_512a_Om02.ps
;
nlev=20
bar='!20!s!u$!n!r!6'
gray=210 & thick=20 & dark=160
!p.charsize=1.8
;!x.margin=[10.6,1.2]
;
loadct,0
!x.title='!8z!6/!8d!6'
!y.title=bar+'u!8!dy!n!6/!9S!6(!8g!6!8d!6)'
!p.title='!6'
restore,'uymz.sav'
!y.range=[1.05*min(fml),1.1*max(fmu)]
plot,zzz,fm
;
;  plot range of unstable region
;
z1=min(zzz)
z2=max(zzz)
polyfill,[z1,0.,0.,z1],[[1,1]*!y.range(0),[1,1]*!y.range(1)],col=gray
polyfill,[z2,1.,1.,z2],[[1,1]*!y.range(0),[1,1]*!y.range(1)],col=gray
;
;  plot extra dash-dotted line and plot thick line at the bottom
;
oplot,[z1,0.,0.,z1],[[1,1]*!y.range(0),[1,1]*!y.range(1)],li=3
oplot,[z2,1.,1.,z2],[[1,1]*!y.range(0),[1,1]*!y.range(1)],li=3
oplot,[z1,0.],[1,1]*!y.range(0),thick=thick
oplot,[1.,z2],[1,1]*!y.range(0),thick=thick
;
;  overplot data and error strip
;
polyfill,[zzz,reverse(zzz)],[fml,reverse(fmu)],col=dark
plot,zzz,fm,/noerase
;
END