#!/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  \
> clat.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> clat.ps
1993.6076	-29.1	7.1
1993.6104	-24.6	6.6
1994.6261	-3.1	7.2
1998.6508	6.7	6.8
1998.6536	11.8	7.3
1998.6563	1.0	7.0
2000.5537	29.3	6.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> clat.ps
1993.6076	-29.1	7.1
1993.6104	-24.6	6.6
1994.6261	-3.1	7.2
1998.6508	6.7	6.8
1998.6536	11.8	7.3
1998.6563	1.0	7.0
2000.5537	29.3	6.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> clat.ps
SBreak	Segment0	Comp0
 1993.6076	 -23.0549	Segment0	Comp0
 2000.5537	  23.4693	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> clat.ps
1993.6076	-12.5	8.0
1993.6104	-21.1	7.2
1994.6261	-17.6	7.8
1998.6508	12.7	7.3
1998.6536	-3.8	9.8
1998.6563	5.5	8.7
2000.5537	13.3	8.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> clat.ps
1993.6076	-12.5	8.0
1993.6104	-21.1	7.2
1994.6261	-17.6	7.8
1998.6508	12.7	7.3
1998.6536	-3.8	9.8
1998.6563	5.5	8.7
2000.5537	13.3	8.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> clat.ps
SBreak	Segment0	Comp1
 1993.6076	 -18.0792	Segment0	Comp1
 2000.5537	  13.1012	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> clat.ps
1993.6076	11.3	22.9
1993.6104	-9.6	21.9
1994.6261	-15.2	23.5
1998.6508	8.8	21.5
1998.6536	-8.1	23.5
1998.6563	15.5	26.4
2000.5537	-6.9	23.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> clat.ps
1993.6076	11.3	22.9
1993.6104	-9.6	21.9
1994.6261	-15.2	23.5
1998.6508	8.8	21.5
1998.6536	-8.1	23.5
1998.6563	15.5	26.4
2000.5537	-6.9	23.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> clat.ps
SBreak	Segment0	Comp2
 1993.6076	  -0.7531	Segment0	Comp2
 2000.5537	   0.7679	Segment0	Comp2
EOT