; $Id: rslice3.pro,v 1.1 2007/05/12 07:43:40 brandenb Exp $
;
;  show two different slices
;
print,'for 1-d fields in the x-direction'
;
file_slice1='data/proc0/slice_uu3.xz'
file_slice2='data/proc0/slice_lnrho.xz'
file_slice3='data/proc0/slice_ss.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,.2
!p.charsize=2
!p.multi=[0,1,3]
!x.margin=[15,2]
xxx=x(l1:l2)
zzz=z(n1:n2)
;
Ma=1.0 & yru=1.5 & yrlr=[-12,-2] & yrs=[.2,.4]
Ma=1.0 & yru=15. & yrlr=[-12,-2] & yrs=[.2,.4]
Ma=1.0 & yru=.015 & yrlr=[-12,-2] & yrs=[.2,.4]
Ma=1.0 & yru=.15 & yrlr=[-1,1]*7 & yrs=[-1,1]*3
;
i=0
while not eof(1) do begin
  readu,1,xz_slice1,t
  readu,2,xz_slice2,t
  readu,3,xz_slice3,t
  print,t,max(abs(xz_slice1))
  ;plot_io,exp(xz_slice1(*,n));,yr=[-1,1]*1e+1
  plot,zzz,xz_slice1(l,*),ytit='u',yst=3,yr=[-1,1]*yru
  plot,zzz,xz_slice2(l,*),ytit='lnrho',yst=3,yr=yrlr
  plot,zzz,xz_slice3(l,*),ytit='s',yst=3,yr=yrs
  if i eq 0 then begin
    tt=t
    slice1=reform(xz_slice1(l,*))
    slice2=reform(xz_slice2(l,*))
    slice3=reform(xz_slice3(l,*))
  endif else begin
    tt=[tt,t]
    slice1=[slice1,reform(xz_slice1(l,*))]
    slice2=[slice2,reform(xz_slice2(l,*))]
    slice3=[slice3,reform(xz_slice3(l,*))]
  endelse
  wait,w
  i=i+1
end
slice1=reform(slice1,nz,i)
slice2=reform(slice2,nz,i)
slice3=reform(slice3,nz,i)
;
!p.multi=0
close,1
close,2
close,3
END