#
# prepare the loveless_meade strain map
#
grep -v NaN < ../loveless_meade/exx.dat | awk '{print $1-360.,$2,$3*-1.e6}' > exx.dat
xyz2grd exx.dat -Gjunk.grd -I0.1 `grdinfo grad.grd -I0.1/0.1` -r -diNaN
grdsample junk.grd -nl+t.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -r -Gmask.grd
surface exx.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5  -r -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -r -Gloveless_meade_exx.grd
#
grep -v NaN < ../loveless_meade/eyy.dat | awk '{print $1-360.,$2,$3*-1.e6}' > eyy.dat
xyz2grd eyy.dat -Gjunk.grd -I0.1 `grdinfo grad.grd -I0.1/0.1` -r -diNaN
grdsample junk.grd -nl+t.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -r -Gmask.grd
surface eyy.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5  -r -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -r -Gloveless_meade_eyy.grd
#
grep -v NaN < ../loveless_meade/exy.dat | awk '{print $1-360.,$2,-$3*1.e6}' > exy.dat
xyz2grd exy.dat -Gjunk.grd -I0.1 `grdinfo grad.grd -I0.1/0.1` -r -diNaN
grdsample junk.grd -nl+t.1 -I.01 `grdinfo grad.grd -I0.1/0.1` -r -Gmask.grd
surface exy.dat -I.01 `grdinfo grad.grd -I0.1/0.1` -T.5  -r -Gsurf.grd -V
grdmath surf.grd mask.grd OR = tmp.grd
grdsample tmp.grd `grdinfo grad.grd -I0.1/0.1` -I0.00833333/0.00833333 -r -Gloveless_meade_exy.grd
#
#  clean up
#
rm surf.grd mask.grd junk.grd exx.dat eyy.dat exy.dat tmp.grd