#!/bin/csh
#
# make a pretty basemap
#
set file = ps.hay_eq_ob
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 map corners (ll = lower left, ur = upper right)
# Some Great Valley
 set ba_ll = ( -124.2 38.8 )
 set ba_ur = ( -119.6 37.5 )
# Bay Area only
 set ba_ll = ( -123.66 38.16 )
 set ba_ur = ( -120.4 37.05 )
# Hayward only
 set ha_ll = ( -122.55  38.0 )
 set ha_ur = ( -121.5 37.3 )
#
# draw the basemap
#
psbasemap -JOc0/0/282/49/6.5 -R$ha_ll[1]/$ha_ll[2]/$ha_ur[1]/$ha_ur[2]r -Ba100.0f7.5mWSne -Y2.0 -K -P >! $file
#
# add the state outlines
#
psxy $cpath/caoutline.asc -JO -R -Ms -W5/$black -: -O -K >> $file
#
# draw the Bridges
#
psxy $data/bridge.asc -JO -R -MB -O -K -W15/175 -: >> $file
#
# add the faults
#
# - old faults
#psxy $cpath/cafaults.asc -JO -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 -JO -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 -JO -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 -JO -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 -JO -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 -JO -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 -JO -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 -JO -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 -JO -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 -JO -R -H2 -M$ -W3/0/0/0to -: -O -K >> $file
#
# identify some towns
#
#-122.350 37.910 10 0 6 2 Richmond
#-121.958 37.522 10 0 6 7 Fremont
#-122.067 37.647 10 0 6 7 Hayward
#-122.267 37.780 10 0 6 2 Oakland
#-121.900 37.330 10 0 6 7 San Jose
#-121.975 37.763 10 0 6 7 San Ramon
#-122.280 37.880 10 0 6 7 Berkeley
#-122.230 38.113 10 0 6 9 Vallejo
#
psxy -JO -R -St0.10 -L -W2/$black -O -K >> $file << eof
-122.350 37.910 
-121.958 37.522 
-122.067 37.647 
-122.267 37.780 
-121.900 37.330 
-121.975 37.763 
-122.280 37.880 
-122.230 38.113 
eof
pstext -JO -R -O -K >> $file << eof
-122.350 37.930 10 0 6 2 Richmond 
-121.968 37.522 10 0 6 7 Fremont 
-122.077 37.647 10 0 6 7 Hayward 
-122.267 37.800 10 0 6 2 Oakland 
-121.910 37.330 10 0 6 7 San Jose
-121.985 37.763 10 0 6 7 San Ramon 
-122.280 37.900 10 0 6 2 Berkeley 
-122.220 38.113 10 0 6 9 Vallejo 
eof
#
# plot seismicity
#
cat search.out.1998 | awk '{print $7, $8, $10*0.035+0.02}' >! temp
psxy temp -JO -R -Sc -L -W3/$blue -: -O -K >> $file
#
# key
#
cat search.eq.key | awk '{print $7, $8, $10*0.035+0.02}' >! temp
psxy temp -Jl-122.2/37.4/38.2/38/1:650000 -R-122.82/-121.6/37.80/38.30 -Sc -L -W3/$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
#
# title
#
pstext -R0/8.5/0/11 -JX8.5/11 -N -O << eof >>  $file
3.5 2 18 0.0 25 6 Jan 1998 Earthquakes
eof
#
#
ghostview $file &