#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mudd.ps
1993.6542	-13.5	7.1
1993.6569	-17.0	6.6
1995.5953	-13.2	7.2
1995.5981	-13.6	6.8
1997.6351	-8.3	7.3
1997.6378	-1.3	7.0
1999.6502	0.4	6.8
1999.6529	4.2	6.5
2001.1781	14.2	6.3
2001.1808	13.6	6.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> mudd.ps
1993.6542	-13.5	7.1
1993.6569	-17.0	6.6
1995.5953	-13.2	7.2
1995.5981	-13.6	6.8
1997.6351	-8.3	7.3
1997.6378	-1.3	7.0
1999.6502	0.4	6.8
1999.6529	4.2	6.5
2001.1781	14.2	6.3
2001.1808	13.6	6.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> mudd.ps
SBreak	Segment0	Comp0
 1993.6542	 -17.5392	Segment0	Comp0
 2001.1808	  11.9785	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mudd.ps
1993.6542	-10.7	8.0
1993.6569	-0.3	7.2
1995.5953	-10.7	7.8
1995.5981	0.6	7.3
1997.6351	4.5	9.8
1997.6378	7.6	8.7
1999.6502	-5.6	8.0
1999.6529	0.6	7.5
2001.1781	3.9	7.4
2001.1808	1.9	7.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> mudd.ps
1993.6542	-10.7	8.0
1993.6569	-0.3	7.2
1995.5953	-10.7	7.8
1995.5981	0.6	7.3
1997.6351	4.5	9.8
1997.6378	7.6	8.7
1999.6502	-5.6	8.0
1999.6529	0.6	7.5
2001.1781	3.9	7.4
2001.1808	1.9	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> mudd.ps
SBreak	Segment0	Comp1
 1993.6542	  -5.8568	Segment0	Comp1
 2001.1808	   3.4191	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mudd.ps
1993.6542	5.4	22.9
1993.6569	15.5	21.9
1995.5953	36.7	23.5
1995.5981	24.7	21.5
1997.6351	17.3	23.5
1997.6378	2.6	26.4
1999.6502	2.8	23.6
1999.6529	5.7	22.3
2001.1781	-22.2	22.8
2001.1808	-28.3	22.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> mudd.ps
1993.6542	5.4	22.9
1993.6569	15.5	21.9
1995.5953	36.7	23.5
1995.5981	24.7	21.5
1997.6351	17.3	23.5
1997.6378	2.6	26.4
1999.6502	2.8	23.6
1999.6529	5.7	22.3
2001.1781	-22.2	22.8
2001.1808	-28.3	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> mudd.ps
SBreak	Segment0	Comp2
 1993.6542	  22.2042	Segment0	Comp2
 2001.1808	 -15.6204	Segment0	Comp2
EOT