#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2004.5/-60/80 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"DBLO North (mm)":WeSn  \
> dblo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dblo.ps
1993.0547	-55.2	5.6
1994.3167	-45.6	5.8
1995.4858	-37.2	5.8
1995.8473	-36.9	5.9
1995.8501	-34.7	5.9
1996.4414	-22.9	6.7
1997.9391	-12.6	6.6
1997.9583	-7.6	6.8
1999.0507	18.3	5.2
1999.0534	8.1	3.3
2001.0987	40.2	6.1
2001.1014	42.4	5.2
2004.2445	74.6	5.5
2004.2636	74.9	5.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dblo.ps
1993.0547	-55.2	5.6
1994.3167	-45.6	5.8
1995.4858	-37.2	5.8
1995.8473	-36.9	5.9
1995.8501	-34.7	5.9
1996.4414	-22.9	6.7
1997.9391	-12.6	6.6
1997.9583	-7.6	6.8
1999.0507	18.3	5.2
1999.0534	8.1	3.3
2001.0987	40.2	6.1
2001.1014	42.4	5.2
2004.2445	74.6	5.5
2004.2636	74.9	5.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> dblo.ps
SBreak	Segment0	Comp0
 1993.0547	 -62.7426	Segment0	Comp0
 2004.2636	  74.0377	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-80/90 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"DBLO East (mm)":WeSn  \
-O >> dblo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dblo.ps
1993.0547	84.5	7.3
1994.3167	38.2	6.9
1995.4858	16.0	7.1
1995.8473	33.8	7.6
1995.8501	32.7	7.6
1996.4414	38.8	9.0
1997.9391	12.6	9.9
1997.9583	0.1	8.1
1999.0507	3.6	7.5
1999.0534	5.7	9.0
2001.0987	-23.9	7.7
2001.1014	-31.9	8.0
2004.2445	-70.1	7.5
2004.2636	-70.2	7.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dblo.ps
1993.0547	84.5	7.3
1994.3167	38.2	6.9
1995.4858	16.0	7.1
1995.8473	33.8	7.6
1995.8501	32.7	7.6
1996.4414	38.8	9.0
1997.9391	12.6	9.9
1997.9583	0.1	8.1
1999.0507	3.6	7.5
1999.0534	5.7	9.0
2001.0987	-23.9	7.7
2001.1014	-31.9	8.0
2004.2445	-70.1	7.5
2004.2636	-70.2	7.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> dblo.ps
SBreak	Segment0	Comp1
 1993.0547	  71.9067	Segment0	Comp1
 2004.2636	 -68.5085	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-70/120 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"DBLO Up (mm)":WeSn  \
-O >> dblo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dblo.ps
1993.0547	51.8	22.8
1994.3167	-17.6	20.4
1995.4858	-25.0	21.0
1995.8473	86.3	20.3
1995.8501	79.5	20.5
1996.4414	112.2	24.0
1997.9391	-29.2	23.8
1997.9583	-46.2	22.5
1999.0507	-64.7	19.7
1999.0534	-8.8	16.0
2001.0987	-58.8	21.1
2001.1014	-44.4	19.6
2004.2445	-18.5	21.3
2004.2636	-17.6	23.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dblo.ps
1993.0547	51.8	22.8
1994.3167	-17.6	20.4
1995.4858	-25.0	21.0
1995.8473	86.3	20.3
1995.8501	79.5	20.5
1996.4414	112.2	24.0
1997.9391	-29.2	23.8
1997.9583	-46.2	22.5
1999.0507	-64.7	19.7
1999.0534	-8.8	16.0
2001.0987	-58.8	21.1
2001.1014	-44.4	19.6
2004.2445	-18.5	21.3
2004.2636	-17.6	23.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> dblo.ps
SBreak	Segment0	Comp2
 1993.0547	  45.4444	Segment0	Comp2
 2004.2636	 -45.3389	Segment0	Comp2
EOT