#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-20/20 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"AKUV North (mm)":WeSn  \
> akuv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akuv.ps
2001.6270	11.5	8.4
2002.4812	7.8	14.4
2002.4839	7.2	20.6
2002.4867	5.3	15.4
2002.4894	8.1	21.8
2004.4552	-10.9	17.9
2004.4580	-13.4	9.2
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akuv.ps
2001.6270	11.5	8.4
2002.4812	7.8	14.4
2002.4839	7.2	20.6
2002.4867	5.3	15.4
2002.4894	8.1	21.8
2004.4552	-10.9	17.9
2004.4580	-13.4	9.2
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akuv.ps
SBreak	Segment0	Comp0
 2001.6270	  13.8873	Segment0	Comp0
 2004.4580	 -11.8828	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/10 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"AKUV East (mm)":WeSn  \
-O >> akuv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akuv.ps
2001.6270	2.2	7.9
2002.4812	3.2	17.0
2002.4839	0.6	23.2
2002.4867	-3.8	17.6
2002.4894	-4.9	21.0
2004.4552	4.8	16.5
2004.4580	1.2	10.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akuv.ps
2001.6270	2.2	7.9
2002.4812	3.2	17.0
2002.4839	0.6	23.2
2002.4867	-3.8	17.6
2002.4894	-4.9	21.0
2004.4552	4.8	16.5
2004.4580	1.2	10.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akuv.ps
SBreak	Segment0	Comp1
 2001.6270	  -2.8017	Segment0	Comp1
 2004.4580	   2.4100	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-10/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"AKUV Up (mm)":WeSn  \
-O >> akuv.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akuv.ps
2001.6270	1.4	23.6
2002.4812	38.3	40.2
2002.4839	-1.9	59.9
2002.4867	5.0	42.0
2002.4894	-1.8	65.8
2004.4552	1.3	44.0
2004.4580	-10.0	29.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akuv.ps
2001.6270	1.4	23.6
2002.4812	38.3	40.2
2002.4839	-1.9	59.9
2002.4867	5.0	42.0
2002.4894	-1.8	65.8
2004.4552	1.3	44.0
2004.4580	-10.0	29.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> akuv.ps
SBreak	Segment0	Comp2
 2001.6270	   4.8857	Segment0	Comp2
 2004.4580	  -4.2567	Segment0	Comp2
EOT