#
# first vompute the gradient of the velocity
#
awk '{print $1, $2, $3}' < Vels.dat > Ve.dat
awk '{print $1, $2, $4}' < Vels.dat > Vn.dat
surface Ve.dat `minmax Vels.dat -I4` -I10 -Gve.grd -V
surface Vn.dat `minmax Vels.dat -I4` -I10 -Gvn.grd -V

grdgradient ve.grd -Gvxx.grd -A270 -V
grdgradient ve.grd -Gvxy.grd -A180 -V
grdgradient vn.grd -Gvyx.grd -A270 -V
grdgradient vn.grd -Gvyy.grd -A180 -V
#
grdmath vxx.grd 1000000. DIV = exx.grd
grdmath vyy.grd 1000000. DIV = eyy.grd
grdmath vxy.grd vyx.grd ADD 2000000. DIV = exy.grd
#