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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> brod.ps
2001.6167	-3.3	8.1
2001.6194	-5.3	4.9
2001.6222	1.1	3.5
2002.4764	6.7	3.7
2002.4791	1.5	3.7
EOT

#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> brod.ps
2001.6167	-3.3	8.1
2001.6194	-5.3	4.9
2001.6222	1.1	3.5
2002.4764	6.7	3.7
2002.4791	1.5	3.7
EOT

#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> brod.ps
SBreak	Segment0	Comp0
 2001.6167	  -2.6369	Segment0	Comp0
 2002.4791	   4.0931	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 >> brod.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> brod.ps
2001.6167	1.2	9.0
2001.6194	1.5	5.1
2001.6222	7.3	3.9
2002.4764	-3.7	4.2
2002.4791	-6.1	4.3
EOT

#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> brod.ps
2001.6167	1.2	9.0
2001.6194	1.5	5.1
2001.6222	7.3	3.9
2002.4764	-3.7	4.2
2002.4791	-6.1	4.3
EOT

#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> brod.ps
SBreak	Segment0	Comp1
 2001.6167	   3.5220	Segment0	Comp1
 2002.4791	  -4.8804	Segment0	Comp1
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:"Up (mm)":WeSn  \
-O >> brod.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> brod.ps
2001.6167	5.0	23.1
2001.6194	-7.6	15.1
2001.6222	7.6	11.2
2002.4764	1.6	12.0
2002.4791	-6.1	12.0
EOT

#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> brod.ps
2001.6167	5.0	23.1
2001.6194	-7.6	15.1
2001.6222	7.6	11.2
2002.4764	1.6	12.0
2002.4791	-6.1	12.0
EOT

#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> brod.ps
SBreak	Segment0	Comp2
 2001.6167	   1.4740	Segment0	Comp2
 2002.4791	  -2.2371	Segment0	Comp2
EOT