#
# prepare the kreemer strain map
#
#   convert to nanostrain per year
#   3.1557e7 sec/yr *10-6 nonostrain per 
#
awk '{print $1,$2,sqrt($3*$3+$4*$4+2*$5*$5)*31.557}' ../kreemer/strain_per_sec.gmt > secinv.gmt
xyz2grd secinv.gmt -Gsecinv.grd -I0.02 -R-143/-97/21.0/57.2 -NNaN
grdsample secinv.grd `grdinfo grad.grd -I0.2/0.2` -I0.00833333/0.00833333 -F -Gkreemer_strain.grd