#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 2271.ps
2003.1110	13.8	5.7
2003.1138	13.2	5.7
2006.1446	0.2	5.2
2006.1474	4.2	5.1
2009.2055	-14.4	5.1
2009.2082	-11.7	5.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 2271.ps
2003.1110	13.8	5.7
2003.1138	13.2	5.7
2006.1446	0.2	5.2
2006.1474	4.2	5.1
2009.2055	-14.4	5.1
2009.2082	-11.7	5.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 2271.ps
SBreak	Segment0	Comp0
 2003.1110	  14.2966	Segment0	Comp0
 2009.2082	 -12.4739	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2003/2010/-60/100 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"2271 East (mm)":WeSn  \
-O >> 2271.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 2271.ps
2003.1110	88.6	9.5
2003.1138	98.1	9.4
2006.1446	20.7	5.9
2006.1474	14.9	5.9
2009.2055	-52.1	5.9
2009.2082	-52.8	5.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 2271.ps
2003.1110	88.6	9.5
2003.1138	98.1	9.4
2006.1446	20.7	5.9
2006.1474	14.9	5.9
2009.2055	-52.1	5.9
2009.2082	-52.8	5.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 2271.ps
SBreak	Segment0	Comp1
 2003.1110	  91.3819	Segment0	Comp1
 2009.2082	 -53.2490	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2003/2010/-30/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"2271 Up (mm)":WeSn  \
-O >> 2271.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 2271.ps
2003.1110	36.1	22.6
2003.1138	38.2	22.6
2006.1446	6.4	17.8
2006.1474	-9.4	17.4
2009.2055	-15.8	16.4
2009.2082	-22.1	17.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 2271.ps
2003.1110	36.1	22.6
2003.1138	38.2	22.6
2006.1446	6.4	17.8
2006.1474	-9.4	17.4
2009.2055	-15.8	16.4
2009.2082	-22.1	17.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> 2271.ps
SBreak	Segment0	Comp2
 2003.1110	  31.6229	Segment0	Comp2
 2009.2082	 -21.8507	Segment0	Comp2
EOT