;$Id: puy_ytslice.pro,v 1.14 2006/12/21 06:49:06 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=14,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
end
;
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_512a_Om02.ps
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_512a_Om02_nocross.ps
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_128a_Om02_theta90_fresh_up3.ps
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_128a_Om01_theta90_fresh_up3b.ps
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_128a_Om02_theta90_fresh_up3b.ps
;  mv idl.ps ~/tex/hydro/shearlayer/fig/puy_ytslice_256_128a_Om005_theta90_fresh_up3c.ps
;
;  Note: the savefile data are produced with ryzslice.pro
;  reads uy_ytslice.sav files
;
!x.title='!8t!6 !9S!6(!8g!6/!8d!6)'
!y.title='!8y!6/!8d!6'
!p.charsize=1.9
;
file='uy_ytslice.sav'
file='uu2_ytslice.sav'
restore,file
levmax=.4
lev=grange(-levmax,levmax,16)
x1=.12 & x2=.77 & y1=.15 & y2=.95
!p.position=[x1,y1,x2,y2]
contour,transpose(yt),tt,yyy,lev=lev,/fil
;
;  overplot while cross
;
default,docross,0.
if docross eq 1. then begin
  oplot,[1,1]*254,[-2,2],col=255,thick=5
  oplot,tt,tt-tt+.35,col=255,thick=5
endif
;
d1=.13 & d2=.16
colorbar,/vert,pos=[x2+d1,y1,x2+d2,y2],ran=[-levmax,levmax+1e-5], $
  xtit='!8u!dy!n!6/!9S!6(!8g!6!8d!6)',ytit='!6',form='(f4.1)', $
  charsize=1.6, DIV=8
;
END