#!/bin/sh
O="lab3.ps"
D=`sed -n 'p' ./pathfile`
DATE=`date '+%m/%d/%y'`
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
-1	110	110	110	-.1	205	205	205
-.1	205	205	205	.8	255	255	255
END

F="-R-118:40/-117:45/33:40/34:20 -JM6.5	-V -K -P"
if [ ! -f $O.tmp ]; then
grdimage $F -X1 -Y3 -Ccpt		$D/topo/03i.grd		> $O.tmp; fi
cat					$O.tmp			> $O
pscoast	$F -O -W3 -S255 -Df					>> $O
psxy	$F -O -W1 -M			$D/scal_roads.xy	>> $O
#psxy	$F -O -W1 -M			$D/ca.xy	>> $O
psxy	$F -O -W8 -M			$D/cal_faults.xy	>> $O
psxy	$F -O -W1 -L -G255 -Sd.05	$D/amp/ralph30		>> $O
psxy	$F -O -W3 -Sc			$D/amp/ralph30		>> $O
psbasemap $F -O -Bf1ma10mNESW -Lfx3.25/-0.8/34/20			>> $O

F="-R0/11/0/11 -Jx1 -V -K -O -N"
pstext	$F -Y-1				<< END			>> $O # title
3.25 0.50 12 0  0 2 Response Spectral Amplification: 3s Period
0.10 1.10  8 0 30 1 Geoffrey P. Ely 
0.85 1.10  8 0 30 1 regenerated by J.Mount  $DATE
5.75 0.50  9 0  0 2 Amplification
5.45 0.00  8 0  0 2 1
5.60 0.00  8 0  0 2 2
5.75 0.00  8 0  0 2 4
5.90 0.00  8 0  0 2 5
6.05 0.00  8 0  0 2 6
0.34 0.35 16 0  0 6 N
END
psxy	$F -L -W3/0 -Sc			<< END			>> $O
5.45 0.180 0.10
5.60 0.205 0.15
5.75 0.230 0.20
5.90 0.255 0.25
6.05 0.280 0.30
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 -Sc1		<< END	| sed 's/C4/ICSbox/'		>> $O
6.40 6.57
END

F="-R-125/-113.5/32.5/43 -JM.7 -V -O"
psxy	$F -K -W3 -L -G255 -X0.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
-117.75	33.67
-117.75	34.33
-118.67	34.33
-118.67	33.67
END