; $Id: rslice3z.pro,v 1.2 2006/05/05 05:27:34 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_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,.1 !p.charsize=2 !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));,yr=[-1,1]*1e+1 plot,zzz,xz_slice1(l,*),ytit='u',yst=3 plot,zzz,xz_slice2(l,*),ytit='lnrho',yst=3;,yr=[-1,1]*1e+1 plot,zzz,xz_slice3(l,*),ytit='s',yst=3;,yr=[-1,1]*1e+1 wait,w end ; !p.multi=0 close,1 close,2 close,3 END