#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1992.5/2004/-50/40 \
-G230/225/200 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"North (mm)":WeSn  \
> rose.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rose.ps
1992.6603	-49.3	7.1
1993.6432	-34.3	6.6
1995.5871	-19.8	7.2
1995.6036	-20.5	6.8
1996.6111	-9.0	7.3
1996.6138	-10.1	7.0
1998.6289	-0.2	6.8
1998.6317	-8.7	6.5
1998.6344	-1.9	6.3
1998.6399	-11.3	6.4
1998.6426	-11.2	6.6
2003.5599	31.1	6.7
2003.5627	30.2	6.1
2003.5654	33.2	6.7
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rose.ps
1992.6603	-49.3	7.1
1993.6432	-34.3	6.6
1995.5871	-19.8	7.2
1995.6036	-20.5	6.8
1996.6111	-9.0	7.3
1996.6138	-10.1	7.0
1998.6289	-0.2	6.8
1998.6317	-8.7	6.5
1998.6344	-1.9	6.3
1998.6399	-11.3	6.4
1998.6426	-11.2	6.6
2003.5599	31.1	6.7
2003.5627	30.2	6.1
2003.5654	33.2	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> rose.ps
SBreak	Segment0	Comp0
 1992.6603	 -40.4648	Segment0	Comp0
 2003.5654	  30.5770	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 >> rose.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rose.ps
1992.6603	-5.8	8.0
1993.6432	-4.8	7.2
1995.5871	6.3	7.8
1995.6036	4.0	7.3
1996.6111	9.3	9.8
1996.6138	2.9	8.7
1998.6289	-4.2	8.0
1998.6317	1.8	7.5
1998.6344	-2.8	7.4
1998.6399	-0.3	7.4
1998.6426	7.5	7.9
2003.5599	-8.4	7.9
2003.5627	-3.4	7.1
2003.5654	0.2	8.0
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rose.ps
1992.6603	-5.8	8.0
1993.6432	-4.8	7.2
1995.5871	6.3	7.8
1995.6036	4.0	7.3
1996.6111	9.3	9.8
1996.6138	2.9	8.7
1998.6289	-4.2	8.0
1998.6317	1.8	7.5
1998.6344	-2.8	7.4
1998.6399	-0.3	7.4
1998.6426	7.5	7.9
2003.5599	-8.4	7.9
2003.5627	-3.4	7.1
2003.5654	0.2	8.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> rose.ps
SBreak	Segment0	Comp1
 1992.6603	   3.0013	Segment0	Comp1
 2003.5654	  -1.8942	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> rose.ps
1992.6603	15.8	22.9
1993.6432	-12.5	21.9
1995.5871	12.2	23.5
1995.6036	8.8	21.5
1996.6111	-2.9	23.5
1996.6138	-12.7	26.4
1998.6289	-22.0	23.6
1998.6317	0.8	22.3
1998.6344	-27.2	22.8
1998.6399	1.3	22.6
1998.6426	26.4	24.1
2003.5599	6.3	24.1
2003.5627	7.3	21.3
2003.5654	4.5	23.1
EOT


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.06 -G255/0/0 -O \
-K >> rose.ps
1992.6603	15.8	22.9
1993.6432	-12.5	21.9
1995.5871	12.2	23.5
1995.6036	8.8	21.5
1996.6111	-2.9	23.5
1996.6138	-12.7	26.4
1998.6289	-22.0	23.6
1998.6317	0.8	22.3
1998.6344	-27.2	22.8
1998.6399	1.3	22.6
1998.6426	26.4	24.1
2003.5599	6.3	24.1
2003.5627	7.3	21.3
2003.5654	4.5	23.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> rose.ps
SBreak	Segment0	Comp2
 1992.6603	  -6.1930	Segment0	Comp2
 2003.5654	   4.7601	Segment0	Comp2
EOT