#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/10 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"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	-2.9	8.1
2002.4730	0.6	4.9
2002.4757	0.4	3.5
2002.4784	-1.8	3.7
2002.4812	3.3	3.7
2002.4839	2.8	3.6
2002.4867	-1.9	3.7
2002.4894	4.9	3.8
2002.4921	3.1	3.8
2002.4949	-3.3	3.8
2002.4976	-5.3	3.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	-2.9	8.1
2002.4730	0.6	4.9
2002.4757	0.4	3.5
2002.4784	-1.8	3.7
2002.4812	3.3	3.7
2002.4839	2.8	3.6
2002.4867	-1.9	3.7
2002.4894	4.9	3.8
2002.4921	3.1	3.8
2002.4949	-3.3	3.8
2002.4976	-5.3	3.8
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/20 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"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	10.1	9.0
2002.4730	6.8	5.1
2002.4757	-2.5	3.9
2002.4784	-2.0	4.2
2002.4812	-3.4	4.3
2002.4839	2.7	4.3
2002.4867	-1.1	4.5
2002.4894	0.1	4.6
2002.4921	-3.2	4.1
2002.4949	-2.1	4.2
2002.4976	-2.5	4.2
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	10.1	9.0
2002.4730	6.8	5.1
2002.4757	-2.5	3.9
2002.4784	-2.0	4.2
2002.4812	-3.4	4.3
2002.4839	2.7	4.3
2002.4867	-1.1	4.5
2002.4894	0.1	4.6
2002.4921	-3.2	4.1
2002.4949	-2.1	4.2
2002.4976	-2.5	4.2
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-20/20 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"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	-13.3	23.1
2002.4730	7.0	15.1
2002.4757	-4.0	11.2
2002.4784	6.0	12.0
2002.4812	-4.6	12.0
2002.4839	-2.5	11.3
2002.4867	14.4	12.0
2002.4894	7.2	12.7
2002.4921	1.7	12.3
2002.4949	-5.4	12.3
2002.4976	-7.7	12.7
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> lav2.ps
2001.6160	-13.3	23.1
2002.4730	7.0	15.1
2002.4757	-4.0	11.2
2002.4784	6.0	12.0
2002.4812	-4.6	12.0
2002.4839	-2.5	11.3
2002.4867	14.4	12.0
2002.4894	7.2	12.7
2002.4921	1.7	12.3
2002.4949	-5.4	12.3
2002.4976	-7.7	12.7
EOT


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