#!/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:"FLOW North (mm)":WeSn  \
> flow.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> flow.ps
2001.6215	20.8	8.4
2001.6242	24.2	14.4
2002.4730	2.5	20.6
2002.4757	1.4	15.4
2004.4826	-51.6	21.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> flow.ps
2001.6215	20.8	8.4
2001.6242	24.2	14.4
2002.4730	2.5	20.6
2002.4757	1.4	15.4
2004.4826	-51.6	21.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> flow.ps
SBreak	Segment0	Comp0
 2001.6215	  23.0397	Segment0	Comp0
 2004.4826	 -51.0374	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"FLOW East (mm)":WeSn  \
-O >> flow.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> flow.ps
2001.6215	4.9	7.9
2001.6242	5.1	17.0
2002.4730	1.8	23.2
2002.4757	-4.4	17.6
2004.4826	-8.4	21.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> flow.ps
2001.6215	4.9	7.9
2001.6242	5.1	17.0
2002.4730	1.8	23.2
2002.4757	-4.4	17.6
2004.4826	-8.4	21.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> flow.ps
SBreak	Segment0	Comp1
 2001.6215	   4.1176	Segment0	Comp1
 2004.4826	  -9.2588	Segment0	Comp1
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:"FLOW Up (mm)":WeSn  \
-O >> flow.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> flow.ps
2001.6215	-1.0	23.6
2001.6242	-2.1	40.2
2002.4730	11.5	59.9
2002.4757	-10.0	42.0
2004.4826	3.3	65.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> flow.ps
2001.6215	-1.0	23.6
2001.6242	-2.1	40.2
2002.4730	11.5	59.9
2002.4757	-10.0	42.0
2004.4826	3.3	65.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> flow.ps
SBreak	Segment0	Comp2
 2001.6215	  -1.4650	Segment0	Comp2
 2004.4826	   3.3257	Segment0	Comp2
EOT