#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1998.5/2004/-10/10 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"North (mm)":WeSn  \
> e040.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> e040.ps
1998.6645	-6.5	7.1
1998.6673	-3.6	6.6
1998.7111	2.4	7.2
1998.7138	3.9	6.8
2003.8420	3.8	7.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> e040.ps
1998.6645	-6.5	7.1
1998.6673	-3.6	6.6
1998.7111	2.4	7.2
1998.7138	3.9	6.8
2003.8420	3.8	7.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> e040.ps
SBreak	Segment0	Comp0
 1998.6645	  -1.0473	Segment0	Comp0
 2003.8420	   3.8469	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> e040.ps
1998.6645	-11.0	8.0
1998.6673	3.2	7.2
1998.7111	-4.4	7.8
1998.7138	0.3	7.3
2003.8420	7.7	9.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> e040.ps
1998.6645	-11.0	8.0
1998.6673	3.2	7.2
1998.7111	-4.4	7.8
1998.7138	0.3	7.3
2003.8420	7.7	9.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> e040.ps
SBreak	Segment0	Comp1
 1998.6645	  -2.5320	Segment0	Comp1
 2003.8420	   7.7399	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> e040.ps
1998.6645	-3.8	22.9
1998.6673	2.8	21.9
1998.7111	11.5	23.5
1998.7138	11.1	21.5
2003.8420	-16.0	23.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> e040.ps
1998.6645	-3.8	22.9
1998.6673	2.8	21.9
1998.7111	11.5	23.5
1998.7138	11.1	21.5
2003.8420	-16.0	23.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> e040.ps
SBreak	Segment0	Comp2
 1998.6645	   5.5162	Segment0	Comp2
 2003.8420	 -15.9186	Segment0	Comp2
EOT