#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993.5/2001/-30/50 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f10a20:"North (mm)":WeSn  \
> nhd4.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nhd4.ps
1993.6268	-22.7	7.1
1993.6295	-19.4	6.6
1993.6487	-26.5	7.2
1993.6514	-24.9	6.8
1994.6289	-18.8	7.3
1994.6316	-12.5	7.0
2000.5510	44.9	6.8
2000.5537	45.1	6.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> nhd4.ps
1993.6268	-22.7	7.1
1993.6295	-19.4	6.6
1993.6487	-26.5	7.2
1993.6514	-24.9	6.8
1994.6289	-18.8	7.3
1994.6316	-12.5	7.0
2000.5510	44.9	6.8
2000.5537	45.1	6.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> nhd4.ps
SBreak	Segment0	Comp0
 1993.6268	 -24.0998	Segment0	Comp0
 2000.5537	  44.8113	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001/-30/50 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"East (mm)":WeSn  \
-O >> nhd4.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nhd4.ps
1993.6268	-8.5	8.0
1993.6295	-25.0	7.2
1993.6487	-24.3	7.8
1993.6514	-28.9	7.3
1994.6289	-23.2	9.8
1994.6316	-26.9	8.7
2000.5510	41.4	8.0
2000.5537	42.8	7.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> nhd4.ps
1993.6268	-8.5	8.0
1993.6295	-25.0	7.2
1993.6487	-24.3	7.8
1993.6514	-28.9	7.3
1994.6289	-23.2	9.8
1994.6316	-26.9	8.7
2000.5510	41.4	8.0
2000.5537	42.8	7.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> nhd4.ps
SBreak	Segment0	Comp1
 1993.6268	 -27.2346	Segment0	Comp1
 2000.5537	  41.3328	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001/-30/40 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"Up (mm)":WeSn  \
-O >> nhd4.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> nhd4.ps
1993.6268	-8.8	22.9
1993.6295	-12.9	21.9
1993.6487	-13.1	23.5
1993.6514	-29.3	21.5
1994.6289	8.5	23.5
1994.6316	-17.8	26.4
2000.5510	31.7	23.6
2000.5537	28.5	22.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> nhd4.ps
1993.6268	-8.8	22.9
1993.6295	-12.9	21.9
1993.6487	-13.1	23.5
1993.6514	-29.3	21.5
1994.6289	8.5	23.5
1994.6316	-17.8	26.4
2000.5510	31.7	23.6
2000.5537	28.5	22.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> nhd4.ps
SBreak	Segment0	Comp2
 1993.6268	 -14.7724	Segment0	Comp2
 2000.5537	  30.3493	Segment0	Comp2
EOT