#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nicc.ps
1993.0383	-290.3	3.7
1994.2948	-253.8	4.0
1995.7158	-213.6	5.9
1996.8275	-179.3	5.0
1997.8789	-146.0	6.2
1997.9172	-147.4	6.1
1998.9138	-123.8	6.4
1998.9330	-119.3	6.1
2001.0220	-46.4	6.2
2001.0247	-49.6	6.0
2006.0871	83.1	2.8
2006.0899	83.6	2.8
2009.1700	171.0	5.1
2009.1727	170.7	2.9
2009.1755	172.7	3.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> nicc.ps
2009.1700	171.0	5.1
2009.1727	170.7	2.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	-290.3	3.7
1994.2948	-253.8	4.0
1995.7158	-213.6	5.9
1996.8275	-179.3	5.0
1997.8789	-146.0	6.2
1997.9172	-147.4	6.1
1998.9138	-123.8	6.4
1998.9330	-119.3	6.1
2001.0220	-46.4	6.2
2001.0247	-49.6	6.0
2006.0871	83.1	2.8
2006.0899	83.6	2.8
2009.1755	172.7	3.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> nicc.ps
SBreak	Segment0	Comp0
 1993.0383	-287.7111	Segment0	Comp0
 2009.1755	 172.1734	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nicc.ps
1993.0383	220.5	10.7
1994.2948	199.6	6.1
1995.7158	181.2	7.6
1996.8275	156.9	6.3
1997.8789	138.6	7.6
1997.9172	136.5	7.6
1998.9138	119.9	8.3
1998.9330	121.0	7.5
2001.0220	77.4	8.0
2001.0247	74.9	6.6
2006.0871	-29.7	3.1
2006.0899	-32.9	3.1
2009.1700	-98.9	5.6
2009.1727	-98.8	3.3
2009.1755	-95.0	3.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> nicc.ps
2009.1700	-98.9	5.6
2009.1727	-98.8	3.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	220.5	10.7
1994.2948	199.6	6.1
1995.7158	181.2	7.6
1996.8275	156.9	6.3
1997.8789	138.6	7.6
1997.9172	136.5	7.6
1998.9138	119.9	8.3
1998.9330	121.0	7.5
2001.0220	77.4	8.0
2001.0247	74.9	6.6
2006.0871	-29.7	3.1
2006.0899	-32.9	3.1
2009.1755	-95.0	3.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> nicc.ps
SBreak	Segment0	Comp1
 1993.0383	 234.4862	Segment0	Comp1
 2009.1755	 -95.5748	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nicc.ps
1993.0383	78.9	20.0
1994.2948	47.6	17.3
1995.7158	30.6	20.3
1996.8275	-23.8	18.4
1997.8789	-0.2	22.6
1997.9172	3.5	21.3
1998.9138	2.8	21.4
1998.9330	-31.7	21.3
2001.0220	-3.0	22.5
2001.0247	-29.6	20.3
2006.0871	-5.8	9.8
2006.0899	4.5	9.8
2009.1700	2.9	17.4
2009.1727	-12.6	10.5
2009.1755	-12.2	10.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G0/0/255 -O \
-K >> nicc.ps
2009.1700	2.9	17.4
2009.1727	-12.6	10.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	78.9	20.0
1994.2948	47.6	17.3
1995.7158	30.6	20.3
1996.8275	-23.8	18.4
1997.8789	-0.2	22.6
1997.9172	3.5	21.3
1998.9138	2.8	21.4
1998.9330	-31.7	21.3
2001.0220	-3.0	22.5
2001.0247	-29.6	20.3
2006.0871	-5.8	9.8
2006.0899	4.5	9.8
2009.1755	-12.2	10.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> nicc.ps
SBreak	Segment0	Comp2
 1993.0383	  25.2644	Segment0	Comp2
 2009.1755	 -12.4044	Segment0	Comp2
EOT