Making Map-View Horizon Data

#!/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/aa2.area.xyz | \

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

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

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

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

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

bhpwritecube init=yes key2=cdp,2225,2,214 key1=ep,825,2,257 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