#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> olym.ps
1993.6487	-23.3	7.1
1993.6514	-15.8	6.6
1993.6542	-19.7	7.2
1993.6569	-18.1	6.8
1995.5926	-24.6	7.3
1995.5953	-20.7	7.0
1995.5981	-21.4	6.8
1995.6008	-17.2	6.5
1997.6323	-5.9	6.3
1997.6351	-4.3	6.4
1999.6502	-0.1	6.6
1999.6529	8.7	6.7
2001.1753	17.0	6.1
2001.1781	16.6	6.7
2001.1808	16.6	6.1
2001.1835	14.1	6.1
2001.1863	17.1	6.5
2003.6366	30.9	6.5
2003.6393	34.0	4.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> olym.ps
1993.6487	-23.3	7.1
1993.6514	-15.8	6.6
1993.6542	-19.7	7.2
1993.6569	-18.1	6.8
1995.5926	-24.6	7.3
1995.5953	-20.7	7.0
1995.5981	-21.4	6.8
1995.6008	-17.2	6.5
1997.6323	-5.9	6.3
1997.6351	-4.3	6.4
1999.6502	-0.1	6.6
1999.6529	8.7	6.7
2001.1753	17.0	6.1
2001.1781	16.6	6.7
2001.1808	16.6	6.1
2001.1835	14.1	6.1
2001.1863	17.1	6.5
2003.6366	30.9	6.5
2003.6393	34.0	4.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> olym.ps
SBreak	Segment0	Comp0
 1993.6487	 -23.2648	Segment0	Comp0
 2003.6393	  28.0606	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> olym.ps
1993.6487	-21.4	8.0
1993.6514	-21.3	7.2
1993.6542	-20.8	7.8
1993.6569	-13.2	7.3
1995.5926	-12.2	9.8
1995.5953	-20.3	8.7
1995.5981	-13.1	8.0
1995.6008	-22.4	7.5
1997.6323	0.9	7.4
1997.6351	-0.8	7.4
1999.6502	-0.2	7.9
1999.6529	5.1	7.9
2001.1753	16.1	7.1
2001.1781	12.8	8.0
2001.1808	10.8	7.0
2001.1835	6.4	7.1
2001.1863	11.2	7.1
2003.6366	22.4	7.3
2003.6393	25.8	5.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> olym.ps
1993.6487	-21.4	8.0
1993.6514	-21.3	7.2
1993.6542	-20.8	7.8
1993.6569	-13.2	7.3
1995.5926	-12.2	9.8
1995.5953	-20.3	8.7
1995.5981	-13.1	8.0
1995.6008	-22.4	7.5
1997.6323	0.9	7.4
1997.6351	-0.8	7.4
1999.6502	-0.2	7.9
1999.6529	5.1	7.9
2001.1753	16.1	7.1
2001.1781	12.8	8.0
2001.1808	10.8	7.0
2001.1835	6.4	7.1
2001.1863	11.2	7.1
2003.6366	22.4	7.3
2003.6393	25.8	5.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> olym.ps
SBreak	Segment0	Comp1
 1993.6487	 -21.1218	Segment0	Comp1
 2003.6393	  21.7101	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> olym.ps
1993.6487	27.9	22.9
1993.6514	11.2	21.9
1993.6542	11.3	23.5
1993.6569	25.6	21.5
1995.5926	37.5	23.5
1995.5953	20.5	26.4
1995.5981	9.4	23.6
1995.6008	2.6	22.3
1997.6323	0.0	22.8
1997.6351	-2.5	22.6
1999.6502	8.1	24.1
1999.6529	0.3	24.1
2001.1753	14.4	21.3
2001.1781	-14.8	23.1
2001.1808	-19.5	20.7
2001.1835	-21.4	21.0
2001.1863	-16.7	21.3
2003.6366	-9.7	20.8
2003.6393	-4.5	15.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> olym.ps
1993.6487	27.9	22.9
1993.6514	11.2	21.9
1993.6542	11.3	23.5
1993.6569	25.6	21.5
1995.5926	37.5	23.5
1995.5953	20.5	26.4
1995.5981	9.4	23.6
1995.6008	2.6	22.3
1997.6323	0.0	22.8
1997.6351	-2.5	22.6
1999.6502	8.1	24.1
1999.6529	0.3	24.1
2001.1753	14.4	21.3
2001.1781	-14.8	23.1
2001.1808	-19.5	20.7
2001.1835	-21.4	21.0
2001.1863	-16.7	21.3
2003.6366	-9.7	20.8
2003.6393	-4.5	15.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> olym.ps
SBreak	Segment0	Comp2
 1993.6487	  17.2413	Segment0	Comp2
 2003.6393	 -22.0150	Segment0	Comp2
EOT