Exercise 6 Answer

#!/bin/sh

project_path=/hou/data/D_152_001/SU_Class/class_project

cd ${project_path}

# make pathlist files and delete outputs

cp int_datasets/template.dat int_datasets/hilb_wavelet_cwt.dat

cp int_datasets/template.dat int_datasets/hilb_wavelet_cwt_slice.dat

cp int_datasets/template.dat int_datasets/hilb_wavelet_cwt_slice_transpose_cdp.dat

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

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

bhpio delete=yes filename=hilb_wavelet_cwt_slice_transpose_cdp \

pathlist=int_datasets/hilb_wavelet_cwt_slice_transpose_cdp.dat

# make a wavelet for cwt

bhpwavelet > seismic/wavelet.su

# take transform and smooth output

bhpread filename=hilb_seismic pathlist=int_datasets/hilb_seismic.dat keys=ep,cdp \

 keylist=7025-9263[80]:18858-21494[80] | \

bhpcwt first=0 inc=0.05 last=1.0 wavelet=seismic/wavelet.su | \

bhpsmooth npass=3 length=20 key=cdp nxmax=21 | \

bhpwritecube init=yes key1=ep,7025,80,29 key2=cdp,18858,80,34 key3=cdpt,1,1,21 \

 filename=hilb_wavelet_cwt pathlist=int_datasets/hilb_wavelet_cwt.dat

# load horizons for proportional slice, flatten horizons at 6400,7300,8200,8500

bhpread filename=hilb_wavelet_cwt pathlist=int_datasets/hilb_wavelet_cwt.dat | \

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,80,34 key2=cdpt,1,1,51 key3=ep,7025,80,29 \

 filename=hilb_wavelet_cwt_slice pathlist=int_datasets/hilb_wavelet_cwt_slice.dat

# make transpose of slice data

bhptranspose filename=hilb_wavelet_cwt_slice pathlist=int_datasets/hilb_wavelet_cwt_slice.dat