#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993.5/2004/-80/50 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"North (mm)":WeSn  \
> grev.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> grev.ps
1993.6487	-79.7	7.1
1995.5981	-52.0	6.6
1997.6159	-33.2	7.2
1997.6187	-33.4	6.8
1999.6392	-12.3	7.3
1999.6420	-4.2	7.0
1999.6447	-9.5	6.8
2001.1808	12.8	6.5
2003.6338	47.1	6.3
2003.6366	46.7	6.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> grev.ps
1993.6487	-79.7	7.1
1995.5981	-52.0	6.6
1997.6159	-33.2	7.2
1997.6187	-33.4	6.8
1999.6392	-12.3	7.3
1999.6420	-4.2	7.0
1999.6447	-9.5	6.8
2001.1808	12.8	6.5
2003.6338	47.1	6.3
2003.6366	46.7	6.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> grev.ps
SBreak	Segment0	Comp0
 1993.6487	 -84.2198	Segment0	Comp0
 2003.6366	  45.7537	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> grev.ps
1993.6487	-102.7	8.0
1995.5981	-68.9	7.2
1997.6159	-30.7	7.8
1997.6187	-37.5	7.3
1999.6392	-3.8	9.8
1999.6420	-5.9	8.7
1999.6447	-13.5	8.0
2001.1808	10.0	7.5
2003.6338	47.6	7.4
2003.6366	46.7	7.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> grev.ps
1993.6487	-102.7	8.0
1995.5981	-68.9	7.2
1997.6159	-30.7	7.8
1997.6187	-37.5	7.3
1999.6392	-3.8	9.8
1999.6420	-5.9	8.7
1999.6447	-13.5	8.0
2001.1808	10.0	7.5
2003.6338	47.6	7.4
2003.6366	46.7	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> grev.ps
SBreak	Segment0	Comp1
 1993.6487	 -91.5246	Segment0	Comp1
 2003.6366	  46.8178	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> grev.ps
1993.6487	-37.9	22.9
1995.5981	-7.7	21.9
1997.6159	-1.3	23.5
1997.6187	-8.5	21.5
1999.6392	-10.1	23.5
1999.6420	5.4	26.4
1999.6447	6.9	23.6
2001.1808	-6.1	22.3
2003.6338	-0.2	22.8
2003.6366	8.8	22.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> grev.ps
1993.6487	-37.9	22.9
1995.5981	-7.7	21.9
1997.6159	-1.3	23.5
1997.6187	-8.5	21.5
1999.6392	-10.1	23.5
1999.6420	5.4	26.4
1999.6447	6.9	23.6
2001.1808	-6.1	22.3
2003.6338	-0.2	22.8
2003.6366	8.8	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> grev.ps
SBreak	Segment0	Comp2
 1993.6487	  -9.1278	Segment0	Comp2
 2003.6366	   4.9588	Segment0	Comp2
EOT