#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-210/300 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f80a160:"1395 North (mm)":WeSn  \
> 1395.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 1395.ps
1993.0328	-195.1	4.2
1993.0355	-202.5	4.2
1994.3003	-153.2	4.0
1995.7131	-107.2	6.2
1996.8248	-79.3	6.7
1997.8789	-60.7	6.3
1997.9172	-50.7	6.2
1998.9138	-11.2	6.8
1998.9357	-18.3	5.9
2001.0247	55.0	6.2
2001.0275	47.9	6.7
2006.0899	198.7	5.1
2006.1063	213.2	5.0
2009.1727	290.6	4.9
2009.1755	293.6	5.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	-195.1	4.2
1993.0355	-202.5	4.2
1994.3003	-153.2	4.0
1995.7131	-107.2	6.2
1996.8248	-79.3	6.7
1997.8789	-60.7	6.3
1997.9172	-50.7	6.2
1998.9138	-11.2	6.8
1998.9357	-18.3	5.9
2001.0247	55.0	6.2
2001.0275	47.9	6.7
2006.0899	198.7	5.1
2006.1063	213.2	5.0
2009.1755	293.6	5.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 1395.ps
SBreak	Segment0	Comp0
 1993.0328	-195.9986	Segment0	Comp0
 2009.1755	 296.0103	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 1395.ps
1993.0328	193.8	7.9
1993.0355	185.6	9.7
1994.3003	154.9	6.1
1995.7131	129.6	8.3
1996.8248	100.1	8.2
1997.8789	77.2	7.9
1997.9172	79.5	7.7
1998.9138	64.1	8.8
1998.9357	51.0	7.6
2001.0247	16.7	7.1
2001.0275	4.9	7.2
2006.0899	-108.2	5.9
2006.1063	-101.1	5.9
2009.1727	-176.7	5.4
2009.1755	-170.9	5.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	193.8	7.9
1993.0355	185.6	9.7
1994.3003	154.9	6.1
1995.7131	129.6	8.3
1996.8248	100.1	8.2
1997.8789	77.2	7.9
1997.9172	79.5	7.7
1998.9138	64.1	8.8
1998.9357	51.0	7.6
2001.0247	16.7	7.1
2001.0275	4.9	7.2
2006.0899	-108.2	5.9
2006.1063	-101.1	5.9
2009.1755	-170.9	5.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 1395.ps
SBreak	Segment0	Comp1
 1993.0328	 187.7054	Segment0	Comp1
 2009.1755	-173.5163	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-40/50 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"1395 Up (mm)":WeSn  \
-O >> 1395.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 1395.ps
1993.0328	24.9	21.7
1993.0355	46.0	23.3
1994.3003	39.6	20.2
1995.7131	-29.9	22.4
1996.8248	-5.3	24.0
1997.8789	-17.6	22.7
1997.9172	8.5	21.8
1998.9138	-12.1	24.2
1998.9357	-7.1	21.7
2001.0247	-32.7	22.6
2001.0275	-35.2	23.0
2006.0899	-1.2	19.2
2006.1063	-0.4	18.4
2009.1727	-3.5	17.1
2009.1755	10.9	18.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	24.9	21.7
1993.0355	46.0	23.3
1994.3003	39.6	20.2
1995.7131	-29.9	22.4
1996.8248	-5.3	24.0
1997.8789	-17.6	22.7
1997.9172	8.5	21.8
1998.9138	-12.1	24.2
1998.9357	-7.1	21.7
2001.0247	-32.7	22.6
2001.0275	-35.2	23.0
2006.0899	-1.2	19.2
2006.1063	-0.4	18.4
2009.1755	10.9	18.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> 1395.ps
SBreak	Segment0	Comp2
 1993.0328	   8.1279	Segment0	Comp2
 2009.1755	  -8.6176	Segment0	Comp2
EOT