#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> stea.ps
1992.6603	-43.8	7.1
1992.6740	-31.3	6.6
1996.6166	-14.5	7.2
1996.6193	-20.5	6.8
1998.6454	-1.9	7.3
1998.6481	-16.6	7.0
1999.6283	-6.1	6.8
1999.6310	-4.3	6.5
1999.6338	-1.5	6.3
2003.5599	24.7	6.4
2003.5627	9.5	6.6
2003.5654	20.6	6.7
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> stea.ps
1992.6603	-43.8	7.1
1992.6740	-31.3	6.6
1996.6166	-14.5	7.2
1996.6193	-20.5	6.8
1998.6454	-1.9	7.3
1998.6481	-16.6	7.0
1999.6283	-6.1	6.8
1999.6310	-4.3	6.5
1999.6338	-1.5	6.3
2003.5599	24.7	6.4
2003.5627	9.5	6.6
2003.5654	20.6	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> stea.ps
SBreak	Segment0	Comp0
 1992.6603	 -39.4473	Segment0	Comp0
 2003.5654	  17.2809	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> stea.ps
1992.6603	-7.1	8.0
1992.6740	9.9	7.2
1996.6166	15.5	7.8
1996.6193	9.2	7.3
1998.6454	-2.4	9.8
1998.6481	10.7	8.7
1999.6283	-4.6	8.0
1999.6310	-1.7	7.5
1999.6338	4.7	7.4
2003.5599	3.3	7.4
2003.5627	-4.8	7.9
2003.5654	-8.6	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> stea.ps
1992.6603	-7.1	8.0
1992.6740	9.9	7.2
1996.6166	15.5	7.8
1996.6193	9.2	7.3
1998.6454	-2.4	9.8
1998.6481	10.7	8.7
1999.6283	-4.6	8.0
1999.6310	-1.7	7.5
1999.6338	4.7	7.4
2003.5599	3.3	7.4
2003.5627	-4.8	7.9
2003.5654	-8.6	7.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> stea.ps
SBreak	Segment0	Comp1
 1992.6603	   8.0287	Segment0	Comp1
 2003.5654	  -2.9516	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> stea.ps
1992.6603	-3.2	22.9
1992.6740	-51.6	21.9
1996.6166	15.5	23.5
1996.6193	5.6	21.5
1998.6454	-20.0	23.5
1998.6481	25.0	26.4
1999.6283	-0.9	23.6
1999.6310	-6.8	22.3
1999.6338	-7.5	22.8
2003.5599	19.7	22.6
2003.5627	10.4	24.1
2003.5654	-19.3	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> stea.ps
1992.6603	-3.2	22.9
1992.6740	-51.6	21.9
1996.6166	15.5	23.5
1996.6193	5.6	21.5
1998.6454	-20.0	23.5
1998.6481	25.0	26.4
1999.6283	-0.9	23.6
1999.6310	-6.8	22.3
1999.6338	-7.5	22.8
2003.5599	19.7	22.6
2003.5627	10.4	24.1
2003.5654	-19.3	24.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> stea.ps
SBreak	Segment0	Comp2
 1992.6603	 -15.4261	Segment0	Comp2
 2003.5654	   7.0613	Segment0	Comp2
EOT