#!/bin/sh
O="sbc3.ps"
D=`sed -n 'p' ./pathfile`
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
gmtset	DOTS_PR_INCH 300
cat 					<< END			> cpt
-1	150	150	150	-.1	225	225	225
-.1	225	225	225	.8	255	255	255
END

F="-R-121/-119/33:45/34:45 -JM6.5 -V -K -P"
if [ ! -f $O.tmp ]; then
grdimage $F -X1 -Y2 $D/topo/06i.grd -Ccpt			> $O.tmp; fi
cat					$O.tmp			> $O
pscoast	$F -O -S255 -Df -Na/3 -W3				>> $O
psxy	$F -O -L -Sc.02 		$D/sta/wells		>> $O
psbasemap $F -O -Bf1ma15mNESW -Lfx5.3/-0.6/34.25/30		>> $O

F="-R0/11/0/11 -Jx1 -V -K -O"
pstext	$F -Y-1				<< END			>> $O # title
0 0.54 10 0 29 1 Figure 3.  Map of the Santa Barbara Channel and western
0 0.36 10 0 29 1 Ventura Basin showing coverage of seismic relocation,
0 0.18 10 0 29 1 seismic refraction surveys, and offshore well data.
6.24 0.35 14 0 0 6 N
END
psxy 	$F -W3 -L -G0			<< END			>> $O
6.24 0.52
6.28 0.50
6.26 0.58
6.24 0.70
6.24 0.02
6.20 0.00
6.21 0.18
6.24 0.20 
END
cat					$D/ics.def		>> $O
psxy	$F -Sc.5 << END | sed 's/C4/ICSbox/'		>> $O
8.85 6.25
END

F="-R-125/-113.5/32.5/43 -JM.7 -V -O"
psxy	$F -K -W3 -L -G255 -X3.9	$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	33.75
-119	34.75
-121	34.75
-121	33.75
END