#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-70/80 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"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	-65.5	4.3
1994.2948	-59.4	5.8
1995.7240	-33.3	6.2
1996.8330	-51.6	5.7
1997.8871	-28.7	6.0
1997.9254	-33.0	6.1
1998.9384	-23.0	5.8
1998.9412	-22.0	6.6
2001.0439	-1.3	6.7
2001.0466	6.5	6.4
2003.1001	25.1	5.8
2003.1110	25.3	5.6
2006.1063	54.5	5.2
2006.1255	44.6	5.1
2009.1918	77.0	5.3
2009.1945	73.6	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	-65.5	4.3
1994.2948	-59.4	5.8
1995.7240	-33.3	6.2
1996.8330	-51.6	5.7
1997.8871	-28.7	6.0
1997.9254	-33.0	6.1
1998.9384	-23.0	5.8
1998.9412	-22.0	6.6
2001.0439	-1.3	6.7
2001.0466	6.5	6.4
2003.1001	25.1	5.8
2003.1110	25.3	5.6
2006.1063	54.5	5.2
2006.1255	44.6	5.1
2009.1918	77.0	5.3
2009.1945	73.6	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	 -70.8143	Segment0	Comp0
 2009.1945	  76.0616	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-70/90 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"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	81.2	9.5
1994.2948	73.1	7.7
1995.7240	49.7	7.7
1996.8330	23.1	7.5
1997.8871	34.8	7.7
1997.9254	28.6	7.7
1998.9384	31.6	7.5
1998.9412	32.2	9.7
2001.0439	12.9	7.9
2001.0466	18.4	8.6
2003.1001	-9.0	7.3
2003.1110	8.7	8.5
2006.1063	-33.0	6.1
2006.1255	-34.5	6.1
2009.1918	-66.2	5.6
2009.1945	-65.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	81.2	9.5
1994.2948	73.1	7.7
1995.7240	49.7	7.7
1996.8330	23.1	7.5
1997.8871	34.8	7.7
1997.9254	28.6	7.7
1998.9384	31.6	7.5
1998.9412	32.2	9.7
2001.0439	12.9	7.9
2001.0466	18.4	8.6
2003.1001	-9.0	7.3
2003.1110	8.7	8.5
2006.1063	-33.0	6.1
2006.1255	-34.5	6.1
2009.1918	-66.2	5.6
2009.1945	-65.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	  77.5039	Segment0	Comp1
 2009.1945	 -61.8916	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.4	22.8
1994.2948	73.1	27.5
1995.7240	-27.9	21.7
1996.8330	24.4	21.0
1997.8871	-14.7	22.2
1997.9254	-9.1	20.7
1998.9384	-6.9	20.1
1998.9412	15.2	22.8
2001.0439	-14.4	23.7
2001.0466	-19.1	22.1
2003.1001	-12.6	20.6
2003.1110	55.2	22.0
2006.1063	1.1	18.8
2006.1255	-19.3	18.2
2009.1918	-18.9	17.5
2009.1945	-23.5	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.4	22.8
1994.2948	73.1	27.5
1995.7240	-27.9	21.7
1996.8330	24.4	21.0
1997.8871	-14.7	22.2
1997.9254	-9.1	20.7
1998.9384	-6.9	20.1
1998.9412	15.2	22.8
2001.0439	-14.4	23.7
2001.0466	-19.1	22.1
2003.1001	-12.6	20.6
2003.1110	55.2	22.0
2006.1063	1.1	18.8
2006.1255	-19.3	18.2
2009.1918	-18.9	17.5
2009.1945	-23.5	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.4614	Segment0	Comp2
 2009.1945	 -20.6842	Segment0	Comp2
EOT