#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1992.5/2004/-50/50 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"North (mm)":WeSn  \
> reub.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> reub.ps
1992.6521	-42.6	7.1
1992.6548	-40.9	6.6
1992.6576	-44.9	7.2
1993.6405	-35.9	6.8
1995.5817	-35.3	7.3
1995.5844	-25.1	7.0
1995.5871	-27.2	6.8
1996.6166	-11.1	6.5
1996.6193	-14.8	6.3
1996.6221	-12.8	6.4
1996.6248	-16.9	6.6
1998.6344	1.0	6.7
1998.6372	-3.0	6.1
2003.5627	38.8	6.7
2003.5654	43.6	6.1
2003.5681	44.5	6.1
2003.5709	43.0	6.5
2003.5736	45.2	6.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> reub.ps
1992.6521	-42.6	7.1
1992.6548	-40.9	6.6
1992.6576	-44.9	7.2
1993.6405	-35.9	6.8
1995.5817	-35.3	7.3
1995.5844	-25.1	7.0
1995.5871	-27.2	6.8
1996.6166	-11.1	6.5
1996.6193	-14.8	6.3
1996.6221	-12.8	6.4
1996.6248	-16.9	6.6
1998.6344	1.0	6.7
1998.6372	-3.0	6.1
2003.5627	38.8	6.7
2003.5654	43.6	6.1
2003.5681	44.5	6.1
2003.5709	43.0	6.5
2003.5736	45.2	6.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> reub.ps
SBreak	Segment0	Comp0
 1992.6521	 -46.9587	Segment0	Comp0
 2003.5736	  42.6634	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> reub.ps
1992.6521	-11.0	8.0
1992.6548	-13.7	7.2
1992.6576	-28.3	7.8
1993.6405	9.1	7.3
1995.5817	-2.9	9.8
1995.5844	1.7	8.7
1995.5871	2.7	8.0
1996.6166	1.8	7.5
1996.6193	-0.2	7.4
1996.6221	4.9	7.4
1996.6248	3.0	7.9
1998.6344	5.5	7.9
1998.6372	6.9	7.1
2003.5627	-7.9	8.0
2003.5654	-1.2	7.0
2003.5681	0.2	7.1
2003.5709	-1.9	7.1
2003.5736	-3.3	7.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> reub.ps
1992.6521	-11.0	8.0
1992.6548	-13.7	7.2
1992.6576	-28.3	7.8
1993.6405	9.1	7.3
1995.5817	-2.9	9.8
1995.5844	1.7	8.7
1995.5871	2.7	8.0
1996.6166	1.8	7.5
1996.6193	-0.2	7.4
1996.6221	4.9	7.4
1996.6248	3.0	7.9
1998.6344	5.5	7.9
1998.6372	6.9	7.1
2003.5627	-7.9	8.0
2003.5654	-1.2	7.0
2003.5681	0.2	7.1
2003.5709	-1.9	7.1
2003.5736	-3.3	7.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> reub.ps
SBreak	Segment0	Comp1
 1992.6521	   1.7609	Segment0	Comp1
 2003.5736	  -1.1290	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> reub.ps
1992.6521	-30.5	22.9
1992.6548	-35.7	21.9
1992.6576	-7.7	23.5
1993.6405	5.6	21.5
1995.5817	28.0	23.5
1995.5844	19.5	26.4
1995.5871	18.4	23.6
1996.6166	26.6	22.3
1996.6193	17.5	22.8
1996.6221	14.2	22.6
1996.6248	14.4	24.1
1998.6344	-19.9	24.1
1998.6372	-17.8	21.3
2003.5627	-4.3	23.1
2003.5654	-14.8	20.7
2003.5681	-28.5	21.0
2003.5709	-9.7	21.3
2003.5736	-12.4	20.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> reub.ps
1992.6521	-30.5	22.9
1992.6548	-35.7	21.9
1992.6576	-7.7	23.5
1993.6405	5.6	21.5
1995.5817	28.0	23.5
1995.5844	19.5	26.4
1995.5871	18.4	23.6
1996.6166	26.6	22.3
1996.6193	17.5	22.8
1996.6221	14.2	22.6
1996.6248	14.4	24.1
1998.6344	-19.9	24.1
1998.6372	-17.8	21.3
2003.5627	-4.3	23.1
2003.5654	-14.8	20.7
2003.5681	-28.5	21.0
2003.5709	-9.7	21.3
2003.5736	-12.4	20.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> reub.ps
SBreak	Segment0	Comp2
 1992.6521	  12.5589	Segment0	Comp2
 2003.5736	 -11.2813	Segment0	Comp2
EOT