#
# prepare the parsons strain map
#
awk '{print $1,$2,$5,$6,$7,$2-31.5+$1+124.5 }' < ../02_models/parsons/strain.dat > tmp.dat
awk '{ if($6 > 4.99 && $6 < 15.5) print $1,$2,sqrt($3*$3+$4*$4+2*$5*$5) }' < tmp.dat > secinv.gmt
xyz2grd secinv.gmt -Gjunk.grd -I0.4 `minmax secinv.gmt -I.2` -r -diNaN
grdsample junk.grd  -nl+t.1 -I.02 `minmax secinv.gmt -I.2` -r -Gmask.grd
surface secinv.gmt -I.02 `minmax secinv.gmt -I.2` -T.5 -Ll0. -r  -Gsurf.grd -V
grdmath surf.grd mask.grd OR = secinv.grd
grdsample secinv.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -r -Gparsons_strain.grd
#
#  clean up
#
rm eyy exy all.dat secinv.gmt secinv.grd  surf.grd mask.grd junk.grd tmp.dat