#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-160/90 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-152.1	3.7
1994.2948	-132.9	4.0
1995.7158	-112.0	5.9
1996.8275	-92.9	5.0
1997.8789	-74.1	6.2
1997.9172	-76.0	6.1
1998.9138	-65.9	6.4
1998.9330	-61.8	6.1
2001.0220	-17.6	6.2
2001.0247	-20.8	6.0
2006.0871	42.7	2.8
2006.0899	43.1	2.8
2009.1700	88.4	5.1
2009.1727	88.0	2.9
2009.1755	90.0	3.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	-152.1	3.7
1994.2948	-132.9	4.0
1995.7158	-112.0	5.9
1996.8275	-92.9	5.0
1997.8789	-74.1	6.2
1997.9172	-76.0	6.1
1998.9138	-65.9	6.4
1998.9330	-61.8	6.1
2001.0220	-17.6	6.2
2001.0247	-20.8	6.0
2006.0871	42.7	2.8
2006.0899	43.1	2.8
2009.1755	90.0	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	-149.5035	Segment0	Comp0
 2009.1755	  89.4666	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-160/360 \
-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	354.2	10.7
1994.2948	318.8	6.1
1995.7158	283.7	7.6
1996.8275	246.4	6.3
1997.8789	216.0	7.6
1997.9172	213.3	7.6
1998.9138	185.1	8.3
1998.9330	186.0	7.5
2001.0220	118.0	8.0
2001.0247	115.5	6.6
2006.0871	-48.2	3.1
2006.0899	-51.4	3.1
2009.1700	-153.2	5.6
2009.1727	-153.3	3.3
2009.1755	-149.6	3.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	354.2	10.7
1994.2948	318.8	6.1
1995.7158	283.7	7.6
1996.8275	246.4	6.3
1997.8789	216.0	7.6
1997.9172	213.3	7.6
1998.9138	185.1	8.3
1998.9330	186.0	7.5
2001.0220	118.0	8.0
2001.0247	115.5	6.6
2006.0871	-48.2	3.1
2006.0899	-51.4	3.1
2009.1755	-149.6	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	 368.2389	Segment0	Comp1
 2009.1755	-150.0914	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	79.4	20.0
1994.2948	48.1	17.3
1995.7158	30.9	20.3
1996.8275	-23.5	18.4
1997.8789	-0.0	22.6
1997.9172	3.7	21.3
1998.9138	3.0	21.4
1998.9330	-31.4	21.3
2001.0220	-2.9	22.5
2001.0247	-29.5	20.3
2006.0871	-5.8	9.8
2006.0899	4.4	9.8
2009.1700	2.7	17.4
2009.1727	-12.9	10.5
2009.1755	-12.5	10.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> nicc.ps
1993.0383	79.4	20.0
1994.2948	48.1	17.3
1995.7158	30.9	20.3
1996.8275	-23.5	18.4
1997.8789	-0.0	22.6
1997.9172	3.7	21.3
1998.9138	3.0	21.4
1998.9330	-31.4	21.3
2001.0220	-2.9	22.5
2001.0247	-29.5	20.3
2006.0871	-5.8	9.8
2006.0899	4.4	9.8
2009.1755	-12.5	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.7404	Segment0	Comp2
 2009.1755	 -12.6381	Segment0	Comp2
EOT