#
#   first grid the velocity data
#
#cat *.info > vel.dat
minmax vels.dat
awk '{print($1-360,$2,$3)}' < vels.dat > ve.xyz
awk '{print($1-360,$2,$4)}' < vels.dat > vn.xyz
#
#   grid the velocities on a 0.01 degree grid
#
surface ve.xyz -R-125/-113/31/42 -I0.01 -T.35 -V -Gve.grd
surface vn.xyz -R-125/-113/31/42 -I0.01 -T.35 -V -Gvn.grd
#
#  now process just like konter_smith
#
# first take the gradidient of the velocity grids
#
grdgradient ve.grd -Gvxx.grd -A270 -fg -V
grdgradient ve.grd -Gvxy.grd -A180 -fg -V
grdgradient vn.grd -Gvyx.grd -A270 -fg -V
grdgradient vn.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 | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 17.0) print $1,$2,$3}' > exx.dat
grd2xyz eyy.grd | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 17.0) print $1,$2,$3}' > eyy.dat
grd2xyz exy.grd | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 17.0) print $1,$2,$3}' > exy.dat