#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> kels.ps
1998.6454	-30.0	7.1
1998.6536	-21.8	6.6
1998.6563	-16.1	7.2
1998.8180	-29.2	6.8
1998.8207	-22.7	7.3
1999.6255	-14.7	7.0
2000.5510	-6.0	6.8
2000.5537	-9.5	6.5
2000.5564	-12.7	6.3
2000.5592	-4.4	6.4
2000.5619	-8.3	6.6
2000.5646	-7.6	6.7
2000.5674	-6.8	6.1
2000.5701	-10.4	6.7
2000.5729	-9.1	6.1
2000.5756	-13.5	6.1
2000.5783	-12.0	6.5
2000.5811	-4.6	6.5
2000.5838	-5.7	4.8
2000.7043	-6.2	6.9
2000.7070	-7.7	7.3
2002.5250	3.9	7.0
2002.5277	2.8	6.8
2003.6338	12.5	6.5
2003.6366	11.7	6.7
2003.6393	14.4	6.0
2003.6421	11.3	6.9
2003.6448	14.1	6.6
2003.6475	12.8	6.5
2003.6503	12.2	6.5
2003.6530	10.4	6.5
2003.6557	13.3	6.6
2003.7707	11.4	6.3
2003.7735	6.4	6.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> kels.ps
1998.6454	-30.0	7.1
1998.6536	-21.8	6.6
1998.6563	-16.1	7.2
1998.8180	-29.2	6.8
1998.8207	-22.7	7.3
1999.6255	-14.7	7.0
2000.5510	-6.0	6.8
2000.5537	-9.5	6.5
2000.5564	-12.7	6.3
2000.5592	-4.4	6.4
2000.5619	-8.3	6.6
2000.5646	-7.6	6.7
2000.5674	-6.8	6.1
2000.5701	-10.4	6.7
2000.5729	-9.1	6.1
2000.5756	-13.5	6.1
2000.5783	-12.0	6.5
2000.5811	-4.6	6.5
2000.5838	-5.7	4.8
2000.7043	-6.2	6.9
2000.7070	-7.7	7.3
2002.5250	3.9	7.0
2002.5277	2.8	6.8
2003.6338	12.5	6.5
2003.6366	11.7	6.7
2003.6393	14.4	6.0
2003.6421	11.3	6.9
2003.6448	14.1	6.6
2003.6475	12.8	6.5
2003.6503	12.2	6.5
2003.6530	10.4	6.5
2003.6557	13.3	6.6
2003.7707	11.4	6.3
2003.7735	6.4	6.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> kels.ps
SBreak	Segment0	Comp0
 1998.6454	 -22.3528	Segment0	Comp0
 2003.7735	  12.8493	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1998.5/2004/-20/10 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"East (mm)":WeSn  \
-O >> kels.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> kels.ps
1998.6454	-16.8	8.0
1998.6536	-15.5	7.2
1998.6563	-15.4	7.8
1998.8180	-15.5	7.3
1998.8207	-2.6	9.8
1999.6255	-13.1	8.7
2000.5510	-7.1	8.0
2000.5537	-7.6	7.5
2000.5564	-3.6	7.4
2000.5592	-8.3	7.4
2000.5619	-9.4	7.9
2000.5646	-9.0	7.9
2000.5674	-6.8	7.1
2000.5701	-9.3	8.0
2000.5729	-6.1	7.0
2000.5756	-2.7	7.1
2000.5783	-9.4	7.1
2000.5811	-7.3	7.3
2000.5838	-11.0	5.5
2000.7043	-5.1	8.9
2000.7070	-2.9	9.9
2002.5250	7.5	9.3
2002.5277	1.6	7.3
2003.6338	9.5	7.0
2003.6366	7.9	7.3
2003.6393	9.2	9.0
2003.6421	8.5	9.6
2003.6448	8.7	7.0
2003.6475	6.8	6.9
2003.6503	8.7	7.0
2003.6530	8.8	7.1
2003.6557	6.5	7.1
2003.7707	-1.2	7.0
2003.7735	2.4	6.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> kels.ps
1998.6454	-16.8	8.0
1998.6536	-15.5	7.2
1998.6563	-15.4	7.8
1998.8180	-15.5	7.3
1998.8207	-2.6	9.8
1999.6255	-13.1	8.7
2000.5510	-7.1	8.0
2000.5537	-7.6	7.5
2000.5564	-3.6	7.4
2000.5592	-8.3	7.4
2000.5619	-9.4	7.9
2000.5646	-9.0	7.9
2000.5674	-6.8	7.1
2000.5701	-9.3	8.0
2000.5729	-6.1	7.0
2000.5756	-2.7	7.1
2000.5783	-9.4	7.1
2000.5811	-7.3	7.3
2000.5838	-11.0	5.5
2000.7043	-5.1	8.9
2000.7070	-2.9	9.9
2002.5250	7.5	9.3
2002.5277	1.6	7.3
2003.6338	9.5	7.0
2003.6366	7.9	7.3
2003.6393	9.2	9.0
2003.6421	8.5	9.6
2003.6448	8.7	7.0
2003.6475	6.8	6.9
2003.6503	8.7	7.0
2003.6530	8.8	7.1
2003.6557	6.5	7.1
2003.7707	-1.2	7.0
2003.7735	2.4	6.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> kels.ps
SBreak	Segment0	Comp1
 1998.6454	 -15.9076	Segment0	Comp1
 2003.7735	   8.1855	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1998.5/2004/-110/20 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"Up (mm)":WeSn  \
-O >> kels.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> kels.ps
1998.6454	-104.4	22.9
1998.6536	-41.1	21.9
1998.6563	-107.4	23.5
1998.8180	-5.6	21.5
1998.8207	-6.9	23.5
1999.6255	-3.1	26.4
2000.5510	2.3	23.6
2000.5537	1.1	22.3
2000.5564	-11.1	22.8
2000.5592	-5.5	22.6
2000.5619	-18.6	24.1
2000.5646	-5.6	24.1
2000.5674	1.9	21.3
2000.5701	-13.5	23.1
2000.5729	-0.4	20.7
2000.5756	8.2	21.0
2000.5783	-0.1	21.3
2000.5811	9.1	20.8
2000.5838	-0.3	15.9
2000.7043	5.6	27.3
2000.7070	-7.7	32.2
2002.5250	-1.5	27.6
2002.5277	-1.7	21.9
2003.6338	-1.9	21.1
2003.6366	7.0	21.7
2003.6393	4.4	21.3
2003.6421	6.8	24.5
2003.6448	6.7	21.8
2003.6475	5.7	21.1
2003.6503	-5.3	20.9
2003.6530	2.3	21.4
2003.6557	-5.2	21.4
2003.7707	13.2	21.0
2003.7735	5.9	21.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> kels.ps
1998.6454	-104.4	22.9
1998.6536	-41.1	21.9
1998.6563	-107.4	23.5
1998.8180	-5.6	21.5
1998.8207	-6.9	23.5
1999.6255	-3.1	26.4
2000.5510	2.3	23.6
2000.5537	1.1	22.3
2000.5564	-11.1	22.8
2000.5592	-5.5	22.6
2000.5619	-18.6	24.1
2000.5646	-5.6	24.1
2000.5674	1.9	21.3
2000.5701	-13.5	23.1
2000.5729	-0.4	20.7
2000.5756	8.2	21.0
2000.5783	-0.1	21.3
2000.5811	9.1	20.8
2000.5838	-0.3	15.9
2000.7043	5.6	27.3
2000.7070	-7.7	32.2
2002.5250	-1.5	27.6
2002.5277	-1.7	21.9
2003.6338	-1.9	21.1
2003.6366	7.0	21.7
2003.6393	4.4	21.3
2003.6421	6.8	24.5
2003.6448	6.7	21.8
2003.6475	5.7	21.1
2003.6503	-5.3	20.9
2003.6530	2.3	21.4
2003.6557	-5.2	21.4
2003.7707	13.2	21.0
2003.7735	5.9	21.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> kels.ps
SBreak	Segment0	Comp2
 1998.6454	  -6.7243	Segment0	Comp2
 2003.7735	   3.7244	Segment0	Comp2
EOT