#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2004.5/-80/120 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-79.0	5.6
1994.3167	-69.6	5.8
1995.4858	-42.2	5.8
1996.4304	-31.7	5.9
1996.4852	-20.8	5.9
1997.9418	-2.2	6.7
1997.9555	-7.2	6.6
1999.0343	0.6	6.8
1999.0370	8.1	5.2
2001.0850	56.4	3.3
2001.0877	57.7	6.1
2004.2445	112.4	5.2
2004.2472	111.2	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	-79.0	5.6
1994.3167	-69.6	5.8
1995.4858	-42.2	5.8
1996.4304	-31.7	5.9
1996.4852	-20.8	5.9
1997.9418	-2.2	6.7
1997.9555	-7.2	6.6
1999.0343	0.6	6.8
1999.0370	8.1	5.2
2001.0850	56.4	3.3
2001.0877	57.7	6.1
2004.2445	112.4	5.2
2004.2472	111.2	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	 -85.6788	Segment0	Comp0
 2004.2472	 105.2745	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-90/120 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"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	111.8	7.3
1994.3167	72.5	6.9
1995.4858	48.5	7.1
1996.4304	34.1	7.6
1996.4852	41.5	7.6
1997.9418	21.0	9.0
1997.9555	1.3	9.9
1999.0343	-7.5	8.1
1999.0370	5.5	7.5
2001.0850	-32.0	9.0
2001.0877	-31.0	7.7
2004.2445	-83.0	8.0
2004.2472	-84.5	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	111.8	7.3
1994.3167	72.5	6.9
1995.4858	48.5	7.1
1996.4304	34.1	7.6
1996.4852	41.5	7.6
1997.9418	21.0	9.0
1997.9555	1.3	9.9
1999.0343	-7.5	8.1
1999.0370	5.5	7.5
2001.0850	-32.0	9.0
2001.0877	-31.0	7.7
2004.2445	-83.0	8.0
2004.2472	-84.5	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	  95.0405	Segment0	Comp1
 2004.2472	 -83.8141	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	37.9	22.8
1994.3167	12.2	20.4
1995.4858	-23.8	21.0
1996.4304	23.9	20.3
1996.4852	2.4	20.5
1997.9418	-8.0	24.0
1997.9555	-20.8	23.8
1999.0343	-12.3	22.5
1999.0370	-5.4	19.7
2001.0850	-40.1	16.0
2001.0877	-11.4	21.1
2004.2445	27.4	19.6
2004.2472	-1.7	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	37.9	22.8
1994.3167	12.2	20.4
1995.4858	-23.8	21.0
1996.4304	23.9	20.3
1996.4852	2.4	20.5
1997.9418	-8.0	24.0
1997.9555	-20.8	23.8
1999.0343	-12.3	22.5
1999.0370	-5.4	19.7
2001.0850	-40.1	16.0
2001.0877	-11.4	21.1
2004.2445	27.4	19.6
2004.2472	-1.7	21.3
EOT


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