#
#
# first take the gradidient of the velocity grids
#
gmt grdgradient ve_tot.grd -Gvxx.grd -A270 -fg -V
gmt grdgradient ve_tot.grd -Gvxy.grd -A180 -fg -V
gmt grdgradient vn_tot.grd -Gvyx.grd -A270 -fg -V
gmt grdgradient vn_tot.grd -Gvyy.grd -A180 -fg -V
#
gmt grdmath vxx.grd 1000. DIV = exx.grd
gmt grdmath vyy.grd 1000. DIV = eyy.grd
gmt grdmath vxy.grd vyx.grd ADD 2000. DIV = exy.grd
#
gmt grd2xyz exx.grd | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 16.0) print $1,$2,$3}' > exx.dat
gmt grd2xyz eyy.grd | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 16.0) print $1,$2,$3}' > eyy.dat
gmt grd2xyz exy.grd | awk '{print $1,$2,$3,$1-31.5+$2+124.5}' | awk '{if ($4 > 2.99 && $4 < 16.0) print $1,$2,$3}' > exy.dat