#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993.5/2001.5/-40/40 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"North (mm)":WeSn  \
> cntr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cntr.ps
1993.6323	-22.1	7.1
1993.6350	-31.8	6.6
1995.5926	-16.3	7.2
1995.5981	-11.3	6.8
1997.6214	-1.1	7.3
1997.6269	-6.9	7.0
1999.6474	8.0	6.8
1999.6529	10.1	6.5
2001.1726	35.0	6.3
2001.1753	25.9	6.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cntr.ps
1993.6323	-22.1	7.1
1993.6350	-31.8	6.6
1995.5926	-16.3	7.2
1995.5981	-11.3	6.8
1997.6214	-1.1	7.3
1997.6269	-6.9	7.0
1999.6474	8.0	6.8
1999.6529	10.1	6.5
2001.1726	35.0	6.3
2001.1753	25.9	6.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cntr.ps
SBreak	Segment0	Comp0
 1993.6323	 -28.9076	Segment0	Comp0
 2001.1753	  25.3569	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001.5/-50/30 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"East (mm)":WeSn  \
-O >> cntr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cntr.ps
1993.6323	-42.2	8.0
1993.6350	-38.6	7.2
1995.5926	-18.4	7.8
1995.5981	-20.7	7.3
1997.6214	-4.5	9.8
1997.6269	2.0	8.7
1999.6474	12.9	8.0
1999.6529	9.3	7.5
2001.1726	18.7	7.4
2001.1753	22.3	7.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cntr.ps
1993.6323	-42.2	8.0
1993.6350	-38.6	7.2
1995.5926	-18.4	7.8
1995.5981	-20.7	7.3
1997.6214	-4.5	9.8
1997.6269	2.0	8.7
1999.6474	12.9	8.0
1999.6529	9.3	7.5
2001.1726	18.7	7.4
2001.1753	22.3	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cntr.ps
SBreak	Segment0	Comp1
 1993.6323	 -36.0435	Segment0	Comp1
 2001.1753	  22.7335	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001.5/-30/40 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"Up (mm)":WeSn  \
-O >> cntr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cntr.ps
1993.6323	-19.1	22.9
1993.6350	-19.7	21.9
1995.5926	22.4	23.5
1995.5981	-18.3	21.5
1997.6214	-5.5	23.5
1997.6269	-29.0	26.4
1999.6474	10.5	23.6
1999.6529	21.2	22.3
2001.1726	10.6	22.8
2001.1753	36.4	22.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cntr.ps
1993.6323	-19.1	22.9
1993.6350	-19.7	21.9
1995.5926	22.4	23.5
1995.5981	-18.3	21.5
1997.6214	-5.5	23.5
1997.6269	-29.0	26.4
1999.6474	10.5	23.6
1999.6529	21.2	22.3
2001.1726	10.6	22.8
2001.1753	36.4	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> cntr.ps
SBreak	Segment0	Comp2
 1993.6323	 -19.4060	Segment0	Comp2
 2001.1753	  19.2548	Segment0	Comp2
EOT