#!/bin/sh
O="sb5.ps"
D=`sed -n 'p' ./pathfile`
YEARS="3259 6076 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97"
rm .gmtcommands .gmtdefaults
gmtset	MEASURE_UNIT inch		#PAPER_WIDTH 8.5
gmtset	BASEMAP_TYPE PLAIN		DEGREE_FORMAT 2
gmtset	LABEL_FONT_SIZE 6		FRAME_PEN 3
gmtset	ANOT_FONT_SIZE 6		ANOT_OFFSET 0.075
gmtset	TICK_PEN 3			TICK_LENGTH -0.08
cat 					<< END			> cpt
B	215	205	155
60	215	205	155	400	100	140	100
F	100	140	100
END

F="-R-120:30/-119:15/33:45/35:00 -JM6.5 -P -V -K"
if [ ! -f $O.tmp ]; then
grdimage $F -X1 -Y2 -Ccpt -I$D/topo/06i.grd $D/topo/06.grd	> $O.tmp; fi
cat					$O.tmp			> $O
pscoast	$F -O -S225/255/255 -Df -W3				>> $O
psxy	$F -O -W6 -M			$D/cdmg/*		>> $O
#psxy	$F -O -W1 -M			$D/scal_roads.xy	>> $O
for	YEAR in $YEARS; do
echo psxy: working on file scsn19$YEAR
awk	'$4<=4.5 && $4>=1.0 {print $1,$2,$3*2}'	$D/eqs/scsn19$YEAR	|
psxy	$F -O -W2/255/0/0 -Sc					>> $O
awk	'$4>4.5  {print $1,$2,$3*2}'	$D/eqs/scsn19$YEAR	|
psxy	$F -O -W8/255/0/0 -Sc				>> $O
awk	'$4>4.5  {print $1,$2,$3*2}'	$D/eqs/scsn19$YEAR	|
psxy	$F -O -W1/0/0/0 -Sc				>> $O; done
psbasemap $F -O -B15mf1mNESW -Lfx3.25/-.8/34.25/30		>> $O

F="-R0/11/0/11 -Jx1 -V -K -O"
pstext	$F -Y-1				<< END			>> $O	# title
3.25 0.60 12 0  0 2 Seismicity of the Santa Barbara Region 1932-1997
0.10 1.10  8 0 30 1 Geoffrey P. Ely
0.90 1.10  8 0 30 1 regenerated by J.Mount 12/12/05
5.65 0.00  8 0  0 2 3
5.80 0.00  8 0  0 2 4
5.95 0.00  8 0  0 2 5
6.10 0.00  8 0  0 2 6
5.90 0.60  8 0  0 2 Magnitude Scale
0.34 0.35 14 0  0 6 N
END
psxy	$F -W2/0/0/0 -Sc		<< END			>> $O
5.65 0.174 0.048
5.80 0.198 0.097
5.95 0.247 0.194
6.10 0.344 0.388
END
psxy 	$F -W3 -L -G0  			<< END			>> $O
0.34 0.52
0.38 0.50
0.36 0.58
0.34 0.70
0.34 0.02
0.30 0.00
0.31 0.18
0.34 0.20
END
cat					$D/ics.def		>> $O
psxy	$F -Sc.7 << END | sed 's/C4/ICSbox/'			>> $O
6.4 8.73
END

F="-R-125/-113.5/32.5/43 -JM.7 -V -O"
psxy	$F -K -W3 -L -G255 -X.6		$D/cal_rough.xy		>> $O
pstext	$F -K				<< END			>> $O
-120.4 38.7 9 -70 0 11 Cal
-120.4 38.8 9 -45 0 9 ifornia
END
psxy	$F -W3 -L			<< END			>> $O
-119.25	33.75
-119.25	35.00
-120.50	35.00
-120.50	33.75
END