#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2004.5/-20/30 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"CAS2 North (mm)":WeSn  \
> cas2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cas2.ps
1993.0547	-11.3	5.6
1994.3167	-18.9	5.8
1995.4858	-7.3	5.8
1996.4304	-9.4	5.9
1996.4852	0.6	5.9
1997.9418	-0.3	6.7
1997.9555	-5.6	6.6
1999.0343	-12.4	6.8
1999.0370	-4.9	5.2
2001.0850	15.7	3.3
2001.0877	17.1	6.1
2004.2445	29.1	5.2
2004.2472	27.9	5.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> cas2.ps
1993.0547	-11.3	5.6
1994.3167	-18.9	5.8
1995.4858	-7.3	5.8
1996.4304	-9.4	5.9
1996.4852	0.6	5.9
1997.9418	-0.3	6.7
1997.9555	-5.6	6.6
1999.0343	-12.4	6.8
1999.0370	-4.9	5.2
2001.0850	15.7	3.3
2001.0877	17.1	6.1
2004.2445	29.1	5.2
2004.2472	27.9	5.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cas2.ps
SBreak	Segment0	Comp0
 1993.0547	 -17.8990	Segment0	Comp0
 2004.2472	  21.9928	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-150/190 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"CAS2 East (mm)":WeSn  \
-O >> cas2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cas2.ps
1993.0547	181.4	7.3
1994.3167	127.2	6.9
1995.4858	89.6	7.1
1996.4304	64.2	7.6
1996.4852	70.9	7.6
1997.9418	33.4	9.0
1997.9555	13.5	9.9
1999.0343	-7.9	8.1
1999.0370	5.1	7.5
2001.0850	-56.3	9.0
2001.0877	-55.4	7.7
2004.2445	-144.3	8.0
2004.2472	-145.7	7.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> cas2.ps
1993.0547	181.4	7.3
1994.3167	127.2	6.9
1995.4858	89.6	7.1
1996.4304	64.2	7.6
1996.4852	70.9	7.6
1997.9418	33.4	9.0
1997.9555	13.5	9.9
1999.0343	-7.9	8.1
1999.0370	5.1	7.5
2001.0850	-56.3	9.0
2001.0877	-55.4	7.7
2004.2445	-144.3	8.0
2004.2472	-145.7	7.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cas2.ps
SBreak	Segment0	Comp1
 1993.0547	 164.5571	Segment0	Comp1
 2004.2472	-145.1193	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-50/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"CAS2 Up (mm)":WeSn  \
-O >> cas2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cas2.ps
1993.0547	38.2	22.8
1994.3167	12.3	20.4
1995.4858	-23.7	21.0
1996.4304	24.0	20.3
1996.4852	2.4	20.5
1997.9418	-8.0	24.0
1997.9555	-20.8	23.8
1999.0343	-12.4	22.5
1999.0370	-5.4	19.7
2001.0850	-40.2	16.0
2001.0877	-11.6	21.1
2004.2445	27.2	19.6
2004.2472	-2.0	21.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> cas2.ps
1993.0547	38.2	22.8
1994.3167	12.3	20.4
1995.4858	-23.7	21.0
1996.4304	24.0	20.3
1996.4852	2.4	20.5
1997.9418	-8.0	24.0
1997.9555	-20.8	23.8
1999.0343	-12.4	22.5
1999.0370	-5.4	19.7
2001.0850	-40.2	16.0
2001.0877	-11.6	21.1
2004.2445	27.2	19.6
2004.2472	-2.0	21.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> cas2.ps
SBreak	Segment0	Comp2
 1993.0547	   9.1566	Segment0	Comp2
 2004.2472	  -9.4895	Segment0	Comp2
EOT