Exercise 5 Answer

#!/bin/sh

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

cd ${project_directory}

# copy pathlist template file

cp int_datasets/template.dat int_datasets/tiger_ampextract.dat

cp int_datasets/template.dat int_datasets/tiger_ampextract_transpose_ep.dat

# delete outputs in case they already exist otherwise, bhpwritecube will quit

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

bhpio delete=yes filename=tiger_ampextract_transpose_ep \ pathlist=int_datasets/tiger_ampextract_transpose_ep.dat

# load ampextract events using picks.su data

bhphorizon < seismic/picks.su phdr=ep shdr=cdp file=horizons/area.xyz | \

bhphorizon phdr=ep shdr=cdp file=horizons/width.xyz | \

bhphorizon phdr=ep shdr=cdp file=horizons/amp.xyz | \

bhphorizon phdr=ep shdr=cdp file=horizons/p1.xyz | \

bhphorizon phdr=ep shdr=cdp file=horizons/p2.xyz | \

bhphorizon phdr=ep shdr=cdp file=horizons/p3.xyz | \

bhpwritecube init=yes key2=cdp,2225,2,214 key1=ep,839,2,11 filename=tiger_ampextract \

 pathlist=int_datasets/tiger_ampextract.dat horizons=area,width,amp,p1,p2,p3

# convert events to map-view   no other parameters are needed

bhptranspose filename=tiger_ampextract pathlist=int_datasets/tiger_ampextract.dat

#verify

bhpio filename=tiger_ampextract pathlist=int_datasets/tiger_ampextract.dat

bhpio filename=tiger_ampextract_transpose_ep pathlist=int_datasets/tiger_ampextract_transpose_ep.dat