process.combined.OLD.sh

#! /bin/sh

### GRID="top_reservoir.twt.grd"
### GRID="JJO_Top_C01_Peak_TWT_msec.dat.trimmed.smooth.grd"
### GRID="JJO_80MaSB_Peak_FINAL_TWT_msec.dat.trimmed.smooth.grd"

### --------------------------- Start Here ------------------------------------- ###

SURFACE="Top_Reservoir_USER_TWT_msec_19Nov2013_picks.dat"
### SURFACE="Top_C_Sand_TWT.reform.dat"
### SURFACE="JJO_80MaSB_Peak_FINAL_TWT_msec.dat"

FNAME=`echo ${SURFACE} | awk -F".dat" '{print $1}'`
GRID="${FNAME}.grd"

RANGE_VELOCITIES="-R511797/527792/152443/169653"
MINMAX=`minmax -C ${SURFACE}`
LOW_TWT=`echo ${MINMAX} | awk '{print $5}'`
HIGH_TWT=`echo ${MINMAX} | awk '{print $6}'`

surface ${SURFACE} -G${GRID} -S2000 -T0.5 -I721+ ${RANGE_VELOCITIES} -Ll${LOW_TWT} -Lu${HIGH_TWT} -V

### FNAME=`echo ${GRID} | awk -F".grd" '{print $1}'`

echo "Intersecting checkshots with horizon = ${FNAME}"

rm -f ${FNAME}.checkshots.intersect.dat
while read -r LINE ; do
   intersect_checkshot_with_grid verbose=1 pfile=./checkshots/${LINE} coeff_x=${GRID} >> ${FNAME}.checkshots.intersect.dat
done < ./checkshots/checkshots.reform.lis
echo
sort ${FNAME}.checkshots.intersect.dat | grep -v -e "nan" > bub
mv bub ${FNAME}.checkshots.intersect.dat
cat ${FNAME}.checkshots.intersect.dat

echo "Extracting seismic average velocities from SEG-Y volume at the TWT horizon"

suextract < Cb_Ph3_TWT_Vavg_062013_32b.su coeff_x1=${GRID} verbose=1 > stuff


make.gmt.grid.sh -s stuff
sleep 10

echo "Compare extracted seismic with reference well average velocities"

GRID1="stuff.trimmed.sample.smooth.grd"

gmtset D_FORMAT %0.4lf
rm -f output.dat
while read -r LINE ; do
   X=`echo ${LINE} | awk '{print $1}'`
   Y=`echo ${LINE} | awk '{print $2}'`
   Vavg_Seismic=`echo "${X} ${Y}" | grdtrack -Qb -G${GRID1} | awk '{print $3}' | sed 's/-//g'`
   echo "${LINE} ${Vavg_Seismic}" >> output.dat
done < ${FNAME}.checkshots.intersect.dat

echo "Compute annd apply a polynomial trend relationship to Seismic average velocities"

awk '{print $6, $5}' output.dat > input.dat
FOO=`supoly < input.dat verbose=0 n=3`
A_COEFF=`echo ${FOO} | awk '{print $1}'`
B_COEFF=`echo ${FOO} | awk '{print $2}'`

gmtset D_FORMAT %0.16lf
grdmath ${GRID1} -1 MUL = stuff.grd
grdmath stuff.grd stuff.grd MUL ${B_COEFF} MUL = part1.grd
grdmath stuff.grd ${A_COEFF} MUL part1.grd ADD = trend.grd

echo "Measure the difference between trend-corrected seismic average velocities and reference checkshot values"

GRID2="trend.grd"

gmtset D_FORMAT %0.4lf
rm -f bub
while read -r LINE ; do
   X=`echo ${LINE} | awk '{print $1}'`
   Y=`echo ${LINE} | awk '{print $2}'`
   Vavg_Seismic=`echo "${X} ${Y}" | grdtrack -Qb -G${GRID2} | awk '{print $3}' | sed 's/-//g'`
   echo "${LINE} ${Vavg_Seismic}" >> bub
done < output.dat

rm -f output1.dat
while read -r LINE ; do
   DIFF=`echo ${LINE} | awk '{printf "%16.8f", $5-$7}'`
   echo "${LINE} ${DIFF}" >> output1.dat
done < bub

echo "Grid the residual velocity corrections and apply to the trend-corrected values"

awk '{print $1, $2, $8}' output1.dat > residual.vavg.dat
DXDY=`grdinfo -I ${GRID2}`
RANGE=`grdinfo -I721+ ${GRID2}`
MINMAX=`minmax -C residual.vavg.dat`
LOW=`echo ${MINMAX} | awk '{print $5}'`
HIGH=`echo ${MINMAX} | awk '{print $6}'`

echo "DXDY = ${DXDY}, RANGE = ${RANGE}, LOW = ${LOW}, HIGH = ${HIGH}"

surface residual.vavg.dat -Gstuff.grd -S16000 -T0 -I721+ ${RANGE} -Ll${LOW} -Lu${HIGH} -V
grdsample stuff.grd -Gresidual.vavg.grd ${DXDY} ${RANGE} -V -Q
grdmath residual.vavg.grd trend.grd ADD = ${FNAME}.vavg.grd

gmtset D_FORMAT %0.2lf
mbm_grdplot -I${FNAME}.vavg.grd -X -V -G2 -A3/30 -B2000
sleep 10

GRID3="${FNAME}.vavg.grd"

echo "Measure the final average-velocity differences after application of residual corrections"

gmtset D_FORMAT %0.4lf
rm -f bub
while read -r LINE ; do
   X=`echo ${LINE} | awk '{print $1}'`
   Y=`echo ${LINE} | awk '{print $2}'`
   Vavg_Seismic=`echo "${X} ${Y}" | grdtrack -Qb -G${GRID3} | awk '{print $3}' | sed 's/-//g'`
   echo "${LINE} ${Vavg_Seismic}" >> bub
done < output1.dat

rm -f output2.dat
while read -r LINE ; do
   DIFF=`echo ${LINE} | awk '{printf "%16.8f", $5-$9}'`
   echo "${LINE} ${DIFF}" >> output2.dat
done < bub

cp output2.dat ${FNAME}.vavg.output.dat
grdmath ${FNAME}.vavg.grd -1 MUL = ${FNAME}.vavg.neg.grd
gmtset D_FORMAT %0.2lf
mbm_grdplot -I${FNAME}.vavg.neg.grd -X -V -G2 -A3/30 -B2000
sleep 10
ps2pdf ${FNAME}.vavg.neg.grd.ps
acroread ${FNAME}.vavg.neg.grd.pdf &
cat output2.dat