#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1995.5/2003/-50/70 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"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	56.2	6.3
1995.6802	52.8	6.6
1995.6830	52.7	4.4
1995.6857	51.5	4.2
1995.6884	49.3	6.4
1995.6912	47.7	5.8
1995.6939	62.8	8.3
1995.7049	63.7	6.5
1999.7871	-9.5	5.8
1999.7898	-7.0	5.9
1999.7926	-7.7	6.7
1999.7953	-6.9	7.3
1999.7980	-13.0	6.2
1999.8008	-7.1	4.4
1999.8035	-7.5	5.3
2001.4353	-31.8	7.1
2001.4380	-30.1	4.5
2001.4408	-28.5	4.6
2001.6023	-35.7	4.9
2001.6051	-37.2	4.5
2001.6078	-34.7	4.3
2002.6318	-46.5	4.9
2002.6345	-46.5	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	56.2	6.3
1995.6802	52.8	6.6
1995.6830	52.7	4.4
1995.6857	51.5	4.2
1995.6884	49.3	6.4
1995.6912	47.7	5.8
1995.6939	62.8	8.3
1995.7049	63.7	6.5
1999.7871	-9.5	5.8
1999.7898	-7.0	5.9
1999.7926	-7.7	6.7
1999.7953	-6.9	7.3
1999.7980	-13.0	6.2
1999.8008	-7.1	4.4
1999.8035	-7.5	5.3
2001.4353	-31.8	7.1
2001.4380	-30.1	4.5
2001.4408	-28.5	4.6
2001.6023	-35.7	4.9
2001.6051	-37.2	4.5
2001.6078	-34.7	4.3
2002.6318	-46.5	4.9
2002.6345	-46.5	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	  53.8334	Segment0	Comp0
 2002.6345	 -49.2885	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1995.5/2003/-40/60 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"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	48.1	7.7
1995.6802	44.7	7.8
1995.6830	50.3	6.7
1995.6857	39.3	6.1
1995.6884	45.8	8.1
1995.6912	42.5	7.4
1995.6939	42.6	10.0
1995.7049	41.7	8.2
1999.7871	-21.6	7.7
1999.7898	-8.3	8.1
1999.7926	-8.0	8.1
1999.7953	-2.7	8.2
1999.7980	-10.9	7.4
1999.8008	-10.4	6.0
1999.8035	-9.9	7.5
2001.4353	-30.8	9.1
2001.4380	-30.4	5.6
2001.4408	-30.8	5.6
2001.6023	-30.5	6.0
2001.6051	-26.3	5.6
2001.6078	-32.9	5.3
2002.6318	-36.3	5.7
2002.6345	-38.8	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	48.1	7.7
1995.6802	44.7	7.8
1995.6830	50.3	6.7
1995.6857	39.3	6.1
1995.6884	45.8	8.1
1995.6912	42.5	7.4
1995.6939	42.6	10.0
1995.7049	41.7	8.2
1999.7871	-21.6	7.7
1999.7898	-8.3	8.1
1999.7926	-8.0	8.1
1999.7953	-2.7	8.2
1999.7980	-10.9	7.4
1999.8008	-10.4	6.0
1999.8035	-9.9	7.5
2001.4353	-30.8	9.1
2001.4380	-30.4	5.6
2001.4408	-30.8	5.6
2001.6023	-30.5	6.0
2001.6051	-26.3	5.6
2001.6078	-32.9	5.3
2002.6318	-36.3	5.7
2002.6345	-38.8	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	  44.0978	Segment0	Comp1
 2002.6345	 -43.7440	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	-3.8	18.1
1995.6802	-13.4	20.8
1995.6830	10.0	12.5
1995.6857	15.6	12.7
1995.6884	1.4	18.8
1995.6912	8.5	17.8
1995.6939	15.1	25.1
1995.7049	1.0	19.0
1999.7871	2.6	17.4
1999.7898	16.8	17.0
1999.7926	15.0	18.2
1999.7953	2.1	20.9
1999.7980	26.6	18.2
1999.8008	13.1	14.4
1999.8035	0.9	15.7
2001.4353	-17.1	21.8
2001.4380	-10.2	13.4
2001.4408	-15.4	14.0
2001.6023	-6.6	15.8
2001.6051	-8.2	13.6
2001.6078	-7.5	13.1
2002.6318	-7.6	14.9
2002.6345	-4.8	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	-3.8	18.1
1995.6802	-13.4	20.8
1995.6830	10.0	12.5
1995.6857	15.6	12.7
1995.6884	1.4	18.8
1995.6912	8.5	17.8
1995.6939	15.1	25.1
1995.7049	1.0	19.0
1999.7871	2.6	17.4
1999.7898	16.8	17.0
1999.7926	15.0	18.2
1999.7953	2.1	20.9
1999.7980	26.6	18.2
1999.8008	13.1	14.4
1999.8035	0.9	15.7
2001.4353	-17.1	21.8
2001.4380	-10.2	13.4
2001.4408	-15.4	14.0
2001.6023	-6.6	15.8
2001.6051	-8.2	13.6
2001.6078	-7.5	13.1
2002.6318	-7.6	14.9
2002.6345	-4.8	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.6095	Segment0	Comp2
 2002.6345	  -6.6471	Segment0	Comp2
EOT