#
#
# first grid the velocities
#
awk 'NR > 1 {print $1, $2, $3}' < ucerf3_zeng_strain.dat | grep -v NaN > llve.dat
awk 'NR > 1 {print $1, $2, $4}' < ucerf3_zeng_strain.dat | grep -v NaN > llvn.dat
surface llve.dat -T.35 -R-125./-114./30.5/42.0 -I.01 -V -r -Gve_tot.grd
surface llvn.dat -T.35 -R-125./-114./30.5/42.0 -I.01 -V -r -Gvn_tot.grd
#
# then take the gradidient of the velocity grids
#
grdgradient ve_tot.grd -Gvxx.grd -A270 -fg -V
grdgradient ve_tot.grd -Gvxy.grd -A180 -fg -V
grdgradient vn_tot.grd -Gvyx.grd -A270 -fg -V
grdgradient vn_tot.grd -Gvyy.grd -A180 -fg -V
#
grdmath vxx.grd 1000. DIV = exx.grd
grdmath vyy.grd 1000. DIV = eyy.grd
grdmath vxy.grd vyx.grd ADD 2000. DIV = exy.grd
#
grd2xyz exx.grd > exx.dat
grd2xyz eyy.grd > eyy.dat
grd2xyz exy.grd > exy.dat