#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994.5/2001/-10/30 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"HAMI North (mm)":WeSn  \
> hami.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hami.ps
1994.6371	20.5	6.9
1998.6372	-7.9	6.4
1998.6454	6.9	7.0
1998.6481	-4.3	6.5
1998.9603	-1.3	7.1
2000.5756	-4.8	6.8
2000.5783	-7.6	6.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hami.ps
1994.6371	20.5	6.9
1998.6372	-7.9	6.4
1998.6454	6.9	7.0
1998.6481	-4.3	6.5
1998.9603	-1.3	7.1
2000.5756	-4.8	6.8
2000.5783	-7.6	6.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hami.ps
SBreak	Segment0	Comp0
 1994.6371	  18.3125	Segment0	Comp0
 2000.5783	  -8.1254	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994.5/2001/-20/30 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"HAMI East (mm)":WeSn  \
-O >> hami.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hami.ps
1994.6371	2.5	7.0
1998.6372	8.4	6.8
1998.6454	-8.4	7.1
1998.6481	22.0	7.0
1998.9603	-4.6	9.5
2000.5756	-9.7	8.5
2000.5783	-12.4	7.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hami.ps
1994.6371	2.5	7.0
1998.6372	8.4	6.8
1998.6454	-8.4	7.1
1998.6481	22.0	7.0
1998.9603	-4.6	9.5
2000.5756	-9.7	8.5
2000.5783	-12.4	7.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hami.ps
SBreak	Segment0	Comp1
 1994.6371	   9.6139	Segment0	Comp1
 2000.5783	  -4.4165	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1994.5/2001/-40/50 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"HAMI Up (mm)":WeSn  \
-O >> hami.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hami.ps
1994.6371	-32.8	22.2
1998.6372	38.2	21.1
1998.6454	-31.6	22.6
1998.6481	48.0	20.1
1998.9603	-29.2	22.4
2000.5756	11.8	24.9
2000.5783	2.9	22.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hami.ps
1994.6371	-32.8	22.2
1998.6372	38.2	21.1
1998.6454	-31.6	22.6
1998.6481	48.0	20.1
1998.9603	-29.2	22.4
2000.5756	11.8	24.9
2000.5783	2.9	22.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> hami.ps
SBreak	Segment0	Comp2
 1994.6371	 -26.1531	Segment0	Comp2
 2000.5783	  12.8975	Segment0	Comp2
EOT