#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2003.5/-10/20 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"North (mm)":WeSn  \
> akhb.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akhb.ps
2001.6133	3.6	8.1
2001.6160	14.5	4.9
2002.4784	4.9	3.5
2002.4812	10.0	3.7
2003.4887	-7.4	3.7
2003.4915	-9.2	3.6
2003.4942	-7.2	3.7
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	3.6	8.1
2001.6160	14.5	4.9
2002.4784	4.9	3.5
2002.4812	10.0	3.7
2003.4887	-7.4	3.7
2003.4915	-9.2	3.6
2003.4942	-7.2	3.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akhb.ps
SBreak	Segment0	Comp0
 2001.6133	  15.1234	Segment0	Comp0
 2003.4942	  -7.2159	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akhb.ps
2001.6133	2.7	9.0
2001.6160	-1.3	5.1
2002.4784	0.9	3.9
2002.4812	-4.2	4.2
2003.4887	0.0	4.3
2003.4915	0.9	4.3
2003.4942	3.1	4.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	2.7	9.0
2001.6160	-1.3	5.1
2002.4784	0.9	3.9
2002.4812	-4.2	4.2
2003.4887	0.0	4.3
2003.4915	0.9	4.3
2003.4942	3.1	4.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akhb.ps
SBreak	Segment0	Comp1
 2001.6133	  -1.6694	Segment0	Comp1
 2003.4942	   0.9190	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akhb.ps
2001.6133	31.1	23.1
2001.6160	-15.5	15.1
2002.4784	24.1	11.2
2002.4812	6.8	12.0
2003.4887	-11.8	12.0
2003.4915	-18.4	11.3
2003.4942	-0.8	12.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	31.1	23.1
2001.6160	-15.5	15.1
2002.4784	24.1	11.2
2002.4812	6.8	12.0
2003.4887	-11.8	12.0
2003.4915	-18.4	11.3
2003.4942	-0.8	12.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> akhb.ps
SBreak	Segment0	Comp2
 2001.6133	  13.6971	Segment0	Comp2
 2003.4942	  -6.8987	Segment0	Comp2
EOT