#
# prepare the parsons strain map
#
awk '{print $1,$2,$5,$6,$7,$2-31.5+$1+124.5 }' < ../parsons/strain.dat > tmp.dat
awk '{ if($6 > 4.99 && $6 < 15.5) print $1, $2, $3, $4, $5 }' < tmp.dat > parsons.dat
awk '{print $1,$2,$3}' < parsons.dat > exx.dat
xyz2grd exx.dat -Gjunk.grd -I0.4 `grdinfo grad.grd -I0.1/0.1` -NNaN
grdsample junk.grd -Ql.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -Gmask.grd
surface exx.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5 -Li0. -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -F -Gparsons_exx.grd
#
awk '{print $1,$2,$4}' parsons.dat > eyy.dat
xyz2grd eyy.dat -Gjunk.grd -I0.4 `grdinfo grad.grd -I0.1/0.1` -NNaN
grdsample junk.grd -Ql.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -Gmask.grd
surface eyy.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5 -Li0. -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -F -Gparsons_eyy.grd
#
awk '{print $1,$2,$5}' parsons.dat > exy.dat
xyz2grd exy.dat -Gjunk.grd -I0.4 `grdinfo grad.grd -I0.1/0.1` -NNaN
grdsample junk.grd -Ql.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -Gmask.grd
surface exy.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5 -Li0. -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -F -Gparsons_exy.grd
#
#  clean up
#
rm surf.grd mask.grd junk.grd exx.dat eyy.dat exy.dat tmp.grd