#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"AKHB 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	5.4	8.4
2001.6160	16.9	14.4
2002.4784	7.9	20.6
2002.4812	11.3	15.4
2003.4887	-2.5	21.8
2003.4915	-5.2	17.9
2003.4942	-3.2	9.2
2004.4799	-7.2	6.7
2004.4826	-9.7	24.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> akhb.ps
2004.4799	-7.2	6.7
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	5.4	8.4
2001.6160	16.9	14.4
2002.4784	7.9	20.6
2002.4812	11.3	15.4
2003.4887	-2.5	21.8
2003.4915	-5.2	17.9
2003.4942	-3.2	9.2
2004.4826	-9.7	24.4
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"AKHB 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	0.3	7.9
2001.6160	-1.9	17.0
2002.4784	-0.9	23.2
2002.4812	-6.1	17.6
2003.4887	-1.1	21.0
2003.4915	0.2	16.5
2003.4942	2.1	10.9
2004.4799	0.7	7.1
2004.4826	6.2	23.2
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> akhb.ps
2004.4799	0.7	7.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	0.3	7.9
2001.6160	-1.9	17.0
2002.4784	-0.9	23.2
2002.4812	-6.1	17.6
2003.4887	-1.1	21.0
2003.4915	0.2	16.5
2003.4942	2.1	10.9
2004.4826	6.2	23.2
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-30/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"AKHB 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	17.1	23.6
2001.6160	-29.3	40.2
2002.4784	16.3	59.9
2002.4812	-0.6	42.0
2003.4887	-18.2	65.8
2003.4915	-24.5	44.0
2003.4942	-7.8	29.1
2004.4799	20.9	22.5
2004.4826	20.9	65.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> akhb.ps
2004.4799	20.9	22.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akhb.ps
2001.6133	17.1	23.6
2001.6160	-29.3	40.2
2002.4784	16.3	59.9
2002.4812	-0.6	42.0
2003.4887	-18.2	65.8
2003.4915	-24.5	44.0
2003.4942	-7.8	29.1
2004.4826	20.9	65.5
EOT


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