#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1992.5/2004/-20/10 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"North (mm)":WeSn  \
> pibu.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pibu.ps
1992.6713	-17.7	7.1
1998.6344	1.6	6.6
1998.6673	-5.7	7.2
1998.6700	-6.5	6.8
1998.6728	-1.5	7.3
1999.1328	7.6	7.0
1999.1356	-5.1	6.8
2002.6345	5.6	6.5
2002.6564	7.6	6.3
2003.5599	7.8	6.4
2003.5709	5.6	6.6
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pibu.ps
1992.6713	-17.7	7.1
1998.6344	1.6	6.6
1998.6673	-5.7	7.2
1998.6700	-6.5	6.8
1998.6728	-1.5	7.3
1999.1328	7.6	7.0
1999.1356	-5.1	6.8
2002.6345	5.6	6.5
2002.6564	7.6	6.3
2003.5599	7.8	6.4
2003.5709	5.6	6.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> pibu.ps
SBreak	Segment0	Comp0
 1992.6713	 -16.1765	Segment0	Comp0
 2003.5709	   8.4177	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1992.5/2004/-10/10 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"East (mm)":WeSn  \
-O >> pibu.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pibu.ps
1992.6713	-2.1	8.0
1998.6344	1.5	7.2
1998.6673	8.5	7.8
1998.6700	-2.9	7.3
1998.6728	-2.3	9.8
1999.1328	-6.5	8.7
1999.1356	2.6	8.0
2002.6345	-2.2	7.5
2002.6564	8.8	7.4
2003.5599	-7.7	7.4
2003.5709	3.2	7.9
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pibu.ps
1992.6713	-2.1	8.0
1998.6344	1.5	7.2
1998.6673	8.5	7.8
1998.6700	-2.9	7.3
1998.6728	-2.3	9.8
1999.1328	-6.5	8.7
1999.1356	2.6	8.0
2002.6345	-2.2	7.5
2002.6564	8.8	7.4
2003.5599	-7.7	7.4
2003.5709	3.2	7.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> pibu.ps
SBreak	Segment0	Comp1
 1992.6713	  -0.0262	Segment0	Comp1
 2003.5709	   0.0155	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1992.5/2004/-50/40 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"Up (mm)":WeSn  \
-O >> pibu.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> pibu.ps
1992.6713	6.3	22.9
1998.6344	-9.7	21.9
1998.6673	13.8	23.5
1998.6700	-5.2	21.5
1998.6728	-15.9	23.5
1999.1328	26.5	26.4
1999.1356	-36.6	23.6
2002.6345	-41.6	22.3
2002.6564	0.5	22.8
2003.5599	32.3	22.6
2003.5709	35.2	24.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> pibu.ps
1992.6713	6.3	22.9
1998.6344	-9.7	21.9
1998.6673	13.8	23.5
1998.6700	-5.2	21.5
1998.6728	-15.9	23.5
1999.1328	26.5	26.4
1999.1356	-36.6	23.6
2002.6345	-41.6	22.3
2002.6564	0.5	22.8
2003.5599	32.3	22.6
2003.5709	35.2	24.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> pibu.ps
SBreak	Segment0	Comp2
 1992.6713	  -6.9347	Segment0	Comp2
 2003.5709	   3.5550	Segment0	Comp2
EOT