#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rust.ps
1992.6576	-30.1	7.1
1993.6432	-29.1	6.6
1995.5817	-26.4	7.2
1995.5844	-23.9	6.8
1999.6283	-4.7	7.3
1999.6310	-3.1	7.0
1999.6338	1.1	6.8
1999.6365	5.4	6.5
2003.5709	20.9	6.3
2003.5736	21.5	6.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rust.ps
1992.6576	-30.1	7.1
1993.6432	-29.1	6.6
1995.5817	-26.4	7.2
1995.5844	-23.9	6.8
1999.6283	-4.7	7.3
1999.6310	-3.1	7.0
1999.6338	1.1	6.8
1999.6365	5.4	6.5
2003.5709	20.9	6.3
2003.5736	21.5	6.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> rust.ps
SBreak	Segment0	Comp0
 1992.6576	 -36.9781	Segment0	Comp0
 2003.5736	  21.0084	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1992.5/2004/-10/10 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"East (mm)":WeSn  \
-O >> rust.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rust.ps
1992.6576	7.6	8.0
1993.6432	7.0	7.2
1995.5817	-1.1	7.8
1995.5844	6.7	7.3
1999.6283	-1.9	9.8
1999.6310	-1.6	8.7
1999.6338	3.4	8.0
1999.6365	-4.9	7.5
2003.5709	-1.2	7.4
2003.5736	-3.7	7.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rust.ps
1992.6576	7.6	8.0
1993.6432	7.0	7.2
1995.5817	-1.1	7.8
1995.5844	6.7	7.3
1999.6283	-1.9	9.8
1999.6310	-1.6	8.7
1999.6338	3.4	8.0
1999.6365	-4.9	7.5
2003.5709	-1.2	7.4
2003.5736	-3.7	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> rust.ps
SBreak	Segment0	Comp1
 1992.6576	   6.6569	Segment0	Comp1
 2003.5736	  -3.0809	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rust.ps
1992.6576	-3.1	22.9
1993.6432	8.3	21.9
1995.5817	21.6	23.5
1995.5844	22.3	21.5
1999.6283	-20.6	23.5
1999.6310	-15.4	26.4
1999.6338	-11.9	23.6
1999.6365	-18.5	22.3
2003.5709	7.8	22.8
2003.5736	10.0	22.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rust.ps
1992.6576	-3.1	22.9
1993.6432	8.3	21.9
1995.5817	21.6	23.5
1995.5844	22.3	21.5
1999.6283	-20.6	23.5
1999.6310	-15.4	26.4
1999.6338	-11.9	23.6
1999.6365	-18.5	22.3
2003.5709	7.8	22.8
2003.5736	10.0	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> rust.ps
SBreak	Segment0	Comp2
 1992.6576	   4.5684	Segment0	Comp2
 2003.5736	  -2.6559	Segment0	Comp2
EOT