Exercise 4 Answer

#!/bin/sh

project_directory=/hou/data/D_152_001/SU_Class/class_project

# make pathlist file

cp ${project_directory}/int_datasets/template.dat ${project_directory}/int_datasets/tiger_inv.dat

# remove in case already there

bhpio delete=yes filename=tiger_inv pathlist=${project_directory}/int_datasets/tiger_inv.dat

# import inversion.sgy to tiger_inv

segyread tape=${project_directory}/seismic/inversion.sgy | segyclean | suchw key1=ep key2=fldr | \

bhpwritecube init=yes filename=tiger_inv pathlist=${project_directory}/int_datasets/tiger_inv.dat \

key1=ep,839,2,11 key2=cdp,2225,2,214

#verify

bhpio filename=tiger_inv pathlist=${project_directory}/int_datasets/tiger_inv.dat

# load aa2_trough and aa2_peak picks and run bhpintegrate

# use bhpstorehdr to save each amplitude extraction event

# don't need to save data, so send to bit bucket

bhpread filename=tiger_inv pathlist=${project_directory}/int_datasets/tiger_inv.dat | \

bhploadhdr key=ep,cdp,selev infile=${project_directory}/horizons/aa2_trough.xyz gridtype=R extrap=no bias=-16 verbose=1 | \

bhploadhdr key=ep,cdp,gelev infile=${project_directory}/horizons/aa2_peak.xyz gridtype=R extrap=no bias=16 | \

bhpintegrate mode=edge type=trough search=5 limit=5 P1=selev P2=gelev | \

bhpstorehdr key=ep,cdp,d2 ofile=${project_directory}/horizons/area.xyz | \

bhpstorehdr key=ep,cdp,f1 ofile=${project_directory}/horizons/width.xyz | \

bhpstorehdr key=ep,cdp,f2 ofile=${project_directory}/horizons/amp.xyz | \

bhpstorehdr key=ep,cdp,selev ofile=${project_directory}/horizons/p1.xyz | \

bhpstorehdr key=ep,cdp,gelev ofile=${project_directory}/horizons/p2.xyz | \

bhpstorehdr key=ep,cdp,corr ofile=${project_directory}/horizons/p3.xyz > /dev/null