#!/bin/bash
#
# make a plot of GPS velocity vectors 
#
rm -f FIG01.ps
#
D=blk.lluv
R=-R123.0W/114W/31.5N/37.5N
U=GPS_u.grd
V=GPS_v.grd
gmtset  FORMAT_GEO_MAP = dddF
gmt select $D $R -fg | awk '{ print($0," 0 ") }' > data.lluvenct
#
#   first make a mask
#
gmt grdlandmask -Gmask.grd -R$U -Df
gmt grdmath $U mask.grd MUL = GPS_u.nc
gmt grdmath $V mask.grd MUL = GPS_v.nc
#
#   make the map
#
gmt pscoast $R -JM7i -P -Glightgray -Ba1f30m -BWSne -K -Df -X0.7i -W1 > FIG01.ps
gmt psxy fault_data.gmt -J -R -W.5 -O -V -K >> FIG01.ps
gmt psvelo data.lluvenct -J $R -Se.008i/0.95/8 -A9p -W.3,red -O -K >> FIG01.ps
gmt grdvector GPS_u.nc GPS_v.nc -Ix12/12 -J $R -O -Q0.03i+e+n0.03i -Gblue -W.3,blue -S1.i --MAP_VECTOR_SHAPE=0.2 >> FIG01.ps