#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-40/40 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"CAML North (mm)":WeSn  \
> caml.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> caml.ps
1993.0355	39.5	4.3
1994.2948	28.6	5.8
1995.7240	35.5	6.2
1996.8330	2.3	5.7
1997.8871	10.9	6.0
1997.9254	6.0	6.1
1998.9384	2.4	5.8
1998.9412	3.4	6.6
2001.0439	-4.2	6.7
2001.0466	3.6	6.4
2003.1001	-5.5	5.8
2003.1110	-5.5	5.6
2006.1063	-16.6	5.2
2006.1255	-26.8	5.1
2009.1918	-35.8	5.3
2009.1945	-39.2	5.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> caml.ps
1993.0355	39.5	4.3
1994.2948	28.6	5.8
1995.7240	35.5	6.2
1996.8330	2.3	5.7
1997.8871	10.9	6.0
1997.9254	6.0	6.1
1998.9384	2.4	5.8
1998.9412	3.4	6.6
2001.0439	-4.2	6.7
2001.0466	3.6	6.4
2003.1001	-5.5	5.8
2003.1110	-5.5	5.6
2006.1063	-16.6	5.2
2006.1255	-26.8	5.1
2009.1918	-35.8	5.3
2009.1945	-39.2	5.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> caml.ps
SBreak	Segment0	Comp0
 1993.0355	  34.1863	Segment0	Comp0
 2009.1945	 -36.7195	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-160/190 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"CAML East (mm)":WeSn  \
-O >> caml.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> caml.ps
1993.0355	187.6	9.5
1994.2948	164.5	7.7
1995.7240	124.2	7.7
1996.8330	84.4	7.5
1997.8871	83.7	7.7
1997.9254	77.1	7.7
1998.9384	68.1	7.5
1998.9412	68.6	9.7
2001.0439	24.3	7.9
2001.0466	29.9	8.6
2003.1001	-21.8	7.3
2003.1110	-4.2	8.5
2006.1063	-81.3	6.1
2006.1255	-83.1	6.1
2009.1918	-151.1	5.6
2009.1945	-150.9	5.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> caml.ps
1993.0355	187.6	9.5
1994.2948	164.5	7.7
1995.7240	124.2	7.7
1996.8330	84.4	7.5
1997.8871	83.7	7.7
1997.9254	77.1	7.7
1998.9384	68.1	7.5
1998.9412	68.6	9.7
2001.0439	24.3	7.9
2001.0466	29.9	8.6
2003.1001	-21.8	7.3
2003.1110	-4.2	8.5
2006.1063	-81.3	6.1
2006.1255	-83.1	6.1
2009.1918	-151.1	5.6
2009.1945	-150.9	5.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> caml.ps
SBreak	Segment0	Comp1
 1993.0355	 183.8594	Segment0	Comp1
 2009.1945	-146.8229	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-30/80 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"CAML Up (mm)":WeSn  \
-O >> caml.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> caml.ps
1993.0355	58.8	22.8
1994.2948	73.5	27.5
1995.7240	-27.6	21.7
1996.8330	24.6	21.0
1997.8871	-14.5	22.2
1997.9254	-8.9	20.7
1998.9384	-6.7	20.1
1998.9412	15.2	22.8
2001.0439	-14.4	23.7
2001.0466	-19.0	22.1
2003.1001	-12.6	20.6
2003.1110	55.1	22.0
2006.1063	0.9	18.8
2006.1255	-19.5	18.2
2009.1918	-19.3	17.5
2009.1945	-23.8	18.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> caml.ps
1993.0355	58.8	22.8
1994.2948	73.5	27.5
1995.7240	-27.6	21.7
1996.8330	24.6	21.0
1997.8871	-14.5	22.2
1997.9254	-8.9	20.7
1998.9384	-6.7	20.1
1998.9412	15.2	22.8
2001.0439	-14.4	23.7
2001.0466	-19.0	22.1
2003.1001	-12.6	20.6
2003.1110	55.1	22.0
2006.1063	0.9	18.8
2006.1255	-19.5	18.2
2009.1918	-19.3	17.5
2009.1945	-23.8	18.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> caml.ps
SBreak	Segment0	Comp2
 1993.0355	  23.8817	Segment0	Comp2
 2009.1945	 -21.0548	Segment0	Comp2
EOT