#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1999/2005.5/-10/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"DIAB North (mm)":WeSn  \
> diab.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> diab.ps
1999.0288	8.8	5.6
1999.0315	2.9	5.8
1999.0343	-6.4	5.8
1999.0370	3.8	5.9
1999.0507	9.8	5.9
1999.0534	7.0	6.7
1999.0562	4.5	6.6
1999.0589	9.7	6.8
2001.0795	4.8	5.2
2001.0822	1.4	3.3
2001.0850	7.6	6.1
2001.0877	8.3	5.2
2001.0987	3.0	5.5
2001.1014	4.9	5.7
2001.1041	4.8	5.5
2001.1069	7.8	5.8
2001.2410	4.0	6.7
2002.1555	6.4	3.0
2002.1582	4.7	5.9
2002.1610	1.2	5.7
2004.1679	-7.2	5.3
2004.1706	-5.0	5.1
2004.1733	-5.3	4.3
2004.2418	-3.7	3.9
2004.2445	-4.4	4.9
2004.2472	-4.3	4.5
2004.2500	-5.3	3.3
2004.2608	-4.8	3.8
2004.2636	-6.6	3.4
2004.3402	-4.9	3.2
2004.3430	-2.5	3.2
2004.6058	-6.2	3.1
2004.6085	-8.3	3.2
2004.6113	-6.8	3.1
2005.0522	0.1	3.1
2005.0549	-1.3	3.2
2005.1261	-3.6	3.0
2005.1289	-1.5	3.3
2005.1891	-5.1	3.2
2005.2028	-1.9	3.2
2005.2055	-3.7	3.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> diab.ps
1999.0288	8.8	5.6
1999.0315	2.9	5.8
1999.0343	-6.4	5.8
1999.0370	3.8	5.9
1999.0507	9.8	5.9
1999.0534	7.0	6.7
1999.0562	4.5	6.6
1999.0589	9.7	6.8
2001.0795	4.8	5.2
2001.0822	1.4	3.3
2001.0850	7.6	6.1
2001.0877	8.3	5.2
2001.0987	3.0	5.5
2001.1014	4.9	5.7
2001.1041	4.8	5.5
2001.1069	7.8	5.8
2001.2410	4.0	6.7
2002.1555	6.4	3.0
2002.1582	4.7	5.9
2002.1610	1.2	5.7
2004.1679	-7.2	5.3
2004.1706	-5.0	5.1
2004.1733	-5.3	4.3
2004.2418	-3.7	3.9
2004.2445	-4.4	4.9
2004.2472	-4.3	4.5
2004.2500	-5.3	3.3
2004.2608	-4.8	3.8
2004.2636	-6.6	3.4
2004.3402	-4.9	3.2
2004.3430	-2.5	3.2
2004.6058	-6.2	3.1
2004.6085	-8.3	3.2
2004.6113	-6.8	3.1
2005.0522	0.1	3.1
2005.0549	-1.3	3.2
2005.1261	-3.6	3.0
2005.1289	-1.5	3.3
2005.1891	-5.1	3.2
2005.2028	-1.9	3.2
2005.2055	-3.7	3.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> diab.ps
SBreak	Segment0	Comp0
 1999.0288	   6.8905	Segment0	Comp0
 2005.2055	  -4.8734	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1999/2005.5/-60/110 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"DIAB East (mm)":WeSn  \
-O >> diab.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> diab.ps
1999.0288	90.8	7.3
1999.0315	93.3	6.9
1999.0343	84.2	7.1
1999.0370	95.1	7.6
1999.0507	100.0	7.6
1999.0534	104.8	9.0
1999.0562	104.8	9.9
1999.0589	99.2	8.1
2001.0795	42.2	7.5
2001.0822	42.9	9.0
2001.0850	42.1	7.7
2001.0877	44.1	8.0
2001.0987	45.4	7.5
2001.1014	45.3	7.4
2001.1041	33.5	6.7
2001.1069	39.9	7.5
2001.2410	49.5	9.3
2002.1555	18.1	4.9
2002.1582	11.7	7.7
2002.1610	11.7	7.3
2004.1679	-27.9	6.4
2004.1706	-27.2	6.2
2004.1733	-28.4	7.1
2004.2418	-30.5	4.9
2004.2445	-32.1	7.4
2004.2472	-31.7	6.2
2004.2500	-31.6	4.9
2004.2608	-33.1	5.2
2004.2636	-29.5	5.0
2004.3402	-36.1	5.3
2004.3430	-32.5	4.6
2004.6058	-41.6	4.5
2004.6085	-40.4	4.7
2004.6113	-43.1	4.9
2005.0522	-47.6	4.7
2005.0549	-48.4	6.0
2005.1261	-53.6	4.3
2005.1289	-54.2	5.3
2005.1891	-50.7	4.6
2005.2028	-47.6	4.6
2005.2055	-53.2	4.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> diab.ps
1999.0288	90.8	7.3
1999.0315	93.3	6.9
1999.0343	84.2	7.1
1999.0370	95.1	7.6
1999.0507	100.0	7.6
1999.0534	104.8	9.0
1999.0562	104.8	9.9
1999.0589	99.2	8.1
2001.0795	42.2	7.5
2001.0822	42.9	9.0
2001.0850	42.1	7.7
2001.0877	44.1	8.0
2001.0987	45.4	7.5
2001.1014	45.3	7.4
2001.1041	33.5	6.7
2001.1069	39.9	7.5
2001.2410	49.5	9.3
2002.1555	18.1	4.9
2002.1582	11.7	7.7
2002.1610	11.7	7.3
2004.1679	-27.9	6.4
2004.1706	-27.2	6.2
2004.1733	-28.4	7.1
2004.2418	-30.5	4.9
2004.2445	-32.1	7.4
2004.2472	-31.7	6.2
2004.2500	-31.6	4.9
2004.2608	-33.1	5.2
2004.2636	-29.5	5.0
2004.3402	-36.1	5.3
2004.3430	-32.5	4.6
2004.6058	-41.6	4.5
2004.6085	-40.4	4.7
2004.6113	-43.1	4.9
2005.0522	-47.6	4.7
2005.0549	-48.4	6.0
2005.1261	-53.6	4.3
2005.1289	-54.2	5.3
2005.1891	-50.7	4.6
2005.2028	-47.6	4.6
2005.2055	-53.2	4.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> diab.ps
SBreak	Segment0	Comp1
 1999.0288	  93.7912	Segment0	Comp1
 2005.2055	 -54.2784	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1999/2005.5/-30/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"DIAB Up (mm)":WeSn  \
-O >> diab.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> diab.ps
1999.0288	5.6	22.8
1999.0315	-0.4	20.4
1999.0343	-7.2	21.0
1999.0370	-5.8	20.3
1999.0507	-6.0	20.5
1999.0534	9.8	24.0
1999.0562	7.2	23.8
1999.0589	18.0	22.5
2001.0795	-13.0	19.7
2001.0822	2.7	16.0
2001.0850	-11.0	21.1
2001.0877	-10.4	19.6
2001.0987	-13.2	21.3
2001.1014	-1.8	23.1
2001.1041	0.3	19.0
2001.1069	-8.2	20.1
2001.2410	-6.7	26.2
2002.1555	-11.9	10.9
2002.1582	-21.3	21.7
2002.1610	0.7	22.0
2004.1679	12.8	18.5
2004.1706	14.6	17.4
2004.1733	11.9	16.1
2004.2418	8.7	14.1
2004.2445	4.2	19.1
2004.2472	6.8	16.3
2004.2500	3.1	11.7
2004.2608	18.8	13.2
2004.2636	10.4	12.6
2004.3402	24.1	12.1
2004.3430	11.2	11.7
2004.6058	-2.8	11.7
2004.6085	-2.9	12.1
2004.6113	-8.7	11.6
2005.0522	-16.9	11.7
2005.0549	-18.7	12.4
2005.1261	-4.6	11.2
2005.1289	-14.2	12.3
2005.1891	3.9	11.9
2005.2028	4.8	11.7
2005.2055	4.2	11.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> diab.ps
1999.0288	5.6	22.8
1999.0315	-0.4	20.4
1999.0343	-7.2	21.0
1999.0370	-5.8	20.3
1999.0507	-6.0	20.5
1999.0534	9.8	24.0
1999.0562	7.2	23.8
1999.0589	18.0	22.5
2001.0795	-13.0	19.7
2001.0822	2.7	16.0
2001.0850	-11.0	21.1
2001.0877	-10.4	19.6
2001.0987	-13.2	21.3
2001.1014	-1.8	23.1
2001.1041	0.3	19.0
2001.1069	-8.2	20.1
2001.2410	-6.7	26.2
2002.1555	-11.9	10.9
2002.1582	-21.3	21.7
2002.1610	0.7	22.0
2004.1679	12.8	18.5
2004.1706	14.6	17.4
2004.1733	11.9	16.1
2004.2418	8.7	14.1
2004.2445	4.2	19.1
2004.2472	6.8	16.3
2004.2500	3.1	11.7
2004.2608	18.8	13.2
2004.2636	10.4	12.6
2004.3402	24.1	12.1
2004.3430	11.2	11.7
2004.6058	-2.8	11.7
2004.6085	-2.9	12.1
2004.6113	-8.7	11.6
2005.0522	-16.9	11.7
2005.0549	-18.7	12.4
2005.1261	-4.6	11.2
2005.1289	-14.2	12.3
2005.1891	3.9	11.9
2005.2028	4.8	11.7
2005.2055	4.2	11.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> diab.ps
SBreak	Segment0	Comp2
 1999.0288	  -2.1732	Segment0	Comp2
 2005.2055	   1.4882	Segment0	Comp2
EOT