#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1996/2004.5/-80/100 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"MSLL North (mm)":WeSn  \
> msll.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> msll.ps
1996.4852	-77.1	5.6
1997.9391	-64.3	5.8
1997.9665	-51.9	5.8
1999.0343	-41.4	5.9
1999.0370	-33.2	5.9
2001.0850	28.1	6.7
2001.0987	17.0	6.6
2004.1706	96.1	6.8
2004.1733	90.9	5.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> msll.ps
1996.4852	-77.1	5.6
1997.9391	-64.3	5.8
1997.9665	-51.9	5.8
1999.0343	-41.4	5.9
1999.0370	-33.2	5.9
2001.0850	28.1	6.7
2001.0987	17.0	6.6
2004.1706	96.1	6.8
2004.1733	90.9	5.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> msll.ps
SBreak	Segment0	Comp0
 1996.4852	 -89.4888	Segment0	Comp0
 2004.1733	  92.2268	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2004.5/-80/80 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"MSLL East (mm)":WeSn  \
-O >> msll.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> msll.ps
1996.4852	70.4	7.3
1997.9391	48.8	6.9
1997.9665	35.6	7.1
1999.0343	10.8	7.6
1999.0370	29.5	7.6
2001.0850	-18.5	9.0
2001.0987	-12.4	9.9
2004.1706	-71.7	8.1
2004.1733	-67.1	7.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> msll.ps
1996.4852	70.4	7.3
1997.9391	48.8	6.9
1997.9665	35.6	7.1
1999.0343	10.8	7.6
1999.0370	29.5	7.6
2001.0850	-18.5	9.0
2001.0987	-12.4	9.9
2004.1706	-71.7	8.1
2004.1733	-67.1	7.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> msll.ps
SBreak	Segment0	Comp1
 1996.4852	  69.0032	Segment0	Comp1
 2004.1733	 -69.8631	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1996/2004.5/-30/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"MSLL Up (mm)":WeSn  \
-O >> msll.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> msll.ps
1996.4852	26.1	22.8
1997.9391	-1.5	20.4
1997.9665	-15.4	21.0
1999.0343	-6.2	20.3
1999.0370	-2.2	20.5
2001.0850	-25.3	24.0
2001.0987	-7.6	23.8
2004.1706	27.1	22.5
2004.1733	6.9	19.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> msll.ps
1996.4852	26.1	22.8
1997.9391	-1.5	20.4
1997.9665	-15.4	21.0
1999.0343	-6.2	20.3
1999.0370	-2.2	20.5
2001.0850	-25.3	24.0
2001.0987	-7.6	23.8
2004.1706	27.1	22.5
2004.1733	6.9	19.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> msll.ps
SBreak	Segment0	Comp2
 1996.4852	  -6.7870	Segment0	Comp2
 2004.1733	   7.1989	Segment0	Comp2
EOT