#
# prepare the mccaffrey strain map
#
awk '{print $3}' <  ../mccaffrey/eyy.dat > eyy 
awk '{print $3}' <  ../mccaffrey/exy.dat > exy 
paste ../mccaffrey/exx.dat eyy exy > all.dat
awk '{print $1,$2,sqrt($3*$3+$4*$4+2*$5*$5)*1.e9 }' < all.dat > secinv.gmt
xyz2grd secinv.gmt -Gjunk.grd -I0.02 `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 secinv.gmt -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5 -Li0. -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 -F -Gmccaffrey_strain.grd
#
#   low low-pass filter with a 20 km gaussian
#
grdfilter mccaffrey_strain.grd -Fg20 -D2 -Gmccaffrey_filtered_strain.grd -V
#
#  clean up
#
#rm eyy exy all.dat secinv.gmt secinv.grd  surf.grd mask.grd junk.grd