#!/bin/bash

#####
#    apply a two station techique to get phase velocity for the
#    two stations on the same great circle path
#####
#    the steps are as follow
#    1. Measure the group velocity dispersion from the two traces
#       and then save the one which is better defined and which 
#       has the greatest range of periods
#       do_mft -G *BHT
#       a) pick the dispersion
#       b) click on the "Match"
#          The screen output will show how sacmat96 is called, e.g.,
#          /Users/rbh/PROGRAMS.310t/PROGRAMS.330/bin/sacmat96 -F 2015.171.021121.G.COYC.BHZ -D disp.d -AUTO
#          we will not use the match dispersion, jus the command line
#          We need to use the same command for both stations. We are really interested in the disp.d file

do_mft -G *BHT
#####
#   when Match is selected in do_mft, a command line is printed on the terminal
#   showing the syntax for executing sacmat96. We apply the disp.d to
#   all traces of the same time so that all are processed with the
#   same range of periods
#####
for TRACE in *BHT
do
	sacmat96 -F ${TRACE} -D disp.d -AUTO
done

#####
#    as a result of this the hase matched traces will have an 's' appended 
#    to the file name
#####

#####
#    Before determining the phase velocities, we cut the waveforms from
#    the original length because this 'hack' is required to get do_pom to 
#    work properly, and because sacpom96 requires the number of points
#    to use the same power of 2 when zero filled
#####
gsac << EOF
cut b b 1000
rh *HT*s
w append  .cut
q
EOF

#####
#    now run do_pom
#####
do_pom *HTs.cut

mv POM96.PLT T_POM96.PLT
mv POM96CMP  T_POM96CMP