#!/bin/sh
O="nr6.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
B	215	205	155
60	215	205	155	400	100	140	100
F	100	140	100
END

F="-R-119/-118/34:00/34:30 -JM9 -V -K"
if [ ! -f $O.tmp ]; then
grdimage $F -X1 -Y2 $D/topo/03.grd -I$D/topo/03i.grd -Ccpt	> $O.tmp; fi
cat					$O.tmp			> $O
echo	"2 setlinejoin"						>> $O
pscoast	$F -O -W3 -S225/255/255 -Df -W3				>> $O
psxy	$F -O -W1 -M			$D/scal_roads.xy	>> $O
psxy	$F -O -W8/0 -M			$D/cal_faults.xy	>> $O
psxy	$F -O -L -St.13 -G255/0/0	$D/sta/usgs		>> $O
psxy	$F -O -L -St.13 -G255/0/0	$D/sta/nort		>> $O
psxy	$F -O -L -Sd.10 -G255/0/255	$D/sta/nr_scsn		>> $O
psxy	$F -O -L -Ss.10 -G0/0/255	$D/sta/nr_terr		>> $O
pstext	$F -O	 			$D/sta/nort.txt		>> $O
pstext	$F -O				$D/sta/nr_scsn.txt	>> $O
pstext	$F -O				$D/sta/nr_terr.txt	>> $O
psxy	$F -O -W6/255/255/0 -L -Sc	$D/amp/usgs2		>> $O
psxy	$F -O -W6/255/255/255 -L -Sc	$D/amp/rs		>> $O
psbasemap $F -O -Bf1ma10mNESW -Lfx4.5/-.8/34.25/20		>> $O

F="-R0/11/0/11 -Jx1 -V -K -O"
pstext	$F -N -Y-1			<< END			>> $O	# title
4.50 0.60 12 0  0 2 Response Spectral Amplification Factors
4.50 0.40 12 0  0 2 at 1.5hz: USGS (yellow) and UCSB (white)
0.10 1.10  8 0 30 1 Geoffrey P. Ely
0.95 1.10  8 0 30 1 regenerated by J.Mount 12/12/05
6.90 0.51  9 0  0 5 Northridge portable deployment
6.90 0.35  9 0  0 5 Southern California Seismic Network
6.90 0.19  9 0  0 5 Terrascope stations
0.75 0.50  9 0  0 2 Amplification
0.45 0.00  8 0  0 2 1
0.60 0.00  8 0  0 2 2
0.75 0.00  8 0  0 2 4
0.90 0.00  8 0  0 2 5
1.05 0.00  8 0  0 2 6
1.54 0.35 14 0  0 6 N
END
psxy	$F -L -St.13 -G255/0/0		<< END			>> $O
6.75 0.51
END
psxy	$F -L -Sd.10 -G255/0/255	<< END			>> $O
6.75 0.35
END
psxy	$F -L -Ss.10 -G0/0/255		<< END			>> $O
6.75 0.19
END
psxy	$F -L -W3 -Sc			<< END			>> $O
0.45 0.180 0.10
0.60 0.205 0.15
0.75 0.230 0.20
0.90 0.255 0.25
1.05 0.280 0.30
END
psxy 	$F -W3 -L -G0 	 		<< END			>> $O
1.54 0.52
1.58 0.50
1.56 0.58
1.54 0.70
1.54 0.02
1.50 0.00
1.51 0.18
1.54 0.20 
END
cat					$D/ics.def		>> $O
psxy	$F -Sc1		<< END | sed 's/C4/ICSbox/'		>> $O
8.9 6.3
END

F="-R-125/-113.5/32.5/43 -JM.7 -V -O"
psxy	$F -K -W3 -L -G255 -X1.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
-118	34
-118	34.5
-119	34.5
-119	34
END