; $Id: rslice3b.pro,v 1.1 2008/12/03 10:30:29 brandenb Exp $
;
;  show two different slices
;
print,'for 1-d fields in the x-direction'
;
file_slice1='data/proc0/slice_chemspec3.xz'
file_slice2='data/proc0/slice_chemspec4.xz'
;
file_slice1='data/proc0/slice_chemspec.xz'
file_slice2='data/proc0/slice_chemspec2.xz'
file_slice3='data/proc0/slice_chemspec9.xz'
;
xz_slice1=one*fltarr(nx,nz)
xz_slice2=one*fltarr(nx,nz)
xz_slice3=one*fltarr(nx,nz)
;
close,1 & openr,1,file_slice1,/f77
close,2 & openr,2,file_slice2,/f77
close,3 & openr,3,file_slice3,/f77
print,'files opened ok'
;
;  set default waiting time w
;  real all snapshots and plot
;
n=0
default,w,.02
!p.charsize=2
!p.multi=[0,1,2]
!x.margin=[15,2]
xxx=x[l1:l2]
;
while not eof(1) do begin
  readu,1,xz_slice1,t
  readu,2,xz_slice2,t
  readu,3,xz_slice3,t
  print,t
  ee=(xz_slice1-xz_slice2)/(xz_slice1+xz_slice2)
  contour,ee,/fil,lev=grange(-1,1,30)*.1
  ;contour,xz_slice2,/fil,lev=grange(0,1,30)
  contour,xz_slice3,/fil,lev=grange(0,1,30)
  wait,w
end
;
!p.multi=0
close,1
close,2
END