#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1998/2010/-130/130 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"PTRB North (mm)":WeSn  \
> ptrb.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ptrb.ps
1998.9138	-114.4	3.0
1998.9165	-117.9	3.0
1998.9193	-117.1	2.9
1998.9220	-121.5	2.9
1998.9330	-114.8	2.9
1998.9357	-113.5	2.9
1998.9384	-116.7	2.9
1998.9412	-113.3	3.0
1998.9576	-113.8	2.9
2001.0220	-69.1	3.0
2001.0247	-67.9	2.9
2001.0275	-70.7	3.0
2001.0302	-67.4	3.1
2001.0439	-63.3	3.1
2001.0466	-63.4	3.0
2001.0494	-69.3	3.0
2001.0685	-62.9	3.1
2001.1041	-63.6	3.0
2001.1069	-60.7	3.0
2001.2383	-59.8	3.1
2003.0590	-19.3	3.1
2003.0617	-22.0	3.2
2003.0754	-18.7	3.2
2003.0782	-20.5	3.2
2003.0919	-20.4	3.1
2003.0946	-18.9	3.1
2003.0973	-18.1	3.1
2003.1001	-18.0	3.1
2003.1110	-15.3	3.1
2003.1138	-18.4	3.2
2003.1165	-16.6	3.2
2006.0844	48.6	2.9
2006.0871	48.2	2.9
2006.0899	48.7	2.8
2006.0926	50.6	2.9
2006.1008	54.0	3.0
2006.1036	54.6	2.9
2006.1063	55.5	2.9
2006.1090	54.9	2.9
2006.1227	51.0	2.9
2006.1255	49.7	2.9
2006.1282	50.4	2.9
2006.1309	50.3	2.9
2006.1446	51.8	2.9
2006.1474	51.4	2.9
2009.1700	120.8	2.9
2009.1727	118.3	3.0
2009.1755	120.2	2.9
2009.1810	121.5	2.9
2009.1836	118.4	2.8
2009.1863	119.7	2.8
2009.1891	121.7	2.8
2009.1918	120.7	2.8
2009.1945	121.9	2.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> ptrb.ps
2009.1700	120.8	2.9
2009.1727	118.3	3.0
2009.1836	118.4	2.8
2009.1863	119.7	2.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ptrb.ps
1998.9138	-114.4	3.0
1998.9165	-117.9	3.0
1998.9193	-117.1	2.9
1998.9220	-121.5	2.9
1998.9330	-114.8	2.9
1998.9357	-113.5	2.9
1998.9384	-116.7	2.9
1998.9412	-113.3	3.0
1998.9576	-113.8	2.9
2001.0220	-69.1	3.0
2001.0247	-67.9	2.9
2001.0275	-70.7	3.0
2001.0302	-67.4	3.1
2001.0439	-63.3	3.1
2001.0466	-63.4	3.0
2001.0494	-69.3	3.0
2001.0685	-62.9	3.1
2001.1041	-63.6	3.0
2001.1069	-60.7	3.0
2001.2383	-59.8	3.1
2003.0590	-19.3	3.1
2003.0617	-22.0	3.2
2003.0754	-18.7	3.2
2003.0782	-20.5	3.2
2003.0919	-20.4	3.1
2003.0946	-18.9	3.1
2003.0973	-18.1	3.1
2003.1001	-18.0	3.1
2003.1110	-15.3	3.1
2003.1138	-18.4	3.2
2003.1165	-16.6	3.2
2006.0844	48.6	2.9
2006.0871	48.2	2.9
2006.0899	48.7	2.8
2006.0926	50.6	2.9
2006.1008	54.0	3.0
2006.1036	54.6	2.9
2006.1063	55.5	2.9
2006.1090	54.9	2.9
2006.1227	51.0	2.9
2006.1255	49.7	2.9
2006.1282	50.4	2.9
2006.1309	50.3	2.9
2006.1446	51.8	2.9
2006.1474	51.4	2.9
2009.1755	120.2	2.9
2009.1810	121.5	2.9
2009.1891	121.7	2.8
2009.1918	120.7	2.8
2009.1945	121.9	2.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> ptrb.ps
SBreak	Segment0	Comp0
 1998.9138	-115.3910	Segment0	Comp0
 2009.1945	 121.6932	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1998/2010/-180/220 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"PTRB East (mm)":WeSn  \
-O >> ptrb.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ptrb.ps
1998.9138	217.0	4.8
1998.9165	217.6	4.7
1998.9193	214.8	4.5
1998.9220	205.9	4.3
1998.9330	217.5	4.6
1998.9357	214.0	4.4
1998.9384	216.9	5.0
1998.9412	218.7	5.7
1998.9576	212.5	4.6
2001.0220	140.1	4.5
2001.0247	142.1	4.0
2001.0275	134.2	4.0
2001.0302	138.0	4.4
2001.0439	130.3	4.3
2001.0466	138.3	4.4
2001.0494	128.6	4.0
2001.0685	136.5	4.6
2001.1041	123.2	4.5
2001.1069	132.1	4.4
2001.2383	134.1	5.1
2003.0590	61.4	4.6
2003.0617	58.4	4.7
2003.0754	65.8	4.8
2003.0782	70.0	4.8
2003.0919	60.2	4.4
2003.0946	60.9	4.2
2003.0973	59.6	4.4
2003.1001	53.7	4.3
2003.1110	61.3	4.7
2003.1138	64.5	4.8
2003.1165	61.0	4.7
2006.0844	-58.3	3.5
2006.0871	-55.6	3.3
2006.0899	-57.0	3.3
2006.0926	-55.3	3.5
2006.1008	-53.5	3.7
2006.1036	-55.9	3.4
2006.1063	-55.6	3.2
2006.1090	-58.7	3.2
2006.1227	-58.0	3.3
2006.1255	-57.5	3.4
2006.1282	-51.7	3.3
2006.1309	-55.3	3.5
2006.1446	-56.1	3.5
2006.1474	-59.2	3.4
2009.1700	-173.0	3.6
2009.1727	-176.1	3.5
2009.1755	-172.2	3.5
2009.1810	-171.2	3.7
2009.1836	-172.3	3.5
2009.1863	-173.4	3.5
2009.1891	-173.4	3.3
2009.1918	-173.3	3.3
2009.1945	-172.8	3.4
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> ptrb.ps
2009.1700	-173.0	3.6
2009.1727	-176.1	3.5
2009.1836	-172.3	3.5
2009.1863	-173.4	3.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ptrb.ps
1998.9138	217.0	4.8
1998.9165	217.6	4.7
1998.9193	214.8	4.5
1998.9220	205.9	4.3
1998.9330	217.5	4.6
1998.9357	214.0	4.4
1998.9384	216.9	5.0
1998.9412	218.7	5.7
1998.9576	212.5	4.6
2001.0220	140.1	4.5
2001.0247	142.1	4.0
2001.0275	134.2	4.0
2001.0302	138.0	4.4
2001.0439	130.3	4.3
2001.0466	138.3	4.4
2001.0494	128.6	4.0
2001.0685	136.5	4.6
2001.1041	123.2	4.5
2001.1069	132.1	4.4
2001.2383	134.1	5.1
2003.0590	61.4	4.6
2003.0617	58.4	4.7
2003.0754	65.8	4.8
2003.0782	70.0	4.8
2003.0919	60.2	4.4
2003.0946	60.9	4.2
2003.0973	59.6	4.4
2003.1001	53.7	4.3
2003.1110	61.3	4.7
2003.1138	64.5	4.8
2003.1165	61.0	4.7
2006.0844	-58.3	3.5
2006.0871	-55.6	3.3
2006.0899	-57.0	3.3
2006.0926	-55.3	3.5
2006.1008	-53.5	3.7
2006.1036	-55.9	3.4
2006.1063	-55.6	3.2
2006.1090	-58.7	3.2
2006.1227	-58.0	3.3
2006.1255	-57.5	3.4
2006.1282	-51.7	3.3
2006.1309	-55.3	3.5
2006.1446	-56.1	3.5
2006.1474	-59.2	3.4
2009.1755	-172.2	3.5
2009.1810	-171.2	3.7
2009.1891	-173.4	3.3
2009.1918	-173.3	3.3
2009.1945	-172.8	3.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> ptrb.ps
SBreak	Segment0	Comp1
 1998.9138	 216.5434	Segment0	Comp1
 2009.1945	-172.9451	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1998/2010/-50/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"PTRB Up (mm)":WeSn  \
-O >> ptrb.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ptrb.ps
1998.9138	-21.8	11.4
1998.9165	-25.7	11.5
1998.9193	-36.1	11.2
1998.9220	-38.2	11.0
1998.9330	-28.8	11.0
1998.9357	-33.4	11.0
1998.9384	-24.2	11.2
1998.9412	-30.0	11.6
1998.9576	-41.3	11.5
2001.0220	16.0	11.4
2001.0247	7.8	10.9
2001.0275	17.1	11.1
2001.0302	15.3	11.6
2001.0439	1.1	11.4
2001.0466	-4.2	11.4
2001.0494	16.3	11.1
2001.0685	-1.3	11.7
2001.1041	9.3	11.3
2001.1069	1.2	11.2
2001.2383	21.0	11.9
2003.0590	16.5	11.5
2003.0617	25.6	11.6
2003.0754	17.4	11.6
2003.0782	19.1	11.6
2003.0919	12.7	11.4
2003.0946	19.5	11.4
2003.0973	29.9	11.5
2003.1001	1.4	11.5
2003.1110	24.2	11.6
2003.1138	25.1	12.0
2003.1165	16.8	11.8
2006.0844	2.1	10.1
2006.0871	6.9	10.2
2006.0899	9.1	9.8
2006.0926	2.8	10.1
2006.1008	-3.1	10.2
2006.1036	-4.5	10.1
2006.1063	-5.1	10.2
2006.1090	-4.6	10.2
2006.1227	15.0	10.1
2006.1255	9.0	10.1
2006.1282	1.8	10.1
2006.1309	3.5	10.3
2006.1446	1.1	10.3
2006.1474	-3.7	10.3
2009.1700	-3.7	10.5
2009.1727	-3.7	10.6
2009.1755	-12.1	10.6
2009.1810	-9.2	10.2
2009.1836	-1.5	10.1
2009.1863	-4.3	10.1
2009.1891	-2.2	10.2
2009.1918	-5.5	10.2
2009.1945	-3.1	10.2
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> ptrb.ps
2009.1700	-3.7	10.5
2009.1727	-3.7	10.6
2009.1836	-1.5	10.1
2009.1863	-4.3	10.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ptrb.ps
1998.9138	-21.8	11.4
1998.9165	-25.7	11.5
1998.9193	-36.1	11.2
1998.9220	-38.2	11.0
1998.9330	-28.8	11.0
1998.9357	-33.4	11.0
1998.9384	-24.2	11.2
1998.9412	-30.0	11.6
1998.9576	-41.3	11.5
2001.0220	16.0	11.4
2001.0247	7.8	10.9
2001.0275	17.1	11.1
2001.0302	15.3	11.6
2001.0439	1.1	11.4
2001.0466	-4.2	11.4
2001.0494	16.3	11.1
2001.0685	-1.3	11.7
2001.1041	9.3	11.3
2001.1069	1.2	11.2
2001.2383	21.0	11.9
2003.0590	16.5	11.5
2003.0617	25.6	11.6
2003.0754	17.4	11.6
2003.0782	19.1	11.6
2003.0919	12.7	11.4
2003.0946	19.5	11.4
2003.0973	29.9	11.5
2003.1001	1.4	11.5
2003.1110	24.2	11.6
2003.1138	25.1	12.0
2003.1165	16.8	11.8
2006.0844	2.1	10.1
2006.0871	6.9	10.2
2006.0899	9.1	9.8
2006.0926	2.8	10.1
2006.1008	-3.1	10.2
2006.1036	-4.5	10.1
2006.1063	-5.1	10.2
2006.1090	-4.6	10.2
2006.1227	15.0	10.1
2006.1255	9.0	10.1
2006.1282	1.8	10.1
2006.1309	3.5	10.3
2006.1446	1.1	10.3
2006.1474	-3.7	10.3
2009.1755	-12.1	10.6
2009.1810	-9.2	10.2
2009.1891	-2.2	10.2
2009.1918	-5.5	10.2
2009.1945	-3.1	10.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> ptrb.ps
SBreak	Segment0	Comp2
 1998.9138	  -5.5853	Segment0	Comp2
 2009.1945	   5.4991	Segment0	Comp2
EOT