#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2005/2009/-20/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"SKEL North (mm)":WeSn  \
> skel.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> skel.ps
2005.8214	11.6	2.9
2005.8242	13.1	2.9
2005.8269	8.2	3.5
2006.6264	-5.0	8.0
2006.6291	0.4	2.9
2006.6319	-3.3	2.9
2006.6346	-2.4	3.5
2007.6230	-5.1	3.0
2007.6257	-3.7	3.0
2007.6285	-3.7	4.5
2008.5347	-12.1	3.1
2008.5374	-8.7	4.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> skel.ps
2005.8214	11.6	2.9
2005.8242	13.1	2.9
2005.8269	8.2	3.5
2006.6264	-5.0	8.0
2006.6291	0.4	2.9
2006.6319	-3.3	2.9
2006.6346	-2.4	3.5
2007.6230	-5.1	3.0
2007.6257	-3.7	3.0
2007.6285	-3.7	4.5
2008.5347	-12.1	3.1
2008.5374	-8.7	4.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> skel.ps
SBreak	Segment0	Comp0
 2005.8214	   8.5129	Segment0	Comp0
 2008.5374	 -12.3621	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2005/2009/-50/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"SKEL East (mm)":WeSn  \
-O >> skel.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> skel.ps
2005.8214	26.1	3.1
2005.8242	24.1	3.0
2005.8269	23.3	3.8
2006.6264	7.6	9.1
2006.6291	7.3	3.1
2006.6319	6.8	3.1
2006.6346	9.5	3.9
2007.6230	-16.1	3.2
2007.6257	-15.6	3.2
2007.6285	-15.6	5.2
2008.5347	-44.1	3.5
2008.5374	-43.8	4.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> skel.ps
2005.8214	26.1	3.1
2005.8242	24.1	3.0
2005.8269	23.3	3.8
2006.6264	7.6	9.1
2006.6291	7.3	3.1
2006.6319	6.8	3.1
2006.6346	9.5	3.9
2007.6230	-16.1	3.2
2007.6257	-15.6	3.2
2007.6285	-15.6	5.2
2008.5347	-44.1	3.5
2008.5374	-43.8	4.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> skel.ps
SBreak	Segment0	Comp1
 2005.8214	  26.1913	Segment0	Comp1
 2008.5374	 -40.3976	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2005/2009/-30/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"SKEL Up (mm)":WeSn  \
-O >> skel.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> skel.ps
2005.8214	5.9	10.4
2005.8242	7.8	9.9
2005.8269	17.5	12.1
2006.6264	-29.7	26.3
2006.6291	-14.4	10.4
2006.6319	-5.4	10.4
2006.6346	-6.0	13.5
2007.6230	-1.6	10.3
2007.6257	2.1	10.5
2007.6285	-0.2	16.9
2008.5347	6.9	10.8
2008.5374	-13.3	15.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> skel.ps
2005.8214	5.9	10.4
2005.8242	7.8	9.9
2005.8269	17.5	12.1
2006.6264	-29.7	26.3
2006.6291	-14.4	10.4
2006.6319	-5.4	10.4
2006.6346	-6.0	13.5
2007.6230	-1.6	10.3
2007.6257	2.1	10.5
2007.6285	-0.2	16.9
2008.5347	6.9	10.8
2008.5374	-13.3	15.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> skel.ps
SBreak	Segment0	Comp2
 2005.8214	   2.7776	Segment0	Comp2
 2008.5374	  -4.1077	Segment0	Comp2
EOT