#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"OPEN 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	-3.9	8.4
2001.6215	-0.1	14.4
2001.6242	4.4	20.6
2002.4675	0.1	15.4
2002.4702	-3.5	21.8
2002.4730	-1.6	17.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> open.ps
2001.6187	-3.9	8.4
2001.6215	-0.1	14.4
2001.6242	4.4	20.6
2002.4675	0.1	15.4
2002.4702	-3.5	21.8
2002.4730	-1.6	17.9
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"OPEN 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.7	7.9
2001.6215	-1.8	17.0
2001.6242	-1.7	23.2
2002.4675	3.1	17.6
2002.4702	2.7	21.0
2002.4730	-5.0	16.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> open.ps
2001.6187	-4.7	7.9
2001.6215	-1.8	17.0
2001.6242	-1.7	23.2
2002.4675	3.1	17.6
2002.4702	2.7	21.0
2002.4730	-5.0	16.5
EOT


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


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"OPEN 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.1	23.6
2001.6215	-8.4	40.2
2001.6242	-4.0	59.9
2002.4675	10.0	42.0
2002.4702	-3.3	65.8
2002.4730	16.9	44.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> open.ps
2001.6187	27.1	23.6
2001.6215	-8.4	40.2
2001.6242	-4.0	59.9
2002.4675	10.0	42.0
2002.4702	-3.3	65.8
2002.4730	16.9	44.0
EOT


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