#!/bin/csh

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

psbasemap -JX6/1.5 \
-R2001.5/2004.5/-30/80 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f20a40:"AKPS North (mm)":WeSn  \
> akps.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akps.ps
2001.6160	72.9	8.4
2001.6187	47.0	14.4
2002.4784	31.1	20.6
2002.4812	31.7	15.4
2004.4580	-28.5	21.8
2004.4607	-26.2	17.9
2004.4634	-27.2	9.2
2004.4662	-23.9	6.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akps.ps
2001.6160	72.9	8.4
2001.6187	47.0	14.4
2002.4784	31.1	20.6
2002.4812	31.7	15.4
2004.4580	-28.5	21.8
2004.4607	-26.2	17.9
2004.4634	-27.2	9.2
2004.4662	-23.9	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akps.ps
SBreak	Segment0	Comp0
 2001.6160	  51.5392	Segment0	Comp0
 2004.4662	 -26.0123	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akps.ps
2001.6160	-28.2	7.9
2001.6187	-1.0	17.0
2002.4784	-3.1	23.2
2002.4812	-3.7	17.6
2004.4580	1.5	21.0
2004.4607	-0.3	16.5
2004.4634	1.2	10.9
2004.4662	3.8	7.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akps.ps
2001.6160	-28.2	7.9
2001.6187	-1.0	17.0
2002.4784	-3.1	23.2
2002.4812	-3.7	17.6
2004.4580	1.5	21.0
2004.4607	-0.3	16.5
2004.4634	1.2	10.9
2004.4662	3.8	7.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> akps.ps
SBreak	Segment0	Comp1
 2001.6160	  -3.5374	Segment0	Comp1
 2004.4662	   1.4195	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> akps.ps
2001.6160	19.8	23.6
2001.6187	17.1	40.2
2002.4784	-5.1	59.9
2002.4812	-11.0	42.0
2004.4580	7.3	65.8
2004.4607	7.0	44.0
2004.4634	-1.8	29.1
2004.4662	-13.9	22.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> akps.ps
2001.6160	19.8	23.6
2001.6187	17.1	40.2
2002.4784	-5.1	59.9
2002.4812	-11.0	42.0
2004.4580	7.3	65.8
2004.4607	7.0	44.0
2004.4634	-1.8	29.1
2004.4662	-13.9	22.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> akps.ps
SBreak	Segment0	Comp2
 2001.6160	   4.1582	Segment0	Comp2
 2004.4662	  -2.1066	Segment0	Comp2
EOT