#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1996/2010/-30/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"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	-15.9	6.2
1997.8816	-25.3	6.1
1997.9199	-26.9	6.0
1998.9193	-13.6	5.8
1998.9357	-14.7	5.8
2001.0275	-0.3	6.0
2001.0302	-4.9	6.2
2003.0919	7.3	5.6
2003.0946	3.8	5.8
2006.1036	16.6	5.3
2006.1227	14.3	5.1
2009.1727	17.2	5.4
2009.1891	19.6	5.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	-15.9	6.2
1997.8816	-25.3	6.1
1997.9199	-26.9	6.0
1998.9193	-13.6	5.8
1998.9357	-14.7	5.8
2001.0275	-0.3	6.0
2001.0302	-4.9	6.2
2003.0919	7.3	5.6
2003.0946	3.8	5.8
2006.1036	16.6	5.3
2006.1227	14.3	5.1
2009.1891	19.6	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	 -21.0935	Segment0	Comp0
 2009.1891	  23.2304	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2010/-160/160 \
-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	158.8	7.6
1997.8816	127.9	7.4
1997.9199	128.9	7.4
1998.9193	107.2	7.3
1998.9357	106.3	7.4
2001.0275	62.0	6.7
2001.0302	57.9	7.2
2003.0919	2.2	7.5
2003.0946	0.3	7.7
2006.1036	-72.6	6.1
2006.1227	-79.9	5.9
2009.1727	-159.0	5.7
2009.1891	-158.5	5.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	158.8	7.6
1997.8816	127.9	7.4
1997.9199	128.9	7.4
1998.9193	107.2	7.3
1998.9357	106.3	7.4
2001.0275	62.0	6.7
2001.0302	57.9	7.2
2003.0919	2.2	7.5
2003.0946	0.3	7.7
2006.1036	-72.6	6.1
2006.1227	-79.9	5.9
2009.1891	-158.5	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	 160.2991	Segment0	Comp1
 2009.1891	-156.6014	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.2	23.1
1997.8816	10.9	23.7
1997.9199	11.9	20.2
1998.9193	3.1	21.0
1998.9357	5.5	20.5
2001.0275	12.3	21.0
2001.0302	14.9	22.3
2003.0919	10.7	21.3
2003.0946	14.1	22.5
2006.1036	4.9	18.7
2006.1227	-13.7	18.1
2009.1727	-10.3	18.8
2009.1891	-34.6	17.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> deal.ps
1996.8686	-5.2	23.1
1997.8816	10.9	23.7
1997.9199	11.9	20.2
1998.9193	3.1	21.0
1998.9357	5.5	20.5
2001.0275	12.3	21.0
2001.0302	14.9	22.3
2003.0919	10.7	21.3
2003.0946	14.1	22.5
2006.1036	4.9	18.7
2006.1227	-13.7	18.1
2009.1891	-34.6	17.6
EOT


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