#!/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  \
> open.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> open.ps
2001.6187	-2.5	8.1
2001.6215	-0.2	4.9
2001.6242	2.9	3.5
2002.4675	-0.2	3.7
2002.4702	-2.1	3.7
2002.4730	0.3	3.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> open.ps
2001.6187	-2.5	8.1
2001.6215	-0.2	4.9
2001.6242	2.9	3.5
2002.4675	-0.2	3.7
2002.4702	-2.1	3.7
2002.4730	0.3	3.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> open.ps
SBreak	Segment0	Comp0
 2001.6187	   1.0057	Segment0	Comp0
 2002.4730	  -1.1818	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/10 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"East (mm)":WeSn  \
-O >> open.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> open.ps
2001.6187	-4.1	9.0
2001.6215	-2.8	5.1
2001.6242	-1.8	3.9
2002.4675	4.0	4.2
2002.4702	3.0	4.3
2002.4730	-4.9	4.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> open.ps
2001.6187	-4.1	9.0
2001.6215	-2.8	5.1
2001.6242	-1.8	3.9
2002.4675	4.0	4.2
2002.4702	3.0	4.3
2002.4730	-4.9	4.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> open.ps
SBreak	Segment0	Comp1
 2001.6187	  -2.5067	Segment0	Comp1
 2002.4730	   3.1334	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/30 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"Up (mm)":WeSn  \
-O >> open.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> open.ps
2001.6187	27.9	23.1
2001.6215	-9.8	15.1
2001.6242	-3.2	11.2
2002.4675	10.5	12.0
2002.4702	-3.4	12.0
2002.4730	18.3	11.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> open.ps
2001.6187	27.9	23.1
2001.6215	-9.8	15.1
2001.6242	-3.2	11.2
2002.4675	10.5	12.0
2002.4702	-3.4	12.0
2002.4730	18.3	11.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> open.ps
SBreak	Segment0	Comp2
 2001.6187	  -3.0392	Segment0	Comp2
 2002.4730	   3.5281	Segment0	Comp2
EOT