#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-20/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"HAGG North (mm)":WeSn  \
> hagg.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hagg.ps
1993.0355	-11.5	3.9
1994.2976	-10.9	4.7
1995.7240	0.2	5.9
1996.8303	-6.4	5.8
1997.8843	-7.5	5.9
1997.9227	-18.2	6.2
1998.9220	-6.9	6.3
1998.9384	-11.9	5.5
2001.0302	2.0	6.0
2001.0439	6.3	6.0
2003.0946	5.5	5.7
2003.0973	8.2	5.5
2006.1063	14.1	5.5
2006.1090	14.8	5.2
2009.1918	12.2	5.2
2009.2110	17.5	5.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hagg.ps
1993.0355	-11.5	3.9
1994.2976	-10.9	4.7
1995.7240	0.2	5.9
1996.8303	-6.4	5.8
1997.8843	-7.5	5.9
1997.9227	-18.2	6.2
1998.9220	-6.9	6.3
1998.9384	-11.9	5.5
2001.0302	2.0	6.0
2001.0439	6.3	6.0
2003.0946	5.5	5.7
2003.0973	8.2	5.5
2006.1063	14.1	5.5
2006.1090	14.8	5.2
2009.1918	12.2	5.2
2009.2110	17.5	5.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hagg.ps
SBreak	Segment0	Comp0
 1993.0355	 -13.7908	Segment0	Comp0
 2009.2110	  16.5823	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-190/210 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"HAGG East (mm)":WeSn  \
-O >> hagg.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hagg.ps
1993.0355	209.3	9.3
1994.2976	167.1	6.3
1995.7240	152.7	7.5
1996.8303	110.5	7.6
1997.8843	90.0	7.4
1997.9227	89.0	8.1
1998.9220	56.8	10.9
1998.9384	75.3	7.2
2001.0302	26.2	7.3
2001.0439	14.7	6.9
2003.0946	-25.1	7.4
2003.0973	-27.4	7.3
2006.1063	-97.4	6.4
2006.1090	-102.2	5.9
2009.1918	-177.5	5.6
2009.2110	-187.7	6.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hagg.ps
1993.0355	209.3	9.3
1994.2976	167.1	6.3
1995.7240	152.7	7.5
1996.8303	110.5	7.6
1997.8843	90.0	7.4
1997.9227	89.0	8.1
1998.9220	56.8	10.9
1998.9384	75.3	7.2
2001.0302	26.2	7.3
2001.0439	14.7	6.9
2003.0946	-25.1	7.4
2003.0973	-27.4	7.3
2006.1063	-97.4	6.4
2006.1090	-102.2	5.9
2009.1918	-177.5	5.6
2009.2110	-187.7	6.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hagg.ps
SBreak	Segment0	Comp1
 1993.0355	 207.3627	Segment0	Comp1
 2009.2110	-177.1162	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-40/70 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"HAGG Up (mm)":WeSn  \
-O >> hagg.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hagg.ps
1993.0355	66.6	21.6
1994.2976	33.3	19.0
1995.7240	27.2	20.3
1996.8303	27.8	22.1
1997.8843	6.8	21.7
1997.9227	-4.0	21.9
1998.9220	-2.2	22.5
1998.9384	-5.3	19.2
2001.0302	-29.9	21.5
2001.0439	-17.5	21.2
2003.0946	-8.6	21.7
2003.0973	3.0	21.6
2006.1063	5.9	19.2
2006.1090	-16.5	18.4
2009.1918	-26.0	17.7
2009.2110	-35.4	17.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hagg.ps
1993.0355	66.6	21.6
1994.2976	33.3	19.0
1995.7240	27.2	20.3
1996.8303	27.8	22.1
1997.8843	6.8	21.7
1997.9227	-4.0	21.9
1998.9220	-2.2	22.5
1998.9384	-5.3	19.2
2001.0302	-29.9	21.5
2001.0439	-17.5	21.2
2003.0946	-8.6	21.7
2003.0973	3.0	21.6
2006.1063	5.9	19.2
2006.1090	-16.5	18.4
2009.1918	-26.0	17.7
2009.2110	-35.4	17.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> hagg.ps
SBreak	Segment0	Comp2
 1993.0355	  33.4226	Segment0	Comp2
 2009.2110	 -32.1490	Segment0	Comp2
EOT