#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> boli.ps
1992.6548	-52.9	7.1
1993.6405	-36.5	6.6
1995.5844	-11.1	7.2
1995.5871	-28.2	6.8
1996.6166	-10.6	7.3
1996.6193	-7.2	7.0
1998.6344	-0.2	6.8
1998.6372	-5.0	6.5
2003.5654	46.3	6.3
2003.5681	41.1	6.4
2003.5709	46.2	6.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> boli.ps
1992.6548	-52.9	7.1
1993.6405	-36.5	6.6
1995.5844	-11.1	7.2
1995.5871	-28.2	6.8
1996.6166	-10.6	7.3
1996.6193	-7.2	7.0
1998.6344	-0.2	6.8
1998.6372	-5.0	6.5
2003.5654	46.3	6.3
2003.5681	41.1	6.4
2003.5709	46.2	6.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> boli.ps
SBreak	Segment0	Comp0
 1992.6548	 -46.6239	Segment0	Comp0
 2003.5709	  43.9598	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> boli.ps
1992.6548	-18.8	8.0
1993.6405	-2.5	7.2
1995.5844	-0.0	7.8
1995.5871	2.1	7.3
1996.6166	-1.3	9.8
1996.6193	1.6	8.7
1998.6344	3.2	8.0
1998.6372	4.1	7.5
2003.5654	-4.0	7.4
2003.5681	0.1	7.4
2003.5709	-0.4	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> boli.ps
1992.6548	-18.8	8.0
1993.6405	-2.5	7.2
1995.5844	-0.0	7.8
1995.5871	2.1	7.3
1996.6166	-1.3	9.8
1996.6193	1.6	8.7
1998.6344	3.2	8.0
1998.6372	4.1	7.5
2003.5654	-4.0	7.4
2003.5681	0.1	7.4
2003.5709	-0.4	7.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> boli.ps
SBreak	Segment0	Comp1
 1992.6548	  -2.9073	Segment0	Comp1
 2003.5709	   3.4008	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> boli.ps
1992.6548	-24.5	22.9
1993.6405	7.0	21.9
1995.5844	-2.2	23.5
1995.5871	9.3	21.5
1996.6166	30.1	23.5
1996.6193	17.1	26.4
1998.6344	-35.9	23.6
1998.6372	-4.0	22.3
2003.5654	-9.5	22.8
2003.5681	16.5	22.6
2003.5709	19.5	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> boli.ps
1992.6548	-24.5	22.9
1993.6405	7.0	21.9
1995.5844	-2.2	23.5
1995.5871	9.3	21.5
1996.6166	30.1	23.5
1996.6193	17.1	26.4
1998.6344	-35.9	23.6
1998.6372	-4.0	22.3
2003.5654	-9.5	22.8
2003.5681	16.5	22.6
2003.5709	19.5	24.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> boli.ps
SBreak	Segment0	Comp2
 1992.6548	   2.4087	Segment0	Comp2
 2003.5709	  -3.0913	Segment0	Comp2
EOT