#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-20/210 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"MRGN North (mm)":WeSn  \
> mrgn.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mrgn.ps
2001.6187	200.4	8.4
2001.6215	5.9	14.4
2001.6242	8.0	20.6
2002.4702	0.3	15.4
2002.4730	0.4	21.8
2002.4757	2.2	17.9
2004.4607	-9.5	9.2
2004.4634	-10.8	6.7
2004.4662	-12.1	24.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> mrgn.ps
2001.6187	200.4	8.4
2001.6215	5.9	14.4
2001.6242	8.0	20.6
2002.4702	0.3	15.4
2002.4730	0.4	21.8
2002.4757	2.2	17.9
2004.4607	-9.5	9.2
2004.4634	-10.8	6.7
2004.4662	-12.1	24.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> mrgn.ps
SBreak	Segment0	Comp0
 2001.6187	   8.4677	Segment0	Comp0
 2004.4662	 -11.0563	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mrgn.ps
2001.6187	-436.4	7.9
2001.6215	2.3	17.0
2001.6242	2.7	23.2
2002.4702	5.5	17.6
2002.4730	3.8	21.0
2002.4757	-1.4	16.5
2004.4607	-2.2	10.9
2004.4634	-2.6	7.1
2004.4662	-4.9	23.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> mrgn.ps
2001.6187	-436.4	7.9
2001.6215	2.3	17.0
2001.6242	2.7	23.2
2002.4702	5.5	17.6
2002.4730	3.8	21.0
2002.4757	-1.4	16.5
2004.4607	-2.2	10.9
2004.4634	-2.6	7.1
2004.4662	-4.9	23.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> mrgn.ps
SBreak	Segment0	Comp1
 2001.6187	   1.8970	Segment0	Comp1
 2004.4662	  -2.1464	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-290/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"MRGN Up (mm)":WeSn  \
-O >> mrgn.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> mrgn.ps
2001.6187	-286.1	23.6
2001.6215	1.4	40.2
2001.6242	-6.2	59.9
2002.4702	6.6	42.0
2002.4730	6.8	65.8
2002.4757	-14.6	44.0
2004.4607	0.6	29.1
2004.4634	15.1	22.5
2004.4662	14.2	65.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> mrgn.ps
2001.6187	-286.1	23.6
2001.6215	1.4	40.2
2001.6242	-6.2	59.9
2002.4702	6.6	42.0
2002.4730	6.8	65.8
2002.4757	-14.6	44.0
2004.4607	0.6	29.1
2004.4634	15.1	22.5
2004.4662	14.2	65.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> mrgn.ps
SBreak	Segment0	Comp2
 2001.6187	  -7.7017	Segment0	Comp2
 2004.4662	  10.0524	Segment0	Comp2
EOT