#
# prepare the SEISM probabilities
#
awk '{print $2,$1,10.^($3)}' < ../02_models/SEISM/mag5_0_data_5yr.txt > secinv.gmt
xyz2grd secinv.gmt -Gjunk.grd -I0.2 `grdinfo grad.grd -I0.1/0.1` -r -diNaN
grdsample junk.grd -nl+t0.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -r -Gmask.grd
surface secinv.gmt -I.01 `grdinfo grad.grd -I0.1/0.1` -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 -GSEISM.grd
#