#
# prepare the loveless_meade strain map
#
# dos2unix exx.dat
# dos2unix eyy.dat
# dos2unix exy.dat
awk '{print $3}' <  ../02_models/loveless_meade/eyy.dat > eyy 
awk '{print $3}' <  ../02_models/loveless_meade/exy.dat > exy 
awk '{print $3}' <  ../02_models/loveless_meade/texx.dat > texx 
awk '{print $3}' <  ../02_models/loveless_meade/teyy.dat > teyy 
awk '{print $3}' <  ../02_models/loveless_meade/texy.dat > texy 
paste ../02_models/loveless_meade/exx.dat eyy exy texx teyy texy > all.dat
awk '{print $1-360.,$2,sqrt(($3+$6)*($3+$6)+($4+$7)*($4+$7)+2*($5+$8)*($5+$8))*1.e6 }' < all.dat > secinv.gmt
xyz2grd secinv.gmt -Gjunk.grd -I0.1 `grdinfo grad.grd -I0.1/0.1` -r -NNaN
grdsample junk.grd  -nl+t0.1 -I.01 `grdinfo grad.grd -I0.2` -r -Gmask.grd
surface secinv.gmt -I.01 `grdinfo grad.grd -I0.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 -Gloveless_meade_strain.grd
#
#  clean up
#
rm eyy exy texx teyy texy all.dat secinv.gmt secinv.grd  surf.grd mask.grd junk.grd