#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2007/-50/40 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"AIRR North (mm)":WeSn  \
> airr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> airr.ps
1993.0355	-49.8	4.2
1993.0383	-47.4	4.0
1994.2976	-34.0	4.6
1995.7185	-13.7	5.9
1996.8275	-36.0	5.4
1997.8816	-15.7	5.5
1997.9199	-21.5	5.4
1998.9165	-1.7	6.2
1998.9357	-9.9	5.5
2001.1041	13.7	5.7
2001.1069	9.3	3.5
2003.0919	16.6	5.1
2003.0946	14.3	3.2
2006.0871	30.7	2.9
2006.0899	31.6	2.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> airr.ps
1993.0355	-49.8	4.2
1993.0383	-47.4	4.0
1994.2976	-34.0	4.6
1995.7185	-13.7	5.9
1996.8275	-36.0	5.4
1997.8816	-15.7	5.5
1997.9199	-21.5	5.4
1998.9165	-1.7	6.2
1998.9357	-9.9	5.5
2001.1041	13.7	5.7
2001.1069	9.3	3.5
2003.0919	16.6	5.1
2003.0946	14.3	3.2
2006.0871	30.7	2.9
2006.0899	31.6	2.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> airr.ps
SBreak	Segment0	Comp0
 1993.0355	 -44.8451	Segment0	Comp0
 2006.0899	  33.1013	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2007/-110/250 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"AIRR East (mm)":WeSn  \
-O >> airr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> airr.ps
1993.0355	248.7	9.6
1993.0383	229.9	11.1
1994.2976	198.3	6.7
1995.7185	180.5	7.8
1996.8275	129.7	7.0
1997.8816	114.0	6.8
1997.9199	110.6	6.9
1998.9165	89.4	10.4
1998.9357	79.5	6.5
2001.1041	17.4	9.0
2001.1069	21.8	4.7
2003.0919	-22.5	6.5
2003.0946	-22.8	4.3
2006.0871	-106.6	3.2
2006.0899	-109.7	3.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> airr.ps
1993.0355	248.7	9.6
1993.0383	229.9	11.1
1994.2976	198.3	6.7
1995.7185	180.5	7.8
1996.8275	129.7	7.0
1997.8816	114.0	6.8
1997.9199	110.6	6.9
1998.9165	89.4	10.4
1998.9357	79.5	6.5
2001.1041	17.4	9.0
2001.1069	21.8	4.7
2003.0919	-22.5	6.5
2003.0946	-22.8	4.3
2006.0871	-106.6	3.2
2006.0899	-109.7	3.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> airr.ps
SBreak	Segment0	Comp1
 1993.0355	 238.8753	Segment0	Comp1
 2006.0899	-107.3119	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2007/-30/120 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"AIRR Up (mm)":WeSn  \
-O >> airr.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> airr.ps
1993.0355	111.8	23.5
1993.0383	98.1	22.6
1994.2976	25.6	19.0
1995.7185	0.1	20.5
1996.8275	-7.1	20.9
1997.8816	-16.9	19.7
1997.9199	-11.9	19.2
1998.9165	14.2	22.1
1998.9357	-16.7	19.1
2001.1041	-23.7	21.9
2001.1069	-15.5	12.7
2003.0919	-10.2	19.1
2003.0946	9.5	11.5
2006.0871	-10.5	9.8
2006.0899	-14.5	10.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> airr.ps
1993.0355	111.8	23.5
1993.0383	98.1	22.6
1994.2976	25.6	19.0
1995.7185	0.1	20.5
1996.8275	-7.1	20.9
1997.8816	-16.9	19.7
1997.9199	-11.9	19.2
1998.9165	14.2	22.1
1998.9357	-16.7	19.1
2001.1041	-23.7	21.9
2001.1069	-15.5	12.7
2003.0919	-10.2	19.1
2003.0946	9.5	11.5
2006.0871	-10.5	9.8
2006.0899	-14.5	10.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> airr.ps
SBreak	Segment0	Comp2
 1993.0355	  34.5046	Segment0	Comp2
 2006.0899	 -18.8179	Segment0	Comp2
EOT