#!/bin/csh
#
# set some rules
#
set file = ps.hay_seis
set cpath = /usr/contrib/data/maps
set data = /data/a1/lind/graphics/gmt/data
set black = 0/0/0
set white = 255/255/255
set grey = 130
set shade = 230
set red = 255/0/0
set redta = 255/0/0ta
set redto = 255/0/0to
set pink = 255/110/255
set pinkta = 255/110/255ta
set pinkto = 255/110/255to
set orange = 255/110/80
set orangeta = 255/110/80ta
set orangeto = 255/110/80to
set green = 0/255/0
set blue = 0/0/255
set tour = 0/255/255
set yellow = 255/255/0
set pblue = 220/240/255
set lblue = 180/200/255
set pyellow = 255/255/160
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.25f7.5mWSne -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
#
# draw the Bridges
#
psxy $data/bridge.asc -Jl -R -MB -O -K -W15/175 -: >> $file
#
# add the faults
#
# - old faults
#psxy $cpath/cafaults.asc -Jl -R -Ms -W4/$red -: -O -K >> $file
# - new faults
cat /usr/contrib/data/maps/jennings.92.hist.asc | grep solid | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W5/$red -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.hist.asc | grep dashed | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/$redta -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.hist.asc | grep dotted | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/$redto -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.holo.asc | grep solid | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W5/$orange -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.holo.asc | grep dashed | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/$orangeta -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.holo.asc | grep dotted | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/$orangeto -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.quat.asc | grep solid | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W1/0/0/0 -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.quat.asc | grep dashed | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/0/0/0ta -: -O -K >> $file
cat /usr/contrib/data/maps/jennings.92.quat.asc | grep dotted | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -H2 -M$ -W3/0/0/0to -: -O -K >> $file
#
# identify some towns
#
cat $data/places.bay | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -St0.10 -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/ncsn.stat | awk '{print $2, $3}' >! temp
psxy temp -Jl -R -L -Sc0.07 -G$blgn -W1/$black -: -O -K  >> $file
#
# UCB BDSN Stations
#
cat $data/bdsn.stat | grep BDSN | grep OPER | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -L -Ss0.10 -G$blue -W1/$black -O -K  >> $file
#
# future BDSN station locations
#
cat $data/bdsn.stat | grep BDSN | grep PLAN | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -L -Ss0.10 -G$white -W1/$blue -O -K  >> $file
#
# UCB & USGS HF Stations
#
cat $data/bdsn.stat | grep HFN | grep OPER | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -Sd0.10 -G$blue -W1/$black -L -O -K>>$file
#
# future UCB HF sites
#
cat $data/bdsn.stat | grep HFN | grep PLAN | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -Sd0.10 -G$white -W1/$blue -L -O -K>> $file
#
# Bridge sites
#
cat $data/bdsn.stat | grep BPN | grep OPER | awk '{print $1, $2}' >! temp
psxy temp -Jl -R -St0.10 -G$blue -W1/$black -L -O -K>> $file
#
# Planned bridge sites
#
cat $data/bdsn.stat | grep BPN | grep PLAN | awk '{print $1, $2}' >! temp
psxy temp -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/$black -L -: -O -K>>$file<<EOF
37.70 -122.76
EOF
#
# UCB-USGS HF Stations
#
psxy -Jl -R -Sd0.10 -G$blue -W1/$black -L -: -O -K>>$file<<EOF
37.65 -122.76
EOF
#
# Bridge sites
#
psxy -Jl -R -St0.10 -G$blue -W1/$black -L -: -O -K>> $file<<EOF
37.60 -122.76
EOF
#
# NCSN sites
#
psxy -Jl -R -Sc0.07 -G$blgn -W1/$black -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 Seismic Monitoring
EOF
#
#
# display
#
ghostview $file