#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 3814.ps
1996.4852	-48.5	5.6
1997.9391	-33.3	5.8
1997.9610	-36.4	5.8
1999.0507	-5.8	5.9
1999.0534	-29.2	5.9
2001.0987	19.7	6.7
2001.1014	16.2	6.6
2004.6058	60.4	6.8
2005.1480	69.5	5.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 3814.ps
1996.4852	-48.5	5.6
1997.9391	-33.3	5.8
1997.9610	-36.4	5.8
1999.0507	-5.8	5.9
1999.0534	-29.2	5.9
2001.0987	19.7	6.7
2001.1014	16.2	6.6
2004.6058	60.4	6.8
2005.1480	69.5	5.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 3814.ps
SBreak	Segment0	Comp0
 1996.4852	 -52.2295	Segment0	Comp0
 2005.1480	  70.0924	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1996/2005.5/-80/60 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"3814 East (mm)":WeSn  \
-O >> 3814.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 3814.ps
1996.4852	53.4	7.3
1997.9391	37.3	6.9
1997.9610	20.5	7.1
1999.0507	18.2	7.6
1999.0534	19.7	7.6
2001.0987	-12.1	9.0
2001.1014	-14.1	9.9
2004.6058	-74.1	8.1
2005.1480	-79.5	7.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 3814.ps
1996.4852	53.4	7.3
1997.9391	37.3	6.9
1997.9610	20.5	7.1
1999.0507	18.2	7.6
1999.0534	19.7	7.6
2001.0987	-12.1	9.0
2001.1014	-14.1	9.9
2004.6058	-74.1	8.1
2005.1480	-79.5	7.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> 3814.ps
SBreak	Segment0	Comp1
 1996.4852	  54.7336	Segment0	Comp1
 2005.1480	 -78.9536	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1996/2005.5/-40/140 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"3814 Up (mm)":WeSn  \
-O >> 3814.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> 3814.ps
1996.4852	130.3	22.8
1997.9391	-11.6	20.4
1997.9610	-9.5	21.0
1999.0507	-38.8	20.3
1999.0534	22.8	20.5
2001.0987	-17.1	24.0
2001.1014	-17.1	23.8
2004.6058	-30.9	22.5
2005.1480	-15.5	19.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> 3814.ps
1996.4852	130.3	22.8
1997.9391	-11.6	20.4
1997.9610	-9.5	21.0
1999.0507	-38.8	20.3
1999.0534	22.8	20.5
2001.0987	-17.1	24.0
2001.1014	-17.1	23.8
2004.6058	-30.9	22.5
2005.1480	-15.5	19.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> 3814.ps
SBreak	Segment0	Comp2
 1996.4852	  30.4809	Segment0	Comp2
 2005.1480	 -39.1989	Segment0	Comp2
EOT