#
#
# first construct the velocity grids
#
awk ' {print $1-360., $2, $14}' < qk07_grid_small.strain.out.txt > llve.dat
awk ' {print $1-360., $2, $15}' < qk07_grid_small.strain.out.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-360., $2, $3}' < qk07_grid_small.strain.out.txt > exx.dat
awk '{print $1-360., $2, $5}' < qk07_grid_small.strain.out.txt > eyy.dat
awk '{print $1-360., $2, $4}' < qk07_grid_small.strain.out.txt > exy.dat
surface exx.dat -T.5 -R-123.0/-114./31.5/37.5 -I.01 -V -r -Gexx.grd
surface exy.dat -T.5 -R-123.0/-114./31.5/37.5 -I.01 -V -r -Gexy.grd
surface eyy.dat -T.5 -R-123.0/-114./31.5/37.5 -I.01 -V -r -Geyy.grd