#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1996/2010/-90/40 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"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	-88.9	8.6
1997.8871	-75.5	6.1
1997.9227	-86.8	6.3
1998.9220	-68.4	6.5
1998.9576	-64.1	6.7
2001.0494	-41.6	6.1
2001.2355	-44.0	5.7
2003.1001	-21.6	5.9
2003.1110	-18.0	5.8
2006.1036	7.7	3.0
2006.1063	6.9	3.0
2006.1090	3.7	3.6
2009.1891	30.4	5.9
2009.1918	28.6	2.9
2009.1945	29.8	2.9
2009.1973	31.4	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	-88.9	8.6
1997.8871	-75.5	6.1
1997.9227	-86.8	6.3
1998.9220	-68.4	6.5
1998.9576	-64.1	6.7
2001.0494	-41.6	6.1
2001.2355	-44.0	5.7
2003.1001	-21.6	5.9
2003.1110	-18.0	5.8
2006.1036	7.7	3.0
2006.1063	6.9	3.0
2006.1090	3.7	3.6
2009.1891	30.4	5.9
2009.1918	28.6	2.9
2009.1945	29.8	2.9
2009.1973	31.4	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	 -84.5408	Segment0	Comp0
 2009.1973	  32.3119	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-40/100 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"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	87.3	10.6
1997.8871	95.1	7.7
1997.9227	85.9	8.0
1998.9220	77.6	8.2
1998.9576	77.6	8.9
2001.0494	50.2	7.1
2001.2355	53.3	8.7
2003.1001	31.4	7.5
2003.1110	41.9	8.8
2006.1036	0.3	3.5
2006.1063	-1.1	3.4
2006.1090	-0.7	3.8
2009.1891	-36.2	6.2
2009.1918	-34.9	3.4
2009.1945	-33.0	3.4
2009.1973	-34.9	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	87.3	10.6
1997.8871	95.1	7.7
1997.9227	85.9	8.0
1998.9220	77.6	8.2
1998.9576	77.6	8.9
2001.0494	50.2	7.1
2001.2355	53.3	8.7
2003.1001	31.4	7.5
2003.1110	41.9	8.8
2006.1036	0.3	3.5
2006.1063	-1.1	3.4
2006.1090	-0.7	3.8
2009.1891	-36.2	6.2
2009.1918	-34.9	3.4
2009.1945	-33.0	3.4
2009.1973	-34.9	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	  99.8551	Segment0	Comp1
 2009.1973	 -34.1298	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	15.7	34.1
1997.8871	-17.6	21.5
1997.9227	-21.2	22.1
1998.9220	-4.6	23.1
1998.9576	-17.8	24.2
2001.0494	-12.0	22.3
2001.2355	5.4	22.9
2003.1001	-6.1	21.0
2003.1110	35.2	22.5
2006.1036	-0.7	10.6
2006.1063	9.3	10.7
2006.1090	-11.3	12.5
2009.1891	-10.6	20.7
2009.1918	2.8	10.7
2009.1945	5.7	10.7
2009.1973	2.3	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	15.7	34.1
1997.8871	-17.6	21.5
1997.9227	-21.2	22.1
1998.9220	-4.6	23.1
1998.9576	-17.8	24.2
2001.0494	-12.0	22.3
2001.2355	5.4	22.9
2003.1001	-6.1	21.0
2003.1110	35.2	22.5
2006.1036	-0.7	10.6
2006.1063	9.3	10.7
2006.1090	-11.3	12.5
2009.1891	-10.6	20.7
2009.1918	2.8	10.7
2009.1945	5.7	10.7
2009.1973	2.3	13.6
EOT


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