#!/bin/csh
#
# set some rules
#
set file = ps.hay_seis
set cpath = /usr/contrib/data/maps
set data = /data/05/lind/graphics/gmt/data
set black = 0/0/0
set white = 255/255/255
set grey = 130
set red = 255/0/0
set green = 0/255/0
set blue = 0/0/255
set yellow = 255/255/0
set blgn = 150/150/250
#
# create the basemap
#
psbasemap -Jl-122.2/37.4/38.2/38/1:650000 -R-122.8/-121.6/37.30/38.30 -Ba0.25WSne -L-122.60/37.4/38./20. -Y2.0 -P -K >! $file
#
# add the California state outline
#
psxy $cpath/caoutline.asc -Jl -R -Ms -W8/$black -: -O -K >> $file
#
# add the faults
#
psxy $cpath/cafaults.asc -Jl -R -Ms -W4/$red -: -O -K >> $file
#
# identify some towns
#
#cat $data/places.bay | awk '{print $1, $2}' >! temp
#psxy temp -Jl -R -Sh0.05 -L -W2/$black -O -K >> $file
#grep '0 6 7' $data/places.bay | awk '{print $1-0.01, $2, $3, $4, $5, $6, $7, $8, $9}' >! temp
#pstext temp -Jl -R -O -K >> $file
#grep '0 6 5' $data/places.bay | awk '{print $1+0.02, $2, $3, $4, $5, $6, $7, $8, $9}' >! temp
#pstext temp -Jl -R -O -K >> $file
#
# NCSN station locations
#
cat $data/menlo.stat | awk '{print $2, $3}' >! temp
psxy temp -Jl -R -L -Sc0.07 -G$blgn -W1/$blgn -: -O -K  >> $file
#
# UCB BDSN Stations
#
cat $data/ncal.stat | awk '{print $2, $3}' >! temp
psxy temp -Jl -R -L -Ss0.10 -G$blue -W1/$blue -: -O -K  >> $file
#
# future BDSN station locations
#
cat $data/ncal.future | awk '{print $2, $3}' >! temp
psxy temp -Jl -R -L -Ss0.10 -G$white -W1/$blue -: -O -K  >> $file
#
# UCB HF Stations
#
psxy $data/hfn.ucb.stat -Jl -R -Sd0.10 -G$blue -W1/$blue -L -: -O -K>>$file
#
# USGS HF Stations
#
psxy $data/hfn.usgs.stat -Jl -R -Sd0.10 -G$blue -W1/$blue -L -: -O -K>> $file
#
# future UCB HF sites
#
psxy $data/hfn.ucb_planned.stat -Jl -R -Sd0.10 -G$white -W2/$blue -L -: -O -K>> $file
#
# Bridge sites
#
psxy $data/hfn.bridge.stat -Jl -R -St0.10 -G$blue -W1/$blue -L -: -O -K>> $file
#
# Planned bridge sites
#
psxy $data/hfn.bridge_planned.stat -Jl -R -St0.10 -G$white -W2/$blue -L -: -O -K>> $file
#
# add labels
#
# UCB BDSN Stations
#
psxy -Jl -R -Ss0.10 -G$blue -W1/$blue -L -: -O -K>>$file<<EOF
37.70 -122.76
EOF
#
# UCB-USGS HF Stations
#
psxy -Jl -R -Sd0.10 -G$blue -W1/$blue -L -: -O -K>>$file<<EOF
37.65 -122.76
EOF
#
# Bridge sites
#
psxy -Jl -R -St0.10 -G$blue -W1/$blue -L -: -O -K>> $file<<EOF
37.60 -122.76
EOF
#
# NCSN sites
#
psxy -Jl -R -Sc0.07 -G$blgn -W1/$blgn -L -: -O -K>> $file<<EOF
37.55 -122.76
EOF

#
# Planned bridge sites
#
psxy -Jl -R -St0.10 -G$white -W2/$blue -L -: -O -K>> $file<<EOF
37.50 -122.72
EOF
#
# Planned UCB HF sites
#
psxy -Jl -R -Sd0.1 -G$white -W2/$blue -L -: -O -K>> $file<<EOF
37.50 -122.76
EOF
#
pstext  -Jl -R -O -K >> $file <<EOF
-122.74 37.70 10 0.0 25 4 BDSN: UCB
-122.74 37.65 10 0.0 25 4 HFN: UCB & USGS
-122.74 37.60 10 0.0 25 4 Bridge Program
-122.74 37.55 10 0.0 25 4 NCSN: USGS
-122.70 37.50 10 0.0 25 4 Planned sites
EOF
#
# add a title
#
pstext  -Jl -R -N -O >> $file <<EOF
-122.2 37.22 18 0.0 25 6 Bay Area Digital Seismic Monitoring
EOF
#
#
# display
#
ghostview $file