#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994/2010/-140/90 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-136.5	4.3
1995.7213	-116.1	6.3
1996.8330	-106.8	5.8
1997.8843	-89.8	6.1
1997.9227	-101.8	6.9
1998.9193	-74.2	5.5
1998.9384	-77.7	5.8
2001.0439	-35.2	6.3
2001.2383	-37.0	5.9
2003.0973	-6.9	5.5
2003.1165	-5.6	5.6
2006.1446	37.7	5.1
2006.1474	40.4	5.1
2009.1891	81.4	2.8
2009.1918	79.9	2.8
2009.1945	77.2	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	-136.5	4.3
1995.7213	-116.1	6.3
1996.8330	-106.8	5.8
1997.8843	-89.8	6.1
1997.9227	-101.8	6.9
1998.9193	-74.2	5.5
1998.9384	-77.7	5.8
2001.0439	-35.2	6.3
2001.2383	-37.0	5.9
2003.0973	-6.9	5.5
2003.1165	-5.6	5.6
2006.1446	37.7	5.1
2006.1474	40.4	5.1
2009.1891	81.4	2.8
2009.1918	79.9	2.8
2009.1945	77.2	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	-140.9705	Segment0	Comp0
 2009.1945	  81.1399	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994/2010/-70/130 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"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	122.0	6.6
1995.7213	111.0	8.0
1996.8330	93.9	7.6
1997.8843	83.6	7.6
1997.9227	88.0	8.8
1998.9193	75.6	7.1
1998.9384	74.9	7.6
2001.0439	42.9	9.3
2001.2383	51.3	19.0
2003.0973	25.4	7.3
2003.1165	28.2	7.8
2006.1446	-14.3	5.8
2006.1474	-18.1	6.0
2009.1891	-59.6	3.3
2009.1918	-60.3	3.4
2009.1945	-60.1	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	122.0	6.6
1995.7213	111.0	8.0
1996.8330	93.9	7.6
1997.8843	83.6	7.6
1997.9227	88.0	8.8
1998.9193	75.6	7.1
1998.9384	74.9	7.6
2001.0439	42.9	9.3
2001.2383	51.3	19.0
2003.0973	25.4	7.3
2003.1165	28.2	7.8
2006.1446	-14.3	5.8
2006.1474	-18.1	6.0
2009.1891	-59.6	3.3
2009.1918	-60.3	3.4
2009.1945	-60.1	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	 130.2568	Segment0	Comp1
 2009.1945	 -58.2322	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	35.8	18.5
1995.7213	17.8	22.2
1996.8330	25.0	21.2
1997.8843	-1.8	22.0
1997.9227	-4.1	25.6
1998.9193	-14.8	19.9
1998.9384	-4.8	20.1
2001.0439	-21.7	24.6
2001.2383	14.8	25.4
2003.0973	4.1	21.6
2003.1165	23.0	21.3
2006.1446	2.2	17.7
2006.1474	-1.7	18.4
2009.1891	-4.3	10.2
2009.1918	-11.4	10.2
2009.1945	-13.6	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	35.8	18.5
1995.7213	17.8	22.2
1996.8330	25.0	21.2
1997.8843	-1.8	22.0
1997.9227	-4.1	25.6
1998.9193	-14.8	19.9
1998.9384	-4.8	20.1
2001.0439	-21.7	24.6
2001.2383	14.8	25.4
2003.0973	4.1	21.6
2003.1165	23.0	21.3
2006.1446	2.2	17.7
2006.1474	-1.7	18.4
2009.1891	-4.3	10.2
2009.1918	-11.4	10.2
2009.1945	-13.6	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.2062	Segment0	Comp2
 2009.1945	  -8.4199	Segment0	Comp2
EOT