; $Id: rslice3z.pro,v 1.3 2006/05/29 22:02:40 brandenb Exp $
;
;  show two different slices
;
;spawn,'sh ../bin/READ.sh'
print,'for 1-d fields in the x-direction'
;
file_slice1='data/proc0/slice_uz.xz'
file_slice2='data/proc0/slice_lnrho.xz'
file_slice3='data/proc0/slice_lnTT.xz'
;
t=zero
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
;
l=0
default,w,.1
!p.charsize=4
!p.multi=[0,1,3]
!x.margin=[15,2]
zzz=z(n1:n2)
;
while not eof(1) do begin
  readu,1,xz_slice1,t
  readu,2,xz_slice2,t
  readu,3,xz_slice3,t
  print,t
  ;plot_io,exp(xz_slice1(*,n))
  plot,zzz,xz_slice1(l,*),ytit='u',yst=3;,yr=[-1,1]
  plot,zzz,xz_slice2(l,*),ytit='lnrho',yst=3;,yr=[-1,1]*1e+1
  plot,zzz,exp(xz_slice3(l,*)),ytit='T',yst=3;,yr=[0,1]*5e5
  wait,w
end
;
!p.multi=0
close,1
close,2
close,3
END