#!/bin/bash
reg=-R-125/-95/20/45
proj=-JH`echo $reg | gawk -f reg2midlon.awk`/7
ifile=gps.vel
ofile=gps.ps
ann=-Ba5f1WeSn

inc=-I0.1

gawk '{if(NR>3)print($1,$2,$3)}' $ifile  > tmp.vx
triangulate -fg $reg $inc tmp.vx -Gtmp.vx.grd -V > /dev/null

pscoast $reg $proj -P -K -Df -A5000 -S128 > $ofile
gawk '{if(NR>3)print($1,$2,$3,$4,$7,$8,$9,$13)}' $ifile | \
    psvelo -Se0.01/0.95/12 $reg $proj -O  -K >> $ofile

psbasemap $ann $reg $proj -O >> $ofile
modifybb $ofile
echo $0: output in $ifile

