#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1995.5/2003/-20/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"WHIT North (mm)":WeSn  \
> whit.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whit.ps
1995.6775	-5.2	6.3
1995.6802	-8.6	6.6
1995.6830	-8.8	4.4
1995.6857	-9.8	4.2
1995.6884	-12.2	6.4
1995.6912	-13.6	5.8
1995.6939	1.5	8.3
1995.7049	2.6	6.5
1999.7871	-1.4	5.8
1999.7898	1.2	5.9
1999.7926	0.5	6.7
1999.7953	1.3	7.3
1999.7980	-4.6	6.2
1999.8008	1.2	4.4
1999.8035	0.9	5.3
2001.4353	4.2	7.1
2001.4380	5.9	4.5
2001.4408	7.6	4.6
2001.6023	3.2	4.9
2001.6051	1.6	4.5
2001.6078	4.2	4.3
2002.6318	9.7	4.9
2002.6345	9.8	7.2
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> whit.ps
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whit.ps
1995.6775	-5.2	6.3
1995.6802	-8.6	6.6
1995.6830	-8.8	4.4
1995.6857	-9.8	4.2
1995.6884	-12.2	6.4
1995.6912	-13.6	5.8
1995.6939	1.5	8.3
1995.7049	2.6	6.5
1999.7871	-1.4	5.8
1999.7898	1.2	5.9
1999.7926	0.5	6.7
1999.7953	1.3	7.3
1999.7980	-4.6	6.2
1999.8008	1.2	4.4
1999.8035	0.9	5.3
2001.4353	4.2	7.1
2001.4380	5.9	4.5
2001.4408	7.6	4.6
2001.6023	3.2	4.9
2001.6051	1.6	4.5
2001.6078	4.2	4.3
2002.6318	9.7	4.9
2002.6345	9.8	7.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whit.ps
SBreak	Segment0	Comp0
 1995.6775	  -7.6727	Segment0	Comp0
 2002.6345	   7.0249	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1995.5/2003/-20/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"WHIT East (mm)":WeSn  \
-O >> whit.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whit.ps
1995.6775	4.8	7.7
1995.6802	1.5	7.8
1995.6830	7.0	6.7
1995.6857	-4.1	6.1
1995.6884	2.5	8.1
1995.6912	-0.7	7.4
1995.6939	-0.5	10.0
1995.7049	-1.3	8.2
1999.7871	-14.0	7.7
1999.7898	-0.6	8.1
1999.7926	-0.3	8.1
1999.7953	5.1	8.2
1999.7980	-3.1	7.4
1999.8008	-2.6	6.0
1999.8035	-2.0	7.5
2001.4353	-2.7	9.1
2001.4380	-2.3	5.6
2001.4408	-2.6	5.6
2001.6023	-0.2	6.0
2001.6051	4.0	5.6
2001.6078	-2.6	5.3
2002.6318	6.6	5.7
2002.6345	4.2	7.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> whit.ps
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whit.ps
1995.6775	4.8	7.7
1995.6802	1.5	7.8
1995.6830	7.0	6.7
1995.6857	-4.1	6.1
1995.6884	2.5	8.1
1995.6912	-0.7	7.4
1995.6939	-0.5	10.0
1995.7049	-1.3	8.2
1999.7871	-14.0	7.7
1999.7898	-0.6	8.1
1999.7926	-0.3	8.1
1999.7953	5.1	8.2
1999.7980	-3.1	7.4
1999.8008	-2.6	6.0
1999.8035	-2.0	7.5
2001.4353	-2.7	9.1
2001.4380	-2.3	5.6
2001.4408	-2.6	5.6
2001.6023	-0.2	6.0
2001.6051	4.0	5.6
2001.6078	-2.6	5.3
2002.6318	6.6	5.7
2002.6345	4.2	7.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whit.ps
SBreak	Segment0	Comp1
 1995.6775	   0.7391	Segment0	Comp1
 2002.6345	  -0.7331	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1995.5/2003/-20/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"WHIT Up (mm)":WeSn  \
-O >> whit.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whit.ps
1995.6775	-4.0	18.1
1995.6802	-13.6	20.8
1995.6830	9.8	12.5
1995.6857	15.4	12.7
1995.6884	1.2	18.8
1995.6912	8.4	17.8
1995.6939	14.9	25.1
1995.7049	0.8	19.0
1999.7871	2.7	17.4
1999.7898	16.9	17.0
1999.7926	15.1	18.2
1999.7953	2.1	20.9
1999.7980	26.6	18.2
1999.8008	13.2	14.4
1999.8035	0.9	15.7
2001.4353	-17.0	21.8
2001.4380	-10.1	13.4
2001.4408	-15.3	14.0
2001.6023	-6.5	15.8
2001.6051	-8.1	13.6
2001.6078	-7.4	13.1
2002.6318	-7.4	14.9
2002.6345	-4.6	19.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> whit.ps
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whit.ps
1995.6775	-4.0	18.1
1995.6802	-13.6	20.8
1995.6830	9.8	12.5
1995.6857	15.4	12.7
1995.6884	1.2	18.8
1995.6912	8.4	17.8
1995.6939	14.9	25.1
1995.7049	0.8	19.0
1999.7871	2.7	17.4
1999.7898	16.9	17.0
1999.7926	15.1	18.2
1999.7953	2.1	20.9
1999.7980	26.6	18.2
1999.8008	13.2	14.4
1999.8035	0.9	15.7
2001.4353	-17.0	21.8
2001.4380	-10.1	13.4
2001.4408	-15.3	14.0
2001.6023	-6.5	15.8
2001.6051	-8.1	13.6
2001.6078	-7.4	13.1
2002.6318	-7.4	14.9
2002.6345	-4.6	19.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> whit.ps
SBreak	Segment0	Comp2
 1995.6775	   6.4233	Segment0	Comp2
 2002.6345	  -6.4598	Segment0	Comp2
EOT