#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> lav2.ps
2001.6160	0.1	8.4
2002.4730	1.2	14.4
2002.4757	2.5	20.6
2002.4784	1.8	15.4
2002.4812	5.4	21.8
2002.4839	4.9	17.9
2002.4867	1.8	9.2
2002.4894	8.9	6.7
2002.4921	6.0	24.4
2002.4949	0.8	42.7
2002.4976	-2.1	47.8
2004.4689	-17.5	28.9
2004.4716	-14.8	16.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	0.1	8.4
2002.4730	1.2	14.4
2002.4757	2.5	20.6
2002.4784	1.8	15.4
2002.4812	5.4	21.8
2002.4839	4.9	17.9
2002.4867	1.8	9.2
2002.4894	8.9	6.7
2002.4921	6.0	24.4
2002.4949	0.8	42.7
2002.4976	-2.1	47.8
2004.4689	-17.5	28.9
2004.4716	-14.8	16.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> lav2.ps
SBreak	Segment0	Comp0
 2001.6160	   8.9317	Segment0	Comp0
 2004.4716	 -14.2516	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"LAV2 East (mm)":WeSn  \
-O >> lav2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> lav2.ps
2001.6160	11.1	7.9
2002.4730	6.2	17.0
2002.4757	-3.7	23.2
2002.4784	-2.3	17.6
2002.4812	-3.9	21.0
2002.4839	1.8	16.5
2002.4867	-2.1	10.9
2002.4894	-0.4	7.1
2002.4921	-1.9	23.2
2002.4949	-2.3	37.1
2002.4976	-1.4	49.7
2004.4689	3.4	37.9
2004.4716	-1.7	17.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	11.1	7.9
2002.4730	6.2	17.0
2002.4757	-3.7	23.2
2002.4784	-2.3	17.6
2002.4812	-3.9	21.0
2002.4839	1.8	16.5
2002.4867	-2.1	10.9
2002.4894	-0.4	7.1
2002.4921	-1.9	23.2
2002.4949	-2.3	37.1
2002.4976	-1.4	49.7
2004.4689	3.4	37.9
2004.4716	-1.7	17.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> lav2.ps
SBreak	Segment0	Comp1
 2001.6160	   0.7773	Segment0	Comp1
 2004.4716	  -1.2375	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-30/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"LAV2 Up (mm)":WeSn  \
-O >> lav2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> lav2.ps
2001.6160	-20.6	23.6
2002.4730	2.6	40.2
2002.4757	-6.2	59.9
2002.4784	4.8	42.0
2002.4812	-5.3	65.8
2002.4839	0.1	44.0
2002.4867	13.5	29.1
2002.4894	3.7	22.5
2002.4921	-0.2	65.5
2002.4949	-9.5	159.6
2002.4976	-7.6	126.0
2004.4689	12.4	78.2
2004.4716	11.0	56.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	-20.6	23.6
2002.4730	2.6	40.2
2002.4757	-6.2	59.9
2002.4784	4.8	42.0
2002.4812	-5.3	65.8
2002.4839	0.1	44.0
2002.4867	13.5	29.1
2002.4894	3.7	22.5
2002.4921	-0.2	65.5
2002.4949	-9.5	159.6
2002.4976	-7.6	126.0
2004.4689	12.4	78.2
2004.4716	11.0	56.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> lav2.ps
SBreak	Segment0	Comp2
 2001.6160	  -8.8906	Segment0	Comp2
 2004.4716	  14.1164	Segment0	Comp2
EOT