#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2006/-10/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"Y382 North (mm)":WeSn  \
> y382.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> y382.ps
2001.6709	3.3	7.1
2001.6873	-1.5	7.1
2001.6900	-3.5	4.9
2005.6901	-2.6	8.3
2005.6929	0.8	7.3
2005.6956	2.0	4.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> y382.ps
2001.6709	3.3	7.1
2001.6873	-1.5	7.1
2001.6900	-3.5	4.9
2005.6901	-2.6	8.3
2005.6929	0.8	7.3
2005.6956	2.0	4.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> y382.ps
SBreak	Segment0	Comp0
 2001.6709	  -0.7600	Segment0	Comp0
 2005.6956	   0.9331	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2006/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"Y382 East (mm)":WeSn  \
-O >> y382.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> y382.ps
2001.6709	-0.3	8.1
2001.6873	-5.0	8.4
2001.6900	-6.2	6.7
2005.6901	-0.9	9.8
2005.6929	4.6	8.6
2005.6956	3.8	5.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> y382.ps
2001.6709	-0.3	8.1
2001.6873	-5.0	8.4
2001.6900	-6.2	6.7
2005.6901	-0.9	9.8
2005.6929	4.6	8.6
2005.6956	3.8	5.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> y382.ps
SBreak	Segment0	Comp1
 2001.6709	  -3.8098	Segment0	Comp1
 2005.6956	   3.6038	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2006/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"Y382 Up (mm)":WeSn  \
-O >> y382.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> y382.ps
2001.6709	4.8	21.9
2001.6873	-2.6	20.8
2001.6900	-0.4	18.5
2005.6901	-8.1	27.8
2005.6929	0.5	24.0
2005.6956	1.5	15.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> y382.ps
2001.6709	4.8	21.9
2001.6873	-2.6	20.8
2001.6900	-0.4	18.5
2005.6901	-8.1	27.8
2005.6929	0.5	24.0
2005.6956	1.5	15.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> y382.ps
SBreak	Segment0	Comp2
 2001.6709	   0.2813	Segment0	Comp2
 2005.6956	  -0.3003	Segment0	Comp2
EOT