#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cast.ps
1993.6076	-18.7	7.1
1993.6104	-17.0	6.6
1994.6261	-0.2	7.2
2000.5510	29.5	6.8
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cast.ps
1993.6076	-18.7	7.1
1993.6104	-17.0	6.6
1994.6261	-0.2	7.2
2000.5510	29.5	6.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cast.ps
SBreak	Segment0	Comp0
 1993.6076	 -14.2013	Segment0	Comp0
 2000.5510	  30.5196	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cast.ps
1993.6076	1.5	8.0
1993.6104	-14.3	7.2
1994.6261	-10.3	7.8
2000.5510	11.8	7.3
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cast.ps
1993.6076	1.5	8.0
1993.6104	-14.3	7.2
1994.6261	-10.3	7.8
2000.5510	11.8	7.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> cast.ps
SBreak	Segment0	Comp1
 1993.6076	  -8.1417	Segment0	Comp1
 2000.5510	  11.3949	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993.5/2001/-20/10 \
-G230/225/200 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"Up (mm)":WeSn  \
-O >> cast.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> cast.ps
1993.6076	9.3	22.9
1993.6104	-13.3	21.9
1994.6261	-3.8	23.5
2000.5510	4.7	21.5
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> cast.ps
1993.6076	9.3	22.9
1993.6104	-13.3	21.9
1994.6261	-3.8	23.5
2000.5510	4.7	21.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> cast.ps
SBreak	Segment0	Comp2
 1993.6076	  -2.1336	Segment0	Comp2
 2000.5510	   4.4028	Segment0	Comp2
EOT