#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-10/10 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f5a10:"AKSO North (mm)":WeSn  \
> akso.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akso.ps
2001.6160	6.0	8.4
2002.4867	-2.4	14.4
2002.4894	3.6	20.6
2002.4921	-8.1	15.4
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akso.ps
2001.6160	6.0	8.4
2002.4867	-2.4	14.4
2002.4894	3.6	20.6
2002.4921	-8.1	15.4
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akso.ps
SBreak	Segment0	Comp0
 2001.6160	   5.9945	Segment0	Comp0
 2002.4921	  -2.1081	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akso.ps
2001.6160	0.2	7.9
2002.4867	-0.2	17.0
2002.4894	0.1	23.2
2002.4921	-0.1	17.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akso.ps
2001.6160	0.2	7.9
2002.4867	-0.2	17.0
2002.4894	0.1	23.2
2002.4921	-0.1	17.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akso.ps
SBreak	Segment0	Comp1
 2001.6160	   0.2124	Segment0	Comp1
 2002.4921	  -0.0578	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2001.5/2002.5/-20/20 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f5a10:"AKSO Up (mm)":WeSn  \
-O >> akso.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akso.ps
2001.6160	8.5	23.6
2002.4867	11.7	40.2
2002.4894	-12.4	59.9
2002.4921	-5.5	42.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akso.ps
2001.6160	8.5	23.6
2002.4867	11.7	40.2
2002.4894	-12.4	59.9
2002.4921	-5.5	42.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> akso.ps
SBreak	Segment0	Comp2
 2001.6160	   8.5112	Segment0	Comp2
 2002.4921	  -2.5001	Segment0	Comp2
EOT