#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> jobe.ps
1993.0547	-45.2	5.6
1994.3167	-38.6	5.8
1995.4886	-20.1	5.8
1996.4441	-22.8	5.9
1997.9610	-12.7	5.9
1997.9637	-6.5	6.7
1999.0562	-2.5	6.6
1999.0589	0.5	6.8
2001.1041	31.2	5.2
2001.1069	25.1	3.3
2004.2608	49.8	6.1
2004.2636	46.3	5.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> jobe.ps
1993.0547	-45.2	5.6
1994.3167	-38.6	5.8
1995.4886	-20.1	5.8
1996.4441	-22.8	5.9
1997.9610	-12.7	5.9
1997.9637	-6.5	6.7
1999.0562	-2.5	6.6
1999.0589	0.5	6.8
2001.1041	31.2	5.2
2001.1069	25.1	3.3
2004.2608	49.8	6.1
2004.2636	46.3	5.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> jobe.ps
SBreak	Segment0	Comp0
 1993.0547	 -46.6621	Segment0	Comp0
 2004.2636	  49.0808	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-60/90 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"JOBE East (mm)":WeSn  \
-O >> jobe.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> jobe.ps
1993.0547	88.7	7.3
1994.3167	56.0	6.9
1995.4886	39.9	7.1
1996.4441	40.1	7.6
1997.9610	13.3	7.6
1997.9637	19.7	9.0
1999.0562	13.6	9.9
1999.0589	11.8	8.1
2001.1041	-23.1	7.5
2001.1069	-21.1	9.0
2004.2608	-56.9	7.7
2004.2636	-53.3	8.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> jobe.ps
1993.0547	88.7	7.3
1994.3167	56.0	6.9
1995.4886	39.9	7.1
1996.4441	40.1	7.6
1997.9610	13.3	7.6
1997.9637	19.7	9.0
1999.0562	13.6	9.9
1999.0589	11.8	8.1
2001.1041	-23.1	7.5
2001.1069	-21.1	9.0
2004.2608	-56.9	7.7
2004.2636	-53.3	8.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> jobe.ps
SBreak	Segment0	Comp1
 1993.0547	  77.1845	Segment0	Comp1
 2004.2636	 -56.4346	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2004.5/-40/110 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"JOBE Up (mm)":WeSn  \
-O >> jobe.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> jobe.ps
1993.0547	50.7	22.8
1994.3167	-2.2	20.4
1995.4886	10.7	21.0
1996.4441	103.9	20.3
1997.9610	-34.6	20.5
1997.9637	-10.6	24.0
1999.0562	-9.1	23.8
1999.0589	0.3	22.5
2001.1041	-16.6	19.7
2001.1069	-12.5	16.0
2004.2608	-8.2	21.1
2004.2636	-21.3	19.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> jobe.ps
1993.0547	50.7	22.8
1994.3167	-2.2	20.4
1995.4886	10.7	21.0
1996.4441	103.9	20.3
1997.9610	-34.6	20.5
1997.9637	-10.6	24.0
1999.0562	-9.1	23.8
1999.0589	0.3	22.5
2001.1041	-16.6	19.7
2001.1069	-12.5	16.0
2004.2608	-8.2	21.1
2004.2636	-21.3	19.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> jobe.ps
SBreak	Segment0	Comp2
 1993.0547	  30.9488	Segment0	Comp2
 2004.2636	 -25.8318	Segment0	Comp2
EOT