#!/bin/sh
O="gv.ps"
D=`sed -n 'p' ./pathfile`
rm .gmtdefaults .gmtcommands
gmtset	MEASURE_UNIT inch		#PAPER_WIDTH 8.5
gmtset	HSV_MIN_VALUE .2		HSV_MAX_VALUE 1
gmtset	BASEMAP_TYPE PLAIN		DEGREE_FORMAT 2
gmtset	LABEL_FONT_SIZE 6		FRAME_PEN 3
gmtset	ANOT_FONT_SIZE 12		ANOT_OFFSET 0.075
gmtset	TICK_PEN 3			TICK_LENGTH 0.08
#gmtset	PAGE_ORIENTATION landscape
cat 					<< END				> cpt
B	0	100	100
10	0	100	100	450	205	255	225
450	255	255	225	1000	255	225	155
1000	255	225	155	1500	150	100	50
1500	100	50	0	3500	255     255     255
END

F="-R-117:30/-116/33:00/34:30 -Js-117/90/115/60 -V -E145/30 -Jz0.000217 -K"
F="-R-117:30/-116/33:00/34:30/-500/3500 -Jm4 -V -E145/30 -Jz0.000217 -K"
#grdcut $D/topo/06i2.grd -R-117:30/-116/33:00/34:30 -G$D/topo/gvdai.grd -V
if [ ! -f $O.tmp ]; then
grdview $F -X1 -Y2 -Qi150 -Z-500 -N-500/0/0/0 -Ccpt $D/topo/06.grd -I$D/topo/gvdai.grd > $O.tmp; fi
cat					$O.tmp				> $O
psbasemap $F -O -Bf1ma10mESZ -Z0					>> $O	# the scale bars
pscoast $F -O -Df -S150/205/255 -Y.09					>> $O
psxyz	$F -O -G255/0/0 -Su.08 -L -Z-500 -Y-.09	<< END			>> $O
-116.673 33.669 1310
END

F="-R0/11/0/11 -Jx1 -V -O -N -P"
cat					$D/ics.def			>> $O
psxy	$F -K -Sc1 << END | sed 's/C4/ICSbox/'				>> $O	# the ICS logo
9 0.9 10
END
pstext	$F -K				<< END				>> $O	# title
0 4.7 18 0  1 1 Garner Valley Downhole
0 4.4 18 0  1 1 Seismic Array
7.5 0.0  8 0 30 3 Geoffrey P. Ely
9 0.0  8 0 30 3 regenerated by J.Mount 1/2/06
0 4.2  9 0  0 1 Vertical Exageration 20x
END
pstext	$F -K -S2			<< END				>> $O
4.65 2.4 10 0 1 5 GVDSA
END
pstext	$F -G255			<< END				>> $O
4.65 2.4 10 0 1 5 GVDSA
END