#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1996/2010/-20/40 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"VAC3 North (mm)":WeSn  \
> vac3.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> vac3.ps
1996.8330	32.0	8.6
1997.8871	31.2	6.1
1997.9227	19.5	6.3
1998.9220	24.2	6.5
1998.9576	28.1	6.7
2001.0494	22.3	6.1
2001.2355	17.4	5.7
2003.1001	14.6	5.9
2003.1110	18.1	5.8
2006.1036	3.3	3.0
2006.1063	2.4	3.0
2006.1090	-0.8	3.6
2009.1891	-15.7	5.9
2009.1918	-17.5	2.9
2009.1945	-16.5	2.9
2009.1973	-14.7	3.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> vac3.ps
1996.8330	32.0	8.6
1997.8871	31.2	6.1
1997.9227	19.5	6.3
1998.9220	24.2	6.5
1998.9576	28.1	6.7
2001.0494	22.3	6.1
2001.2355	17.4	5.7
2003.1001	14.6	5.9
2003.1110	18.1	5.8
2006.1036	3.3	3.0
2006.1063	2.4	3.0
2006.1090	-0.8	3.6
2009.1891	-15.7	5.9
2009.1918	-17.5	2.9
2009.1945	-16.5	2.9
2009.1973	-14.7	3.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> vac3.ps
SBreak	Segment0	Comp0
 1996.8330	  36.3516	Segment0	Comp0
 2009.1973	 -13.8938	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-80/200 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"VAC3 East (mm)":WeSn  \
-O >> vac3.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> vac3.ps
1996.8330	196.2	10.6
1997.8871	191.6	7.7
1997.9227	181.9	8.0
1998.9220	161.9	8.2
1998.9576	161.4	8.9
2001.0494	109.3	7.1
2001.2355	110.2	8.7
2003.1001	66.2	7.5
2003.1110	76.6	8.8
2006.1036	-0.3	3.5
2006.1063	-1.8	3.4
2006.1090	-1.4	3.8
2009.1891	-73.3	6.2
2009.1918	-72.1	3.4
2009.1945	-70.2	3.4
2009.1973	-72.1	4.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> vac3.ps
1996.8330	196.2	10.6
1997.8871	191.6	7.7
1997.9227	181.9	8.0
1998.9220	161.9	8.2
1998.9576	161.4	8.9
2001.0494	109.3	7.1
2001.2355	110.2	8.7
2003.1001	66.2	7.5
2003.1110	76.6	8.8
2006.1036	-0.3	3.5
2006.1063	-1.8	3.4
2006.1090	-1.4	3.8
2009.1891	-73.3	6.2
2009.1918	-72.1	3.4
2009.1945	-70.2	3.4
2009.1973	-72.1	4.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> vac3.ps
SBreak	Segment0	Comp1
 1996.8330	 208.7887	Segment0	Comp1
 2009.1973	 -71.3627	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-30/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"VAC3 Up (mm)":WeSn  \
-O >> vac3.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> vac3.ps
1996.8330	16.1	34.1
1997.8871	-17.3	21.5
1997.9227	-20.8	22.1
1998.9220	-4.3	23.1
1998.9576	-17.5	24.2
2001.0494	-11.8	22.3
2001.2355	5.5	22.9
2003.1001	-6.0	21.0
2003.1110	35.3	22.5
2006.1036	-0.7	10.6
2006.1063	9.3	10.7
2006.1090	-11.4	12.5
2009.1891	-10.7	20.7
2009.1918	2.7	10.7
2009.1945	5.5	10.7
2009.1973	2.1	13.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> vac3.ps
1996.8330	16.1	34.1
1997.8871	-17.3	21.5
1997.9227	-20.8	22.1
1998.9220	-4.3	23.1
1998.9576	-17.5	24.2
2001.0494	-11.8	22.3
2001.2355	5.5	22.9
2003.1001	-6.0	21.0
2003.1110	35.3	22.5
2006.1036	-0.7	10.6
2006.1063	9.3	10.7
2006.1090	-11.4	12.5
2009.1891	-10.7	20.7
2009.1918	2.7	10.7
2009.1945	5.5	10.7
2009.1973	2.1	13.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> vac3.ps
SBreak	Segment0	Comp2
 1996.8330	  -8.6661	Segment0	Comp2
 2009.1973	   3.3235	Segment0	Comp2
EOT