#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2006/-30/30 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"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	23.6	7.1
2001.6873	18.5	7.1
2001.6900	16.6	4.9
2005.6901	-27.4	8.3
2005.6929	-24.0	7.3
2005.6956	-22.8	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	23.6	7.1
2001.6873	18.5	7.1
2001.6900	16.6	4.9
2005.6901	-27.4	8.3
2005.6929	-24.0	7.3
2005.6956	-22.8	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	  19.4820	Segment0	Comp0
 2005.6956	 -23.9184	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2006/-30/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"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	30.3	8.1
2001.6873	25.4	8.4
2001.6900	24.1	6.7
2005.6901	-29.7	9.8
2005.6929	-24.4	8.6
2005.6956	-25.2	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	30.3	8.1
2001.6873	25.4	8.4
2001.6900	24.1	6.7
2005.6901	-29.7	9.8
2005.6929	-24.4	8.6
2005.6956	-25.2	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	  26.8297	Segment0	Comp1
 2005.6956	 -25.3791	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.9	21.9
2001.6873	-2.7	20.8
2001.6900	-0.3	18.5
2005.6901	-8.2	27.8
2005.6929	0.4	24.0
2005.6956	1.4	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.9	21.9
2001.6873	-2.7	20.8
2001.6900	-0.3	18.5
2005.6901	-8.2	27.8
2005.6929	0.4	24.0
2005.6956	1.4	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.3345	Segment0	Comp2
 2005.6956	  -0.3571	Segment0	Comp2
EOT