#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"REF2 North (mm)":WeSn  \
> ref2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ref2.ps
2001.6187	-2.4	8.4
2002.4894	9.3	14.4
2002.4921	4.9	20.6
2002.4949	0.0	15.4
2002.4976	-0.7	21.8
2002.5003	-0.2	17.9
2004.4689	-2.1	9.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ref2.ps
2001.6187	-2.4	8.4
2002.4894	9.3	14.4
2002.4921	4.9	20.6
2002.4949	0.0	15.4
2002.4976	-0.7	21.8
2002.5003	-0.2	17.9
2004.4689	-2.1	9.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> ref2.ps
SBreak	Segment0	Comp0
 2001.6187	   0.5517	Segment0	Comp0
 2004.4689	  -0.8716	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"REF2 East (mm)":WeSn  \
-O >> ref2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ref2.ps
2001.6187	7.1	7.9
2002.4894	36.3	17.0
2002.4921	-1.0	23.2
2002.4949	1.0	17.6
2002.4976	-0.3	21.0
2002.5003	6.0	16.5
2004.4689	-6.7	10.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ref2.ps
2001.6187	7.1	7.9
2002.4894	36.3	17.0
2002.4921	-1.0	23.2
2002.4949	1.0	17.6
2002.4976	-0.3	21.0
2002.5003	6.0	16.5
2004.4689	-6.7	10.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> ref2.ps
SBreak	Segment0	Comp1
 2001.6187	   4.9527	Segment0	Comp1
 2004.4689	  -7.3963	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-240/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"REF2 Up (mm)":WeSn  \
-O >> ref2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> ref2.ps
2001.6187	7.6	23.6
2002.4894	-234.6	40.2
2002.4921	-2.7	59.9
2002.4949	3.7	42.0
2002.4976	-0.7	65.8
2002.5003	-12.5	44.0
2004.4689	1.7	29.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> ref2.ps
2001.6187	7.6	23.6
2002.4894	-234.6	40.2
2002.4921	-2.7	59.9
2002.4949	3.7	42.0
2002.4976	-0.7	65.8
2002.5003	-12.5	44.0
2004.4689	1.7	29.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> ref2.ps
SBreak	Segment0	Comp2
 2001.6187	   0.8297	Segment0	Comp2
 2004.4689	  -1.3463	Segment0	Comp2
EOT