#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-120/140 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-111.0	3.9
1994.2976	-93.4	4.7
1995.7240	-63.0	5.9
1996.8303	-54.6	5.8
1997.8843	-41.4	5.9
1997.9227	-51.4	6.2
1998.9220	-26.7	6.3
1998.9384	-31.4	5.5
2001.0302	10.7	6.0
2001.0439	15.2	6.0
2003.0946	42.3	5.7
2003.0973	45.0	5.5
2006.1063	91.8	5.5
2006.1090	92.5	5.2
2009.1918	131.7	5.2
2009.2110	137.2	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	-111.0	3.9
1994.2976	-93.4	4.7
1995.7240	-63.0	5.9
1996.8303	-54.6	5.8
1997.8843	-41.4	5.9
1997.9227	-51.4	6.2
1998.9220	-26.7	6.3
1998.9384	-31.4	5.5
2001.0302	10.7	6.0
2001.0439	15.2	6.0
2003.0946	42.3	5.7
2003.0973	45.0	5.5
2006.1063	91.8	5.5
2006.1090	92.5	5.2
2009.1918	131.7	5.2
2009.2110	137.2	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	-113.3644	Segment0	Comp0
 2009.2110	 136.3117	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-100/110 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"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	106.4	9.3
1994.2976	79.1	6.3
1995.7240	81.7	7.5
1996.8303	52.5	7.6
1997.8843	44.4	7.4
1997.9227	43.8	8.1
1998.9220	23.3	10.9
1998.9384	42.1	7.2
2001.0302	17.6	7.3
2001.0439	6.3	6.9
2003.0946	-9.4	7.4
2003.0973	-11.6	7.3
2006.1063	-46.2	6.4
2006.1090	-51.0	5.9
2009.1918	-89.9	5.6
2009.2110	-99.8	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	106.4	9.3
1994.2976	79.1	6.3
1995.7240	81.7	7.5
1996.8303	52.5	7.6
1997.8843	44.4	7.4
1997.9227	43.8	8.1
1998.9220	23.3	10.9
1998.9384	42.1	7.2
2001.0302	17.6	7.3
2001.0439	6.3	6.9
2003.0946	-9.4	7.4
2003.0973	-11.6	7.3
2006.1063	-46.2	6.4
2006.1090	-51.0	5.9
2009.1918	-89.9	5.6
2009.2110	-99.8	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	 104.5543	Segment0	Comp1
 2009.2110	 -89.3037	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.2	21.6
1994.2976	33.0	19.0
1995.7240	26.9	20.3
1996.8303	27.6	22.1
1997.8843	6.7	21.7
1997.9227	-4.2	21.9
1998.9220	-2.3	22.5
1998.9384	-5.4	19.2
2001.0302	-29.9	21.5
2001.0439	-17.5	21.2
2003.0946	-8.6	21.7
2003.0973	3.1	21.6
2006.1063	6.1	19.2
2006.1090	-16.3	18.4
2009.1918	-25.6	17.7
2009.2110	-35.0	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.2	21.6
1994.2976	33.0	19.0
1995.7240	26.9	20.3
1996.8303	27.6	22.1
1997.8843	6.7	21.7
1997.9227	-4.2	21.9
1998.9220	-2.3	22.5
1998.9384	-5.4	19.2
2001.0302	-29.9	21.5
2001.0439	-17.5	21.2
2003.0946	-8.6	21.7
2003.0973	3.1	21.6
2006.1063	6.1	19.2
2006.1090	-16.3	18.4
2009.1918	-25.6	17.7
2009.2110	-35.0	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.0248	Segment0	Comp2
 2009.2110	 -31.7663	Segment0	Comp2
EOT