#
# prepare the std strain map
#
grd2xyz ../08_average/exx_std.grd -s | awk '{print $1,$2,$3}' > exx
grd2xyz ../08_average/eyy_std.grd -s | awk '{print $3}' > eyy
grd2xyz ../08_average/exy_std.grd -s | awk '{print $3}' > exy
paste exx eyy exy > all.dat
awk '{print $1,$2,sqrt($3*$3+$4*$4+2*$5*$5) }' < all.dat > secinv.gmt
xyz2grd secinv.gmt -Gjunk.grd -I60s `grdinfo grad.grd -I120s` -r -diNaN
grdsample junk.grd -I60s -nl+t0.1 `grdinfo grad.grd -I120s` -r -Gmask.grd
surface secinv.gmt -I60s `grdinfo grad.grd -I120c` -T.5 -Ll0. -r -Gsurf.grd -V
grdmath surf.grd mask.grd OR = secinv.grd
grdsample secinv.grd -nb+t.5 `grdinfo grad.grd -I90s` -I30s -r -Gstd_strain.grd
#
#  clean up
#
#rm eyy exy all.dat secinv.gmt secinv.grd  surf.grd mask.grd junk.grd