#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994/2001.5/-150/80 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"PALO North (mm)":WeSn  \
> palo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> palo.ps
1994.3140	-144.8	5.6
1995.4831	-113.2	5.8
1996.4277	-79.2	5.8
1997.9418	-29.0	5.9
1997.9446	-39.4	5.9
1999.0288	3.5	6.7
1999.0315	-1.6	6.6
2001.0822	71.0	6.8
2001.0850	76.1	5.2
2001.0877	73.9	3.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> palo.ps
1994.3140	-144.8	5.6
1995.4831	-113.2	5.8
1996.4277	-79.2	5.8
1997.9418	-29.0	5.9
1997.9446	-39.4	5.9
1999.0288	3.5	6.7
1999.0315	-1.6	6.6
2001.0822	71.0	6.8
2001.0850	76.1	5.2
2001.0877	73.9	3.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> palo.ps
SBreak	Segment0	Comp0
 1994.3140	-153.1813	Segment0	Comp0
 2001.0877	  71.8859	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994/2001.5/-60/130 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"PALO East (mm)":WeSn  \
-O >> palo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> palo.ps
1994.3140	123.3	7.3
1995.4831	71.9	6.9
1996.4277	78.0	7.1
1997.9418	27.7	7.6
1997.9446	25.9	7.6
1999.0288	4.1	9.0
1999.0315	0.1	9.9
2001.0822	-58.4	8.1
2001.0850	-59.4	7.5
2001.0877	-58.4	9.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> palo.ps
1994.3140	123.3	7.3
1995.4831	71.9	6.9
1996.4277	78.0	7.1
1997.9418	27.7	7.6
1997.9446	25.9	7.6
1999.0288	4.1	9.0
1999.0315	0.1	9.9
2001.0822	-58.4	8.1
2001.0850	-59.4	7.5
2001.0877	-58.4	9.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> palo.ps
SBreak	Segment0	Comp1
 1994.3140	 120.9220	Segment0	Comp1
 2001.0877	 -56.9569	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1994/2001.5/-50/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"PALO Up (mm)":WeSn  \
-O >> palo.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> palo.ps
1994.3140	36.1	22.8
1995.4831	-47.0	20.4
1996.4277	30.9	21.0
1997.9418	3.4	20.3
1997.9446	-2.4	20.5
1999.0288	21.2	24.0
1999.0315	10.3	23.8
2001.0822	-6.9	22.5
2001.0850	-16.9	19.7
2001.0877	-5.0	16.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> palo.ps
1994.3140	36.1	22.8
1995.4831	-47.0	20.4
1996.4277	30.9	21.0
1997.9418	3.4	20.3
1997.9446	-2.4	20.5
1999.0288	21.2	24.0
1999.0315	10.3	23.8
2001.0822	-6.9	22.5
2001.0850	-16.9	19.7
2001.0877	-5.0	16.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> palo.ps
SBreak	Segment0	Comp2
 1994.3140	  13.9455	Segment0	Comp2
 2001.0877	  -6.6032	Segment0	Comp2
EOT