#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1999/2005.5/-50/40 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"P300 North (mm)":WeSn  \
> p300.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> p300.ps
1999.0562	-41.9	5.6
1999.0589	-35.2	5.8
2001.1041	-11.1	5.8
2001.1069	-1.6	5.9
2005.0522	30.4	5.9
2005.0549	30.2	6.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> p300.ps
1999.0562	-41.9	5.6
1999.0589	-35.2	5.8
2001.1041	-11.1	5.8
2001.1069	-1.6	5.9
2005.0522	30.4	5.9
2005.0549	30.2	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> p300.ps
SBreak	Segment0	Comp0
 1999.0562	 -34.6206	Segment0	Comp0
 2005.0549	  31.7529	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1999/2005.5/-40/50 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"P300 East (mm)":WeSn  \
-O >> p300.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> p300.ps
1999.0562	43.2	7.3
1999.0589	39.4	6.9
2001.1041	-1.8	7.1
2001.1069	9.3	7.6
2005.0522	-34.7	7.6
2005.0549	-34.0	9.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> p300.ps
1999.0562	43.2	7.3
1999.0589	39.4	6.9
2001.1041	-1.8	7.1
2001.1069	9.3	7.6
2005.0522	-34.7	7.6
2005.0549	-34.0	9.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> p300.ps
SBreak	Segment0	Comp1
 1999.0562	  36.4154	Segment0	Comp1
 2005.0549	 -36.4171	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1999/2005.5/-20/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"P300 Up (mm)":WeSn  \
-O >> p300.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> p300.ps
1999.0562	22.7	22.8
1999.0589	24.2	20.4
2001.1041	-5.8	21.0
2001.1069	-7.9	20.3
2005.0522	-13.2	20.5
2005.0549	-13.3	24.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> p300.ps
1999.0562	22.7	22.8
1999.0589	24.2	20.4
2001.1041	-5.8	21.0
2001.1069	-7.9	20.3
2005.0522	-13.2	20.5
2005.0549	-13.3	24.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> p300.ps
SBreak	Segment0	Comp2
 1999.0562	  15.8315	Segment0	Comp2
 2005.0549	 -16.4070	Segment0	Comp2
EOT