#
#
# first construct the velocity grids
#
awk ' {print $1, $2, $3}' < strain1.txt > llve.dat
awk ' {print $1, $2, $4}' < strain1.txt > llvn.dat
surface llve.dat -T.5 -R-123.0/-114./31.5/37.5 -I.01 -V -r -Gve_tot.grd
surface llvn.dat -T.5 -R-123.0/-114./31.5/37.5 -I.01 -V -r -Gvn_tot.grd
#
#   script to extract the strain components
#   (Note for Shen one cannot compute strain from velocity.)
#
awk '{print $1, $2, $6}' < strain1.dat > exx.dat
awk '{print $1, $2, $7}' < strain1.dat > eyy.dat
awk '{print $1, $2, $8}' < strain1.dat > exy.dat