#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-20/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"BROD 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.6160	9.6	8.4
2001.6187	1.4	14.4
2001.6215	2.5	20.6
2001.6242	7.7	15.4
2002.4757	-1.9	21.8
2002.4784	15.2	17.9
2002.4812	8.8	9.2
2004.4716	-13.0	6.7
2004.4744	-11.5	24.4
2004.4771	-11.0	42.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> brod.ps
2001.6160	9.6	8.4
2001.6187	1.4	14.4
2001.6215	2.5	20.6
2001.6242	7.7	15.4
2002.4757	-1.9	21.8
2002.4784	15.2	17.9
2002.4812	8.8	9.2
2004.4716	-13.0	6.7
2004.4744	-11.5	24.4
2004.4771	-11.0	42.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> brod.ps
SBreak	Segment0	Comp0
 2001.6160	   7.9438	Segment0	Comp0
 2004.4771	 -10.0765	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-20/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"BROD 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.6160	-2.3	7.9
2001.6187	0.3	17.0
2001.6215	-0.1	23.2
2001.6242	4.0	17.6
2002.4757	-14.2	21.0
2002.4784	-5.5	16.5
2002.4812	-9.5	10.9
2004.4716	7.9	7.1
2004.4744	1.4	23.2
2004.4771	4.0	37.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> brod.ps
2001.6160	-2.3	7.9
2001.6187	0.3	17.0
2001.6215	-0.1	23.2
2001.6242	4.0	17.6
2002.4757	-14.2	21.0
2002.4784	-5.5	16.5
2002.4812	-9.5	10.9
2004.4716	7.9	7.1
2004.4744	1.4	23.2
2004.4771	4.0	37.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> brod.ps
SBreak	Segment0	Comp1
 2001.6160	  -2.1932	Segment0	Comp1
 2004.4771	   2.7592	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-50/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"BROD 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.6160	-40.4	23.6
2001.6187	-0.7	40.2
2001.6215	-13.7	59.9
2001.6242	0.3	42.0
2002.4757	30.2	65.8
2002.4784	-1.9	44.0
2002.4812	-10.1	29.1
2004.4716	4.8	22.5
2004.4744	14.5	65.5
2004.4771	5.3	159.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> brod.ps
2001.6160	-40.4	23.6
2001.6187	-0.7	40.2
2001.6215	-13.7	59.9
2001.6242	0.3	42.0
2002.4757	30.2	65.8
2002.4784	-1.9	44.0
2002.4812	-10.1	29.1
2004.4716	4.8	22.5
2004.4744	14.5	65.5
2004.4771	5.3	159.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> brod.ps
SBreak	Segment0	Comp2
 2001.6160	  -6.0848	Segment0	Comp2
 2004.4771	   7.7722	Segment0	Comp2
EOT