#!/bin/csh
#
# make a pretty basemap
#
set file = ps.hay_hist
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
#
# draw 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 state outlines
#
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
#
# plot seismicity
#
cat $data/ellsworth.cat | awk '{print $7, $8, $10*0.035+0.02}' >! temp
psxy temp -Jl -R -Sc -L -W5/$blue -: -O -K >> $file
#
# key
#
cat search.eq.key | awk '{print $7, $8, $10*0.035+0.02}' >! temp
psxy temp -Jl -R -Sc -L -W5/$blue -: -N -O -K >> $file
pstext -Jl -R -N -O -K >> $file <<eof
-121.50 38.10 14 0.0 25 5 M 2
-121.50 38.05 14 0.0 25 5 M 3
-121.50 38.00 14 0.0 25 5 M 4
-121.50 37.95 14 0.0 25 5 M 5
-121.50 37.90 14 0.0 25 5 M 6
-121.50 37.85 14 0.0 25 5 M 7
-121.50 37.80 14 0.0 25 5 M 8
eof
#
# 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
#
# title
#
pstext -Jl -R -N -O >> $file << eof
-122.2 37.22 18 0.0 25 6 Bay Area Seismicity M > 5.5 1769-1997
eof
#
#
ghostview $file &