#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> asto.ps
1993.6268	-30.0	7.1
1993.6487	-20.5	6.6
1993.6514	-23.0	7.2
1994.6261	-14.6	6.8
1994.6289	-14.4	7.3
1998.6508	10.8	7.0
1998.6536	14.5	6.8
1998.6563	9.5	6.5
2000.5564	22.0	6.3
2000.5592	28.8	6.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> asto.ps
1993.6268	-30.0	7.1
1993.6487	-20.5	6.6
1993.6514	-23.0	7.2
1994.6261	-14.6	6.8
1994.6289	-14.4	7.3
1998.6508	10.8	7.0
1998.6536	14.5	6.8
1998.6563	9.5	6.5
2000.5564	22.0	6.3
2000.5592	28.8	6.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> asto.ps
SBreak	Segment0	Comp0
 1993.6268	 -23.4875	Segment0	Comp0
 2000.5592	  25.4409	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> asto.ps
1993.6268	-18.8	8.0
1993.6487	-32.6	7.2
1993.6514	-32.6	7.8
1994.6261	-27.2	7.3
1994.6289	-19.7	9.8
1998.6508	18.0	8.7
1998.6536	1.1	8.0
1998.6563	15.2	7.5
2000.5564	26.5	7.4
2000.5592	24.3	7.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> asto.ps
1993.6268	-18.8	8.0
1993.6487	-32.6	7.2
1993.6514	-32.6	7.8
1994.6261	-27.2	7.3
1994.6289	-19.7	9.8
1998.6508	18.0	8.7
1998.6536	1.1	8.0
1998.6563	15.2	7.5
2000.5564	26.5	7.4
2000.5592	24.3	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> asto.ps
SBreak	Segment0	Comp1
 1993.6268	 -29.8554	Segment0	Comp1
 2000.5592	  25.9991	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001/-50/50 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"Up (mm)":WeSn  \
-O >> asto.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> asto.ps
1993.6268	18.2	22.9
1993.6487	40.3	21.9
1993.6514	25.9	23.5
1994.6261	-11.4	21.5
1994.6289	-9.9	23.5
1998.6508	-41.3	26.4
1998.6536	-47.8	23.6
1998.6563	-17.8	22.3
2000.5564	4.2	22.8
2000.5592	10.9	22.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> asto.ps
1993.6268	18.2	22.9
1993.6487	40.3	21.9
1993.6514	25.9	23.5
1994.6261	-11.4	21.5
1994.6289	-9.9	23.5
1998.6508	-41.3	26.4
1998.6536	-47.8	23.6
1998.6563	-17.8	22.3
2000.5564	4.2	22.8
2000.5592	10.9	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> asto.ps
SBreak	Segment0	Comp2
 1993.6268	  11.1599	Segment0	Comp2
 2000.5592	 -12.4898	Segment0	Comp2
EOT