#!/bin/sh
O="g.ps"
D=`sed -n 'p' ./pathfile`
rm .gmtdefaults .gmtcommands
gmtset	MEASURE_UNIT inch		#PAPER_WIDTH 8.5
gmtset	LABEL_FONT_SIZE 6		FRAME_PEN 3
gmtset	ANOT_FONT_SIZE 6		ANOT_OFFSET 0.08
gmtset	TICK_PEN 1			TICK_LENGTH -0.1
gmtset	GRID_PEN 1			GRID_CROSS_SIZE 0
gmtset	DOTS_PR_INCH 300
gmtset	HEADER_FONT_SIZE 6
cat 					<< END		> cpt
B	100	100	255
-11000	185	120	255	-6000	185	150	255
-6000	185	150	255	-5000	150	200	255
-5000	150	200	255	-2000	225	255	255
-2000	225	255	255	-200	225	255	255
-200	225	255	255	0	255	255	255
0	120	180	180	100	165	220	220
100	120	180	180	400	205	255	225
400	255	255	225	1000	255	225	155
1000	255	225	155	1500	150	100	50
1500	100	50	0	2000	180     155     130
2000	180	155	130	4000	255     255     255
4000	255	255	255	9000	255	255	255
END

F="-R0/360/-90/90 -JG-105/45/7 -V -K -P"
if [ ! -f $O.tmp ]; then
grdimage $F -X.75 -Y2.5 -Ccpt -I$D/topo/tbase15mi.grd	$D/topo/tbase15m.grd 	> $O.tmp; fi
cat							$O.tmp			> $O
#psbasemap $F -Bnesw -X.75 -Y2.5 -P						> $O
psscale	-O -V -K -Ccpt -D3.5/-.15/6.5/.15h -I -L -B:."Elevation in Meters":	>> $O
pscoast	$F -O -Dl -A1000 -W1 -B15g15						>> $O
YEARS="1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997"
YEARS="1995 1996 1997"
for YEAR in $YEARS; do
awk '$4>3.0&& $5<70		{$3=$3/4;print}'	$D/eqs/cnss$YEAR |
psxy	$F -O -Sc -W6/255/255/0							>> $O
awk '$4>3.0&& $5>=70&& $5<300	{$3=$3/4;print}'	$D/eqs/cnss$YEAR |
psxy	$F -O -Sc -W6/255/180/0							>> $O
awk '$4>3.0&& $5>=300		{$3=$3/4;print}'	$D/eqs/cnss$YEAR |
psxy	$F -O -Sc -W6/255/0/0							>> $O
done
awk '$3=="D1" || $3=="D2"'				$D/vol/gvp.xy |
psxy	$F -O -W1/150 -St.05							>> $O
awk '$3=="D3" || $3=="D4" || $3=="D5" || $3=="D6"'	$D/vol/gvp.xy |
psxy	$F -O -W1/150 -St.025							>> $O
awk '$3=="D7" || $3=="U" || $3=="Q" || $3=="?"'		$D/vol/gvp.xy |
psxy	$F -O -W1/150 -St.025							>> $O
psbasemap $F -O -Bnesw								>> $O

F="-R0/11/0/11 -Jx1 -V -O"
pstext	$F -K -X-.75 -Y-2.5				<< END			>> $O	# title
4.25 9.95 18 0 31 6 NORTH AMERICAN SEISMICITY 1987-1997
7.00 0.88  8 0 30 3 Geoffrey P. Ely 08/22/97
7.25 0.75  8 0 30 3 regenerated by J. Mount 12/05/05
0.75 0.88  7 0  0 1 Topographic data provided by the National Geophysical Data Center.
0.75 0.75  7 0  0 1 Eathquake locations provided by the Council of the National Seismic System.
1.50 1.35  8 0  0 2 Earthquake Depth in Km
1.00 1.55  8 0  0 2 <70
1.50 1.55  8 0  0 2 70-300
2.00 1.55  8 0  0 2 >300
7.00 1.35  8 0  0 2 Earthquake Magnitude
4.20 1.750  8 0  0 2 Elevation in Meters
6.50 1.55  8 0  0 2 4
6.75 1.55  8 0  0 2 5
7.00 1.55  8 0  0 2 6
7.25 1.55  8 0  0 2 7
7.50 1.55  8 0  0 2 8
END
psxy	$F -K -W6/255/255/0 -Sc				<< END			>> $O
1.00 1.797 0.194
END
psxy	$F -K -W6/255/180/0 -Sc				<< END			>> $O
1.50 1.797 0.194
END
psxy	$F -K -W6/255/0/0 -Sc				<< END			>> $O
2.00 1.797 0.194
6.50 1.706 0.012
6.75 1.712 0.024
7.00 1.724 0.048
7.25 1.748 0.097
7.50 1.797 0.194
END
psxy	$F -K -W3 -L					<< END			>> $O
0.55  0.55
0.55 10.45
7.95 10.45
7.95  0.55
END
cat							$D/ics.def		>> $O
psxy	$F -Sc.5 << END | sed 's/C4/ICSbox/'					>> $O
7.75 1.13
END