#!/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	-4.8	2.9
2005.8242	-3.2	2.9
2005.8269	-8.2	3.5
2006.6291	-5.0	2.9
2006.6319	-8.6	2.9
2006.6346	-7.7	3.5
2007.6230	3.2	3.0
2007.6257	4.6	3.0
2007.6285	4.7	4.5
2008.5320	10.7	3.0
2008.5347	8.8	3.1
2008.5374	12.1	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	-4.8	2.9
2005.8242	-3.2	2.9
2005.8269	-8.2	3.5
2006.6291	-5.0	2.9
2006.6319	-8.6	2.9
2006.6346	-7.7	3.5
2007.6230	3.2	3.0
2007.6257	4.6	3.0
2007.6285	4.7	4.5
2008.5320	10.7	3.0
2008.5347	8.8	3.1
2008.5374	12.1	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.9063	Segment0	Comp0
 2008.5374	   9.1574	Segment0	Comp0
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 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.1	3.1
2005.8242	13.2	3.0
2005.8269	12.4	3.8
2006.6291	5.9	3.1
2006.6319	5.6	3.1
2006.6346	8.3	3.9
2007.6230	-5.5	3.2
2007.6257	-5.0	3.2
2007.6285	-5.1	5.2
2008.5320	-16.5	3.2
2008.5347	-22.7	3.5
2008.5374	-22.3	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.1	3.1
2005.8242	13.2	3.0
2005.8269	12.4	3.8
2006.6291	5.9	3.1
2006.6319	5.6	3.1
2006.6346	8.3	3.9
2007.6230	-5.5	3.2
2007.6257	-5.0	3.2
2007.6285	-5.1	5.2
2008.5320	-16.5	3.2
2008.5347	-22.7	3.5
2008.5374	-22.3	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.0410	Segment0	Comp1
 2008.5374	 -18.2548	Segment0	Comp1
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 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.4	10.4
2005.8242	7.2	9.9
2005.8269	16.9	12.1
2006.6291	-15.0	10.4
2006.6319	-6.0	10.4
2006.6346	-6.5	13.5
2007.6230	-2.1	10.3
2007.6257	1.5	10.5
2007.6285	-0.6	16.9
2008.5320	0.4	10.5
2008.5347	6.4	10.8
2008.5374	-13.9	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.4	10.4
2005.8242	7.2	9.9
2005.8269	16.9	12.1
2006.6291	-15.0	10.4
2006.6319	-6.0	10.4
2006.6346	-6.5	13.5
2007.6230	-2.1	10.3
2007.6257	1.5	10.5
2007.6285	-0.6	16.9
2008.5320	0.4	10.5
2008.5347	6.4	10.8
2008.5374	-13.9	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.5644	Segment0	Comp2
 2008.5374	  -2.9859	Segment0	Comp2
EOT