#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994/2004.5/-20/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"ROC2 North (mm)":WeSn  \
> roc2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> roc2.ps
1994.3167	-15.1	5.6
1995.4886	-2.9	5.8
1996.4414	4.2	5.8
1997.9418	-4.6	5.9
1997.9637	-13.8	5.9
1999.0507	-0.8	6.7
1999.0534	-16.9	6.6
2001.1014	13.8	6.8
2001.1069	7.0	5.2
2004.2472	15.7	3.3
2004.2500	13.2	6.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> roc2.ps
1994.3167	-15.1	5.6
1995.4886	-2.9	5.8
1996.4414	4.2	5.8
1997.9418	-4.6	5.9
1997.9637	-13.8	5.9
1999.0507	-0.8	6.7
1999.0534	-16.9	6.6
2001.1014	13.8	6.8
2001.1069	7.0	5.2
2004.2472	15.7	3.3
2004.2500	13.2	6.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> roc2.ps
SBreak	Segment0	Comp0
 1994.3167	 -12.6653	Segment0	Comp0
 2004.2500	  12.6930	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994/2004.5/-130/130 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"ROC2 East (mm)":WeSn  \
-O >> roc2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> roc2.ps
1994.3167	125.9	7.3
1995.4886	99.4	6.9
1996.4414	84.7	7.1
1997.9418	46.0	7.6
1997.9637	42.6	7.6
1999.0507	17.9	9.0
1999.0534	20.4	9.9
2001.1014	-34.5	8.1
2001.1069	-46.2	7.5
2004.2472	-126.5	9.0
2004.2500	-126.6	7.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> roc2.ps
1994.3167	125.9	7.3
1995.4886	99.4	6.9
1996.4414	84.7	7.1
1997.9418	46.0	7.6
1997.9637	42.6	7.6
1999.0507	17.9	9.0
1999.0534	20.4	9.9
2001.1014	-34.5	8.1
2001.1069	-46.2	7.5
2004.2472	-126.5	9.0
2004.2500	-126.6	7.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> roc2.ps
SBreak	Segment0	Comp1
 1994.3167	 137.1583	Segment0	Comp1
 2004.2500	-124.1731	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1994/2004.5/-50/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"ROC2 Up (mm)":WeSn  \
-O >> roc2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> roc2.ps
1994.3167	21.4	22.8
1995.4886	15.6	20.4
1996.4414	27.1	21.0
1997.9418	-5.3	20.3
1997.9637	0.1	20.5
1999.0507	-44.1	24.0
1999.0534	11.1	23.8
2001.1014	-8.9	22.5
2001.1069	-20.5	19.7
2004.2472	11.1	16.0
2004.2500	0.2	21.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> roc2.ps
1994.3167	21.4	22.8
1995.4886	15.6	20.4
1996.4414	27.1	21.0
1997.9418	-5.3	20.3
1997.9637	0.1	20.5
1999.0507	-44.1	24.0
1999.0534	11.1	23.8
2001.1014	-8.9	22.5
2001.1069	-20.5	19.7
2004.2472	11.1	16.0
2004.2500	0.2	21.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> roc2.ps
SBreak	Segment0	Comp2
 1994.3167	   7.6554	Segment0	Comp2
 2004.2500	  -7.1401	Segment0	Comp2
EOT