#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994/2010/-30/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"MADI North (mm)":WeSn  \
> madi.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> madi.ps
1994.2948	-8.6	4.3
1995.7213	-7.4	6.3
1996.8330	-13.1	5.8
1997.8843	-10.3	6.1
1997.9227	-22.8	6.9
1998.9193	-8.8	5.5
1998.9384	-12.6	5.8
2001.0439	1.4	6.3
2001.2383	-2.9	5.9
2003.0973	2.0	5.5
2003.1165	2.9	5.6
2006.1446	5.3	5.1
2006.1474	8.0	5.1
2009.1891	7.8	2.8
2009.1918	6.2	2.8
2009.1945	3.5	6.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> madi.ps
1994.2948	-8.6	4.3
1995.7213	-7.4	6.3
1996.8330	-13.1	5.8
1997.8843	-10.3	6.1
1997.9227	-22.8	6.9
1998.9193	-8.8	5.5
1998.9384	-12.6	5.8
2001.0439	1.4	6.3
2001.2383	-2.9	5.9
2003.0973	2.0	5.5
2003.1165	2.9	5.6
2006.1446	5.3	5.1
2006.1474	8.0	5.1
2009.1891	7.8	2.8
2009.1918	6.2	2.8
2009.1945	3.5	6.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> madi.ps
SBreak	Segment0	Comp0
 1994.2948	 -12.9718	Segment0	Comp0
 2009.1945	   7.4663	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994/2010/-120/250 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"MADI East (mm)":WeSn  \
-O >> madi.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> madi.ps
1994.2948	243.4	6.6
1995.7213	215.6	8.0
1996.8330	185.4	7.6
1997.8843	162.7	7.6
1997.9227	166.7	8.8
1998.9193	142.4	7.1
1998.9384	141.4	7.6
2001.0439	84.7	9.3
2001.2383	90.7	19.0
2003.0973	42.9	7.3
2003.1165	45.5	7.8
2006.1446	-32.6	5.8
2006.1474	-36.4	6.0
2009.1891	-113.8	3.3
2009.1918	-114.5	3.4
2009.1945	-114.4	6.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> madi.ps
1994.2948	243.4	6.6
1995.7213	215.6	8.0
1996.8330	185.4	7.6
1997.8843	162.7	7.6
1997.9227	166.7	8.8
1998.9193	142.4	7.1
1998.9384	141.4	7.6
2001.0439	84.7	9.3
2001.2383	90.7	19.0
2003.0973	42.9	7.3
2003.1165	45.5	7.8
2006.1446	-32.6	5.8
2006.1474	-36.4	6.0
2009.1891	-113.8	3.3
2009.1918	-114.5	3.4
2009.1945	-114.4	6.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> madi.ps
SBreak	Segment0	Comp1
 1994.2948	 251.6463	Segment0	Comp1
 2009.1945	-112.5002	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1994/2010/-30/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"MADI Up (mm)":WeSn  \
-O >> madi.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> madi.ps
1994.2948	36.2	18.5
1995.7213	18.2	22.2
1996.8330	25.3	21.2
1997.8843	-1.5	22.0
1997.9227	-3.8	25.6
1998.9193	-14.6	19.9
1998.9384	-4.5	20.1
2001.0439	-21.6	24.6
2001.2383	15.0	25.4
2003.0973	4.1	21.6
2003.1165	23.1	21.3
2006.1446	2.1	17.7
2006.1474	-1.9	18.4
2009.1891	-4.5	10.2
2009.1918	-11.7	10.2
2009.1945	-13.8	20.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> madi.ps
1994.2948	36.2	18.5
1995.7213	18.2	22.2
1996.8330	25.3	21.2
1997.8843	-1.5	22.0
1997.9227	-3.8	25.6
1998.9193	-14.6	19.9
1998.9384	-4.5	20.1
2001.0439	-21.6	24.6
2001.2383	15.0	25.4
2003.0973	4.1	21.6
2003.1165	23.1	21.3
2006.1446	2.1	17.7
2006.1474	-1.9	18.4
2009.1891	-4.5	10.2
2009.1918	-11.7	10.2
2009.1945	-13.8	20.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> madi.ps
SBreak	Segment0	Comp2
 1994.2948	  16.6310	Segment0	Comp2
 2009.1945	  -8.6406	Segment0	Comp2
EOT