#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pabh.ps
1999.6392	-23.2	7.1
1999.6420	-19.7	6.6
1999.6447	-23.8	7.2
1999.6474	-19.1	6.8
1999.6502	-26.5	7.3
1999.6529	-21.6	7.0
1999.6557	-22.8	6.8
1999.6584	-23.5	6.5
1999.6611	-23.0	6.3
2001.1726	-2.9	6.4
2001.1753	-4.8	6.6
2001.1781	-0.9	6.7
2001.1808	-0.3	6.1
2001.1835	-1.5	6.7
2001.1863	-1.1	6.1
2001.1890	-3.1	6.1
2001.4736	1.8	6.5
2003.6338	26.1	6.5
2003.6366	26.5	4.8
2003.6393	28.8	6.9
2003.6421	26.0	7.3
2003.6448	28.2	7.0
2003.6475	24.4	6.8
2003.6503	25.1	6.5
2003.6530	25.7	6.7
2003.6557	27.0	6.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pabh.ps
1999.6392	-23.2	7.1
1999.6420	-19.7	6.6
1999.6447	-23.8	7.2
1999.6474	-19.1	6.8
1999.6502	-26.5	7.3
1999.6529	-21.6	7.0
1999.6557	-22.8	6.8
1999.6584	-23.5	6.5
1999.6611	-23.0	6.3
2001.1726	-2.9	6.4
2001.1753	-4.8	6.6
2001.1781	-0.9	6.7
2001.1808	-0.3	6.1
2001.1835	-1.5	6.7
2001.1863	-1.1	6.1
2001.1890	-3.1	6.1
2001.4736	1.8	6.5
2003.6338	26.1	6.5
2003.6366	26.5	4.8
2003.6393	28.8	6.9
2003.6421	26.0	7.3
2003.6448	28.2	7.0
2003.6475	24.4	6.8
2003.6503	25.1	6.5
2003.6530	25.7	6.7
2003.6557	27.0	6.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> pabh.ps
SBreak	Segment0	Comp0
 1999.6392	 -22.0841	Segment0	Comp0
 2003.6557	  27.0036	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pabh.ps
1999.6392	-16.4	8.0
1999.6420	-19.7	7.2
1999.6447	-27.2	7.8
1999.6474	-24.8	7.3
1999.6502	-29.6	9.8
1999.6529	-25.0	8.7
1999.6557	-21.8	8.0
1999.6584	-22.3	7.5
1999.6611	-25.6	7.4
2001.1726	-10.1	7.4
2001.1753	-3.2	7.9
2001.1781	-3.9	7.9
2001.1808	-6.0	7.1
2001.1835	-11.0	8.0
2001.1863	-5.4	7.0
2001.1890	-11.1	7.1
2001.4736	0.0	7.1
2003.6338	32.7	7.3
2003.6366	31.2	5.5
2003.6393	31.3	8.9
2003.6421	31.4	9.9
2003.6448	32.4	9.3
2003.6475	31.6	7.3
2003.6503	31.2	7.0
2003.6530	30.9	7.3
2003.6557	29.2	9.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pabh.ps
1999.6392	-16.4	8.0
1999.6420	-19.7	7.2
1999.6447	-27.2	7.8
1999.6474	-24.8	7.3
1999.6502	-29.6	9.8
1999.6529	-25.0	8.7
1999.6557	-21.8	8.0
1999.6584	-22.3	7.5
1999.6611	-25.6	7.4
2001.1726	-10.1	7.4
2001.1753	-3.2	7.9
2001.1781	-3.9	7.9
2001.1808	-6.0	7.1
2001.1835	-11.0	8.0
2001.1863	-5.4	7.0
2001.1890	-11.1	7.1
2001.4736	0.0	7.1
2003.6338	32.7	7.3
2003.6366	31.2	5.5
2003.6393	31.3	8.9
2003.6421	31.4	9.9
2003.6448	32.4	9.3
2003.6475	31.6	7.3
2003.6503	31.2	7.0
2003.6530	30.9	7.3
2003.6557	29.2	9.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> pabh.ps
SBreak	Segment0	Comp1
 1999.6392	 -25.1231	Segment0	Comp1
 2003.6557	  30.5487	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pabh.ps
1999.6392	6.0	22.9
1999.6420	-3.1	21.9
1999.6447	8.9	23.5
1999.6474	10.5	21.5
1999.6502	6.0	23.5
1999.6529	7.5	26.4
1999.6557	-1.2	23.6
1999.6584	4.3	22.3
1999.6611	-5.2	22.8
2001.1726	-6.6	22.6
2001.1753	0.4	24.1
2001.1781	-9.5	24.1
2001.1808	-12.6	21.3
2001.1835	-26.1	23.1
2001.1863	-7.3	20.7
2001.1890	-11.1	21.0
2001.4736	-15.1	21.3
2003.6338	3.0	20.8
2003.6366	11.4	15.9
2003.6393	12.0	27.3
2003.6421	2.8	32.2
2003.6448	10.1	27.6
2003.6475	4.7	21.9
2003.6503	2.6	21.1
2003.6530	10.1	21.7
2003.6557	-1.8	21.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pabh.ps
1999.6392	6.0	22.9
1999.6420	-3.1	21.9
1999.6447	8.9	23.5
1999.6474	10.5	21.5
1999.6502	6.0	23.5
1999.6529	7.5	26.4
1999.6557	-1.2	23.6
1999.6584	4.3	22.3
1999.6611	-5.2	22.8
2001.1726	-6.6	22.6
2001.1753	0.4	24.1
2001.1781	-9.5	24.1
2001.1808	-12.6	21.3
2001.1835	-26.1	23.1
2001.1863	-7.3	20.7
2001.1890	-11.1	21.0
2001.4736	-15.1	21.3
2003.6338	3.0	20.8
2003.6366	11.4	15.9
2003.6393	12.0	27.3
2003.6421	2.8	32.2
2003.6448	10.1	27.6
2003.6475	4.7	21.9
2003.6503	2.6	21.1
2003.6530	10.1	21.7
2003.6557	-1.8	21.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> pabh.ps
SBreak	Segment0	Comp2
 1999.6392	  -1.7042	Segment0	Comp2
 2003.6557	   2.1266	Segment0	Comp2
EOT