#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2005/2009/-10/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	-6.1	2.9
2005.8242	-4.2	2.9
2005.8269	-7.5	3.5
2006.6291	-4.1	2.9
2006.6319	-8.5	2.9
2006.6346	-7.9	3.5
2007.6230	6.3	3.0
2007.6257	5.7	3.0
2007.6285	4.7	4.5
2008.5320	10.7	3.0
2008.5347	6.5	3.1
2008.5374	10.0	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	-6.1	2.9
2005.8242	-4.2	2.9
2005.8269	-7.5	3.5
2006.6291	-4.1	2.9
2006.6319	-8.5	2.9
2006.6346	-7.9	3.5
2007.6230	6.3	3.0
2007.6257	5.7	3.0
2007.6285	4.7	4.5
2008.5320	10.7	3.0
2008.5347	6.5	3.1
2008.5374	10.0	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	  -7.8736	Segment0	Comp0
 2008.5374	   9.1195	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2005/2009/-20/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"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	15.9	3.1
2005.8242	15.9	3.0
2005.8269	16.2	3.8
2006.6291	4.2	3.1
2006.6319	3.8	3.1
2006.6346	4.8	3.9
2007.6230	-5.4	3.2
2007.6257	-6.0	3.2
2007.6285	-5.5	5.2
2008.5320	-19.9	3.2
2008.5347	-19.9	3.5
2008.5374	-19.6	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	15.9	3.1
2005.8242	15.9	3.0
2005.8269	16.2	3.8
2006.6291	4.2	3.1
2006.6319	3.8	3.1
2006.6346	4.8	3.9
2007.6230	-5.4	3.2
2007.6257	-6.0	3.2
2007.6285	-5.5	5.2
2008.5320	-19.9	3.2
2008.5347	-19.9	3.5
2008.5374	-19.6	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	  15.7298	Segment0	Comp1
 2008.5374	 -19.0907	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2005/2009/-10/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"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	0.8	10.4
2005.8242	1.4	9.9
2005.8269	4.3	12.1
2006.6291	-3.6	10.4
2006.6319	2.6	10.4
2006.6346	-1.5	13.5
2007.6230	-7.7	10.3
2007.6257	-4.5	10.5
2007.6285	-6.4	16.9
2008.5320	-2.2	10.5
2008.5347	12.1	10.8
2008.5374	4.2	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	0.8	10.4
2005.8242	1.4	9.9
2005.8269	4.3	12.1
2006.6291	-3.6	10.4
2006.6319	2.6	10.4
2006.6346	-1.5	13.5
2007.6230	-7.7	10.3
2007.6257	-4.5	10.5
2007.6285	-6.4	16.9
2008.5320	-2.2	10.5
2008.5347	12.1	10.8
2008.5374	4.2	15.3
EOT


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