#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> oakr.ps
1992.6713	-35.4	7.1
1996.6111	-20.2	6.6
1996.6138	-17.1	7.2
1998.6454	-9.7	6.8
1998.6481	-20.7	7.3
1998.6673	-13.1	7.0
1998.6700	-15.7	6.8
1998.6728	-0.1	6.5
2003.5681	11.7	6.3
2003.5709	13.5	6.4
2003.5736	16.1	6.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> oakr.ps
1992.6713	-35.4	7.1
1996.6111	-20.2	6.6
1996.6138	-17.1	7.2
1998.6454	-9.7	6.8
1998.6481	-20.7	7.3
1998.6673	-13.1	7.0
1998.6700	-15.7	6.8
1998.6728	-0.1	6.5
2003.5681	11.7	6.3
2003.5709	13.5	6.4
2003.5736	16.1	6.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> oakr.ps
SBreak	Segment0	Comp0
 1992.6713	 -38.5023	Segment0	Comp0
 2003.5736	  13.2326	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> oakr.ps
1992.6713	-19.8	8.0
1996.6111	6.1	7.2
1996.6138	-5.5	7.8
1998.6454	-9.7	7.3
1998.6481	3.7	9.8
1998.6673	10.9	8.7
1998.6700	-0.6	8.0
1998.6728	1.6	7.5
2003.5681	2.2	7.4
2003.5709	1.3	7.4
2003.5736	-1.1	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> oakr.ps
1992.6713	-19.8	8.0
1996.6111	6.1	7.2
1996.6138	-5.5	7.8
1998.6454	-9.7	7.3
1998.6481	3.7	9.8
1998.6673	10.9	8.7
1998.6700	-0.6	8.0
1998.6728	1.6	7.5
2003.5681	2.2	7.4
2003.5709	1.3	7.4
2003.5736	-1.1	7.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> oakr.ps
SBreak	Segment0	Comp1
 1992.6713	  -5.1816	Segment0	Comp1
 2003.5736	   1.5614	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> oakr.ps
1992.6713	-14.7	22.9
1996.6111	14.9	21.9
1996.6138	16.5	23.5
1998.6454	-19.2	21.5
1998.6481	2.6	23.5
1998.6673	-73.4	26.4
1998.6700	-110.7	23.6
1998.6728	3.9	22.3
2003.5681	16.3	22.8
2003.5709	19.4	22.6
2003.5736	16.0	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> oakr.ps
1992.6713	-14.7	22.9
1996.6111	14.9	21.9
1996.6138	16.5	23.5
1998.6454	-19.2	21.5
1998.6481	2.6	23.5
1998.6673	-73.4	26.4
1998.6700	-110.7	23.6
1998.6728	3.9	22.3
2003.5681	16.3	22.8
2003.5709	19.4	22.6
2003.5736	16.0	24.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> oakr.ps
SBreak	Segment0	Comp2
 1992.6713	 -34.7658	Segment0	Comp2
 2003.5736	  12.0698	Segment0	Comp2
EOT