#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2003/2010/-40/30 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"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	-30.1	5.7
2003.1138	-30.6	5.7
2006.1446	-2.8	5.2
2006.1474	1.3	5.1
2009.2055	23.8	5.1
2009.2082	26.6	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	-30.1	5.7
2003.1138	-30.6	5.7
2006.1446	-2.8	5.2
2006.1474	1.3	5.1
2009.2055	23.8	5.1
2009.2082	26.6	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	 -29.5877	Segment0	Comp0
 2009.2082	  25.8154	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2003/2010/-30/60 \
-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	42.9	9.5
2003.1138	52.3	9.4
2006.1446	11.0	5.9
2006.1474	5.1	5.9
2009.2055	-25.4	5.9
2009.2082	-26.1	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	42.9	9.5
2003.1138	52.3	9.4
2006.1446	11.0	5.9
2006.1474	5.1	5.9
2009.2055	-25.4	5.9
2009.2082	-26.1	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	  45.6255	Segment0	Comp1
 2009.2082	 -26.5864	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.0	22.6
2003.1138	38.0	22.6
2006.1446	6.3	17.8
2006.1474	-9.4	17.4
2009.2055	-15.7	16.4
2009.2082	-22.0	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.0	22.6
2003.1138	38.0	22.6
2006.1446	6.3	17.8
2006.1474	-9.4	17.4
2009.2055	-15.7	16.4
2009.2082	-22.0	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.4820	Segment0	Comp2
 2009.2082	 -21.7533	Segment0	Comp2
EOT