#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gp17.ps
1993.6542	-4.8	7.1
1993.6569	-0.1	6.6
1995.5981	-1.9	7.2
1997.6351	-12.8	6.8
1997.6378	2.4	7.3
1999.6502	17.1	7.0
1999.6529	4.0	6.8
2001.1835	-10.0	6.5
2001.1890	0.7	6.3
2003.6421	-0.8	6.4
2003.6448	21.9	6.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gp17.ps
1993.6542	-4.8	7.1
1993.6569	-0.1	6.6
1995.5981	-1.9	7.2
1997.6351	-12.8	6.8
1997.6378	2.4	7.3
1999.6502	17.1	7.0
1999.6529	4.0	6.8
2001.1835	-10.0	6.5
2001.1890	0.7	6.3
2003.6421	-0.8	6.4
2003.6448	21.9	6.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> gp17.ps
SBreak	Segment0	Comp0
 1993.6542	  -4.2026	Segment0	Comp0
 2003.6448	   6.2780	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2004/-30/20 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"East (mm)":WeSn  \
-O >> gp17.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gp17.ps
1993.6542	-13.1	8.0
1993.6569	0.3	7.2
1995.5981	4.5	7.8
1997.6351	-1.4	7.3
1997.6378	19.7	9.8
1999.6502	6.1	8.7
1999.6529	4.5	8.0
2001.1835	-29.4	7.5
2001.1890	-8.6	7.4
2003.6421	-20.6	7.4
2003.6448	2.3	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gp17.ps
1993.6542	-13.1	8.0
1993.6569	0.3	7.2
1995.5981	4.5	7.8
1997.6351	-1.4	7.3
1997.6378	19.7	9.8
1999.6502	6.1	8.7
1999.6529	4.5	8.0
2001.1835	-29.4	7.5
2001.1890	-8.6	7.4
2003.6421	-20.6	7.4
2003.6448	2.3	7.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> gp17.ps
SBreak	Segment0	Comp1
 1993.6542	   2.6331	Segment0	Comp1
 2003.6448	  -3.1881	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2004/-40/120 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"Up (mm)":WeSn  \
-O >> gp17.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gp17.ps
1993.6542	-32.7	22.9
1993.6569	3.4	21.9
1995.5981	-16.9	23.5
1997.6351	-3.5	21.5
1997.6378	-28.9	23.5
1999.6502	22.4	26.4
1999.6529	-7.8	23.6
2001.1835	-19.7	22.3
2001.1890	-9.2	22.8
2003.6421	112.8	22.6
2003.6448	70.8	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gp17.ps
1993.6542	-32.7	22.9
1993.6569	3.4	21.9
1995.5981	-16.9	23.5
1997.6351	-3.5	21.5
1997.6378	-28.9	23.5
1999.6502	22.4	26.4
1999.6529	-7.8	23.6
2001.1835	-19.7	22.3
2001.1890	-9.2	22.8
2003.6421	112.8	22.6
2003.6448	70.8	24.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> gp17.ps
SBreak	Segment0	Comp2
 1993.6542	 -32.0709	Segment0	Comp2
 2003.6448	  39.2268	Segment0	Comp2
EOT