#!/bin/csh

/bin/rm tmp.ps

## Plot shear wave splitting parameters on a map.
# The splitting parameters are in a file called meansig.out
# The longitude and latitude of the measurements are in
# the 3rd and second column; the fast orientation is in the 
#	4th and the splitting time is in the 6th columns


# Plot a base map:
pscoast -R28/38/-18/-8 -JM5.0 -Ba2f1 \
-W2 -Df -A100 -Na/1 -P -K -X1.5 -Y2.5 >tmp.ps

# Plot the bars as head-less arrows (-SV)

awk '{print $3, $2, $4, 0.10*($6+0.0)}' meansig.out | \
psxy -R -JM -O -K -SV0/0/0 -L -W18/255/0/0 -V >> tmp.ps

# Plot the other half of the bars as head-less arrows (-SV)

awk '{print $3, $2, $4+180, 0.10*($6+0.0)}' meansig.out | \
psxy -R -JM -O -K -SV0/0/0 -L -W18/255/0/0 -V >> tmp.ps

# Plot a dot at the station

awk '{print $3, $2}' meansig.out | \
psxy -R -JM -O -Sc0.035 -G255/0/0 >> tmp.ps