#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-120/170 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-106.9	4.2
1993.0355	-114.4	4.2
1994.3003	-82.4	4.0
1995.7131	-55.7	6.2
1996.8248	-43.0	6.7
1997.8789	-39.0	6.3
1997.9172	-29.4	6.2
1998.9138	-3.7	6.8
1998.9357	-11.1	5.9
2001.0247	33.6	6.2
2001.0275	26.5	6.7
2006.0899	107.8	5.1
2006.1063	122.1	5.0
2009.1727	157.4	4.9
2009.1755	160.5	5.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	-106.9	4.2
1993.0355	-114.4	4.2
1994.3003	-82.4	4.0
1995.7131	-55.7	6.2
1996.8248	-43.0	6.7
1997.8789	-39.0	6.3
1997.9172	-29.4	6.2
1998.9138	-3.7	6.8
1998.9357	-11.1	5.9
2001.0247	33.6	6.2
2001.0275	26.5	6.7
2006.0899	107.8	5.1
2006.1063	122.1	5.0
2009.1755	160.5	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	-107.8169	Segment0	Comp0
 2009.1755	 162.8323	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-270/300 \
-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	291.5	7.9
1993.0355	283.3	9.7
1994.3003	237.9	6.1
1995.7131	196.1	8.3
1996.8248	153.7	8.2
1997.8789	118.5	7.9
1997.9172	120.3	7.7
1998.9138	93.3	8.8
1998.9357	80.0	7.6
2001.0247	21.3	7.1
2001.0275	9.5	7.2
2006.0899	-162.6	5.9
2006.1063	-155.7	5.9
2009.1727	-267.0	5.4
2009.1755	-261.2	5.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	291.5	7.9
1993.0355	283.3	9.7
1994.3003	237.9	6.1
1995.7131	196.1	8.3
1996.8248	153.7	8.2
1997.8789	118.5	7.9
1997.9172	120.3	7.7
1998.9138	93.3	8.8
1998.9357	80.0	7.6
2001.0247	21.3	7.1
2001.0275	9.5	7.2
2006.0899	-162.6	5.9
2006.1063	-155.7	5.9
2009.1755	-261.2	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	 285.4637	Segment0	Comp1
 2009.1755	-263.8849	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	25.2	21.7
1993.0355	46.4	23.3
1994.3003	40.0	20.2
1995.7131	-29.8	22.4
1996.8248	-5.1	24.0
1997.8789	-17.4	22.7
1997.9172	8.7	21.8
1998.9138	-12.0	24.2
1998.9357	-7.0	21.7
2001.0247	-32.7	22.6
2001.0275	-35.1	23.0
2006.0899	-1.5	19.2
2006.1063	-0.7	18.4
2009.1727	-3.8	17.1
2009.1755	10.5	18.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 1395.ps
1993.0328	25.2	21.7
1993.0355	46.4	23.3
1994.3003	40.0	20.2
1995.7131	-29.8	22.4
1996.8248	-5.1	24.0
1997.8789	-17.4	22.7
1997.9172	8.7	21.8
1998.9138	-12.0	24.2
1998.9357	-7.0	21.7
2001.0247	-32.7	22.6
2001.0275	-35.1	23.0
2006.0899	-1.5	19.2
2006.1063	-0.7	18.4
2009.1755	10.5	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.4727	Segment0	Comp2
 2009.1755	  -8.9832	Segment0	Comp2
EOT