Proportional Slice

 

#!/bin/sh

 

cd ${project_directory}

 

cp int_datasets/default.dat int_datasets/hilb_seismic_slice.dat

bhpio delete=yes filename=hilb_seismic_slice pathlist=int_datasets/hilb_seismic_slice.dat

 

bhpread filename=hilb_seismic pathlist=int_datasets/hilb_seismic.dat keys=ep,cdp keylist=7025-9263:18858-21494 | \

bhploadhdr key=ep,cdp,f1 infile=horizons/M2aSB.xyz gridtype=R extrap=no | \

bhploadhdr key=ep,cdp,f2 infile=horizons/M7SB.xyz gridtype=R extrap=no | \

bhploadhdr key=ep,cdp,d1 infile=horizons/top_eocene_hie.xyz gridtype=R extrap=no | \

bhploadhdr key=ep,cdp,d2 infile=horizons/marker_6.xyz gridtype=R extrap=no | \

bhpslice hlist=f1,f2,d1,d2 times=6400,7300,8200,8500 tmin=5000 tmax=9500 dt=4 | \

sushw key=f1 a=5.0 | \

sushw key=f2 a=0 | \

sushw key=d1 a=0 | \

sushw key=d2 a=0 | \

bhpwritecube init=yes key1=cdp,18858,2,1319 key2=ep,7025,2,1120 \

filename=hilb_seismic_slice pathlist=int_datasets/hilb_seismic_slice.dat