# # prepare the smith_konter strain map # awk '{print $3}' < ../02_models/smith_konter/eyy.dat > eyy awk '{print $3}' < ../02_models/smith_konter/exy.dat > exy paste ../02_models/smith_konter/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 -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 -Gsmith_konter_strain.grd # # clean up # rm eyy exy all.dat secinv.gmt secinv.grd surf.grd mask.grd junk.grd