#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1996/2010/-100/110 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"DEAL North (mm)":WeSn  \
> deal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> deal.ps
1996.8686	-95.4	6.2
1997.8816	-91.1	6.1
1997.9199	-92.2	6.0
1998.9193	-65.5	5.8
1998.9357	-66.2	5.8
2001.0275	-23.5	6.0
2001.0302	-28.0	6.2
2003.0919	12.2	5.6
2003.0946	8.7	5.8
2006.1036	62.4	5.3
2006.1227	60.4	5.1
2009.1727	104.6	5.4
2009.1891	107.3	5.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	-95.4	6.2
1997.8816	-91.1	6.1
1997.9199	-92.2	6.0
1998.9193	-65.5	5.8
1998.9357	-66.2	5.8
2001.0275	-23.5	6.0
2001.0302	-28.0	6.2
2003.0919	12.2	5.6
2003.0946	8.7	5.8
2006.1036	62.4	5.3
2006.1227	60.4	5.1
2009.1891	107.3	5.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> deal.ps
SBreak	Segment0	Comp0
 1996.8686	-100.7118	Segment0	Comp0
 2009.1891	 110.9144	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-90/90 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"DEAL East (mm)":WeSn  \
-O >> deal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> deal.ps
1996.8686	85.5	7.6
1997.8816	66.5	7.4
1997.9199	67.9	7.4
1998.9193	58.0	7.3
1998.9357	57.3	7.4
2001.0275	37.6	6.7
2001.0302	33.6	7.2
2003.0919	2.1	7.5
2003.0946	0.2	7.7
2006.1036	-37.2	6.1
2006.1227	-44.3	5.9
2009.1727	-87.6	5.7
2009.1891	-86.9	5.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	85.5	7.6
1997.8816	66.5	7.4
1997.9199	67.9	7.4
1998.9193	58.0	7.3
1998.9357	57.3	7.4
2001.0275	37.6	6.7
2001.0302	33.6	7.2
2003.0919	2.1	7.5
2003.0946	0.2	7.7
2006.1036	-37.2	6.1
2006.1227	-44.3	5.9
2009.1891	-86.9	5.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> deal.ps
SBreak	Segment0	Comp1
 1996.8686	  86.9831	Segment0	Comp1
 2009.1891	 -84.9766	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-40/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"DEAL Up (mm)":WeSn  \
-O >> deal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> deal.ps
1996.8686	-5.5	23.1
1997.8816	10.7	23.7
1997.9199	11.7	20.2
1998.9193	2.9	21.0
1998.9357	5.4	20.5
2001.0275	12.3	21.0
2001.0302	14.8	22.3
2003.0919	10.8	21.3
2003.0946	14.1	22.5
2006.1036	5.1	18.7
2006.1227	-13.5	18.1
2009.1727	-10.1	18.8
2009.1891	-34.3	17.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	-5.5	23.1
1997.8816	10.7	23.7
1997.9199	11.7	20.2
1998.9193	2.9	21.0
1998.9357	5.4	20.5
2001.0275	12.3	21.0
2001.0302	14.8	22.3
2003.0919	10.8	21.3
2003.0946	14.1	22.5
2006.1036	5.1	18.7
2006.1227	-13.5	18.1
2009.1891	-34.3	17.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> deal.ps
SBreak	Segment0	Comp2
 1996.8686	  13.9822	Segment0	Comp2
 2009.1891	 -14.7529	Segment0	Comp2
EOT