#!/bin/csh

#  Draw basemap
#  ------------
gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica

psbasemap -JX6/1.5 \
-R1992.5/2003.5/-30/20 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"North (mm)":WeSn  \
> silv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> silv.ps
1992.6713	-16.8	7.1
1992.6740	-20.2	6.6
1999.1328	-3.4	7.2
1999.1356	-11.9	6.8
1999.1383	-10.4	7.3
1999.1410	-9.0	7.0
1999.7132	1.7	6.8
1999.7159	0.2	6.5
1999.7378	2.8	6.3
2003.4696	14.7	6.4
2003.4723	16.7	6.6
2003.4805	17.6	6.7
2003.4833	18.1	6.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> silv.ps
1992.6713	-16.8	7.1
1992.6740	-20.2	6.6
1999.1328	-3.4	7.2
1999.1356	-11.9	6.8
1999.1383	-10.4	7.3
1999.1410	-9.0	7.0
1999.7132	1.7	6.8
1999.7159	0.2	6.5
1999.7378	2.8	6.3
2003.4696	14.7	6.4
2003.4723	16.7	6.6
2003.4805	17.6	6.7
2003.4833	18.1	6.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> silv.ps
SBreak	Segment0	Comp0
 1992.6713	 -30.1779	Segment0	Comp0
 2003.4833	  13.6850	Segment0	Comp0
EOT


#  Draw basemap
#  ------------
gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica

psbasemap -JX6/1.5 \
-R1992.5/2003.5/-10/30 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"East (mm)":WeSn  \
-O >> silv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> silv.ps
1992.6713	-0.4	8.0
1992.6740	21.3	7.2
1999.1328	1.7	7.8
1999.1356	0.8	7.3
1999.1383	-0.3	9.8
1999.1410	10.6	8.7
1999.7132	-1.1	8.0
1999.7159	1.8	7.5
1999.7378	3.7	7.4
2003.4696	-3.9	7.4
2003.4723	-2.8	7.9
2003.4805	-8.6	7.9
2003.4833	-6.6	7.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> silv.ps
1992.6713	-0.4	8.0
1992.6740	21.3	7.2
1999.1328	1.7	7.8
1999.1356	0.8	7.3
1999.1383	-0.3	9.8
1999.1410	10.6	8.7
1999.7132	-1.1	8.0
1999.7159	1.8	7.5
1999.7378	3.7	7.4
2003.4696	-3.9	7.4
2003.4723	-2.8	7.9
2003.4805	-8.6	7.9
2003.4833	-6.6	7.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> silv.ps
SBreak	Segment0	Comp1
 1992.6713	  13.6159	Segment0	Comp1
 2003.4833	  -4.7457	Segment0	Comp1
EOT


#  Draw basemap
#  ------------
gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica

psbasemap -JX6/1.5 \
-R1992.5/2003.5/-30/40 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"Up (mm)":WeSn  \
-O >> silv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> silv.ps
1992.6713	-21.6	22.9
1992.6740	35.0	21.9
1999.1328	0.7	23.5
1999.1356	-8.9	21.5
1999.1383	-25.2	23.5
1999.1410	-7.3	26.4
1999.7132	4.0	23.6
1999.7159	4.5	22.3
1999.7378	-10.0	22.8
2003.4696	30.1	22.6
2003.4723	11.5	24.1
2003.4805	1.6	24.1
2003.4833	7.5	21.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> silv.ps
1992.6713	-21.6	22.9
1992.6740	35.0	21.9
1999.1328	0.7	23.5
1999.1356	-8.9	21.5
1999.1383	-25.2	23.5
1999.1410	-7.3	26.4
1999.7132	4.0	23.6
1999.7159	4.5	22.3
1999.7378	-10.0	22.8
2003.4696	30.1	22.6
2003.4723	11.5	24.1
2003.4805	1.6	24.1
2003.4833	7.5	21.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> silv.ps
SBreak	Segment0	Comp2
 1992.6713	  -7.9702	Segment0	Comp2
 2003.4833	   3.6277	Segment0	Comp2
EOT