#!/bin/sh
O="cal.ps"
D=`sed -n 'p' ./pathfile`
DATE=`date '+%m/%d/%y'`
rm .gmtdefaults .gmtcommands
gmtset	MEASURE_UNIT inch		#PAPER_WIDTH 8.5
gmtset	BASEMAP_TYPE PLAIN		DEGREE_FORMAT 2
gmtset	LABEL_FONT_SIZE 6		FRAME_PEN 1
gmtset	ANOT_FONT_SIZE 6		ANOT_OFFSET 0.075
gmtset	TICK_PEN 1			TICK_LENGTH -0.05
gmtset  ANOT_MIN_ANGLE 45
cat 					<< END			> cpt
-1	100	100	100	-.01	178	178	178
-.01	178	178	178	1	255	255	255
END
cat					<< END			> cpt
-30000	0	0	0	-50	200	200	200
-50	200	200	200	50	230	230	230
50	230	230	230	30000	255	255	255
END
cat					<< END			> cpt
B	200	200	200
1	200	200	200	20	210	210	210
20	210	210	210	75	220	220	220
F	220	220	220
END
CONTS="	0400 0800 1200 1600 2000 2400 2800 3200 3600 4000 4400 4800 5200"

#grdcut -V -R-126/-113/32/43		$D/topo/usa30.grd	-G$D/topo/cal30.grd
#grdgradient -V -Nt -A270		$D/topo/cal30.grd	-G$D/topo/cal30i.grd
#grdcut -V -R-126/-113/32/43		$D/topo/usa60.grd	-G$D/topo/cal60.grd
#grdgradient -V -Nt -A270		$D/topo/cal60.grd	-G$D/topo/cal60i.grd
#grdmath -V $D/topo/cal60i.grd SIGN $D/topo/cal60i.grd 25000 DIV ABS SQRT MUL = $D/topo/cal60i2.grd
#grdmath -V $D/topo/cal60i.grd .2 ADD 1.25 MUL = $D/topo/cal60i2.grd
#grdmath -V .25 $D/topo/cal60i.grd ABS .5 SUB 2 POW SUB $D/topo/cal60i.grd ADD = $D/topo/cal60i2.grd
F="-R-125:30/-113:30/32:00/42:45 -JM7.5 -V -K -P"
F="-R-125/32:01/-113:15/42:30r -JL-119.5/38/35/40/7.2 -V -K -P"
if [ ! -f $O.tmp ]; then
grdimage $F -O -Ccpt -I$D/topo/cal60i2.grd $D/topo/cal60.grd	> $O.tmp; fi
psclip	$F -X.65 -Y1.65 -M $D/ca_nv.xy | sed ' 9s/^/.9999 .9999 scale/' > $O
echo "2 setlinejoin"						>> $O
cat					$O.tmp			>> $O
psclip	-C -O -K						>> $O
awk '$3=="c"' $D/ca_nv.xy | psxy $F -O -W1 -M			>> $O
psxy	$F -O -L -W1 -M -G210		$D/ca_isl.xy		>> $O
awk '$3=="b"' $D/ca_nv.xy | psxy $F -O -W1t5_5_25_5_5_5:0	>> $O
psxy	$F -O -W1t5_5_25_5_5_5:0	$D/ca_nv_border.xy	>> $O
for CONT in $CONTS
do
psxy	$F -K -O -H34 -M# -: -W1/200	$D/bathy/eez/pbty$CONT.mgf	>> $O
done
psbasemap $F -O -Lfx5.27/-.3/37/200				>> $O	# the km scalebar
gmtset	ANOT_OFFSET -0.05
psbasemap $F -O -Bf1a2NSew					>> $O	# the longitude degrees
gmtset	ANOT_OFFSET -0.05
psbasemap $F -O -Bf1a2nsEW | sed '/2156 284 M V -86.7/d'	>> $O	# the latitude degress
psxy	$F -O -M -W4/0/105/255		$D/cdmg/q		>> $O
psxy	$F -O -M -W4/205/0/205		$D/cdmg/lq		>> $O
psxy	$F -O -M -W4/255/185/0		$D/cdmg/holo		>> $O
psxy	$F -O -M -W4/255/0/0		$D/cdmg/hist		>> $O
awk '$4<=5.5' $D/eqs/cnss19* | psxy $F -O -W2/255/0/0 -Sc	>> $O
awk '$4> 5.5' $D/eqs/cnss19* | psxy $F -O -W5/255/255/0 -Sc	>> $O

F="-R0/11/0/11 -Jx1 -V -O -N"
pstext	$F -K -Y-1			<< END			>> $O	# title
3.60 9.40 16 0 29 2 EARTHQUAKES AND FAULTS IN CALIFORNIA AND NEVADA
1.00 1.40  7 0 30 1 Geoffrey P. Ely 08/22/97
1.00 1.30  7 0 30 1 regenerated by J.Mount $DATE
0.25 0.70  8 0  0 1 Age of Faulting
0.60 0.55  8 0  0 5 Historical
0.60 0.40  8 0  0 5 Holocene
0.60 0.25  8 0  0 5 Late Quaternary
0.60 0.10  8 0  0 5 Quaternary
2.15 0.05  8 0  0 2 4
2.40 0.05  8 0  0 2 5
2.65 0.05  8 0  0 2 6
2.90 0.05  8 0  0 2 7
2.00 0.70  8 0  0 1 Earthquake Magnitude
2.00 0.55  8 0  0 1 1987 to 1997
3.70 0.29  7 0  0 1 Fault map provided by the California Division of Mines and Geology.
3.70 0.17  7 0  0 1 Earthquake locations provided by the Council of the National Seismic System.
3.70 0.05  7 0  0 1 Topographic data provided by the United States Geological Survey.
END
psxy	$F -K -W4/255/0/0		<< END			>> $O
0.25 0.55
0.50 0.55
END
psxy	$F -K -W4/255/185/0		<< END			>> $O
0.25 0.40
0.50 0.40
END
psxy	$F -K -W4/205/0/205		<< END			>> $O
0.25 0.25
0.50 0.25
END
psxy	$F -K -W4/0/105/255		<< END			>> $O
0.25 0.10
0.50 0.10
END
psxy	$F -K -W2 -Sc			<< END			>> $O
2.15 0.224 0.048
2.40 0.248 0.097
END
psxy	$F -K -W5 -Sc			<< END			>> $O
2.65 0.297 0.194
2.90 0.387 0.388
END
cat					$D/ics.def		>> $O
psxy	$F -K -Sc.7 << END	| sed 's/C4/ICSbox/'			>> $O
0.90 1.80
END

psxy	$F -X-.65 -Y-.65 -W3 -L	<< END			>> $O
0.55  0.55
0.55 10.45
7.95 10.45
7.95  0.55
END