#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gray.ps
1993.6323	-52.4	7.1
1993.6350	-51.6	6.6
1993.6514	-48.5	7.2
1995.5953	-23.1	6.8
1995.5981	-22.3	7.3
1997.6159	-9.0	7.0
1997.6214	-8.8	6.8
1999.6392	4.6	6.5
1999.6420	9.7	6.3
2001.1726	34.0	6.4
2001.1753	31.9	6.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gray.ps
1993.6323	-52.4	7.1
1993.6350	-51.6	6.6
1993.6514	-48.5	7.2
1995.5953	-23.1	6.8
1995.5981	-22.3	7.3
1997.6159	-9.0	7.0
1997.6214	-8.8	6.8
1999.6392	4.6	6.5
1999.6420	9.7	6.3
2001.1726	34.0	6.4
2001.1753	31.9	6.6
EOT


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


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gray.ps
1993.6323	-53.0	8.0
1993.6350	-53.0	7.2
1993.6514	-58.5	7.8
1995.5953	-33.1	7.3
1995.5981	-32.7	9.8
1997.6159	-1.5	8.7
1997.6214	-4.0	8.0
1999.6392	18.1	7.5
1999.6420	14.6	7.4
2001.1726	20.8	7.4
2001.1753	25.0	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gray.ps
1993.6323	-53.0	8.0
1993.6350	-53.0	7.2
1993.6514	-58.5	7.8
1995.5953	-33.1	7.3
1995.5981	-32.7	9.8
1997.6159	-1.5	8.7
1997.6214	-4.0	8.0
1999.6392	18.1	7.5
1999.6420	14.6	7.4
2001.1726	20.8	7.4
2001.1753	25.0	7.9
EOT


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


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> gray.ps
1993.6323	16.3	22.9
1993.6350	12.0	21.9
1993.6514	-12.4	23.5
1995.5953	44.7	21.5
1995.5981	-9.3	23.5
1997.6159	-14.5	26.4
1997.6214	-11.3	23.6
1999.6392	-4.7	22.3
1999.6420	-4.7	22.8
2001.1726	-0.4	22.6
2001.1753	-4.5	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> gray.ps
1993.6323	16.3	22.9
1993.6350	12.0	21.9
1993.6514	-12.4	23.5
1995.5953	44.7	21.5
1995.5981	-9.3	23.5
1997.6159	-14.5	26.4
1997.6214	-11.3	23.6
1999.6392	-4.7	22.3
1999.6420	-4.7	22.8
2001.1726	-0.4	22.6
2001.1753	-4.5	24.1
EOT


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