#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-60/30 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"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	23.5	8.4
2002.4730	6.3	14.4
2002.4757	7.7	20.6
2002.4784	6.9	15.4
2002.4812	10.4	21.8
2002.4839	9.8	17.9
2002.4867	6.6	9.2
2002.4894	13.7	6.7
2002.4921	10.7	24.4
2002.4949	5.5	42.7
2002.4976	2.6	47.8
2004.4689	-54.8	28.9
2004.4716	-52.2	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	23.5	8.4
2002.4730	6.3	14.4
2002.4757	7.7	20.6
2002.4784	6.9	15.4
2002.4812	10.4	21.8
2002.4839	9.8	17.9
2002.4867	6.6	9.2
2002.4894	13.7	6.7
2002.4921	10.7	24.4
2002.4949	5.5	42.7
2002.4976	2.6	47.8
2004.4689	-54.8	28.9
2004.4716	-52.2	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	  32.3657	Segment0	Comp0
 2004.4716	 -51.6433	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	14.6	7.9
2002.4730	7.0	17.0
2002.4757	-2.9	23.2
2002.4784	-1.5	17.6
2002.4812	-3.1	21.0
2002.4839	2.6	16.5
2002.4867	-1.4	10.9
2002.4894	0.4	7.1
2002.4921	-1.1	23.2
2002.4949	-1.5	37.1
2002.4976	-0.7	49.7
2004.4689	-2.3	37.9
2004.4716	-7.4	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	14.6	7.9
2002.4730	7.0	17.0
2002.4757	-2.9	23.2
2002.4784	-1.5	17.6
2002.4812	-3.1	21.0
2002.4839	2.6	16.5
2002.4867	-1.4	10.9
2002.4894	0.4	7.1
2002.4921	-1.1	23.2
2002.4949	-1.5	37.1
2002.4976	-0.7	49.7
2004.4689	-2.3	37.9
2004.4716	-7.4	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	   4.3588	Segment0	Comp1
 2004.4716	  -6.9391	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.5	23.6
2002.4730	2.7	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.4	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	10.9	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.5	23.6
2002.4730	2.7	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.4	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	10.9	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.8283	Segment0	Comp2
 2004.4716	  14.0176	Segment0	Comp2
EOT