#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whlc.ps
2001.6270	15.2	8.4
2002.4702	3.6	14.4
2002.4730	6.4	20.6
2002.4757	11.3	15.4
2004.4580	-5.8	21.8
2004.4607	-6.3	17.9
2004.4634	-9.2	9.2
2004.4662	-6.1	6.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whlc.ps
2001.6270	15.2	8.4
2002.4702	3.6	14.4
2002.4730	6.4	20.6
2002.4757	11.3	15.4
2004.4580	-5.8	21.8
2004.4607	-6.3	17.9
2004.4634	-9.2	9.2
2004.4662	-6.1	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whlc.ps
SBreak	Segment0	Comp0
 2001.6270	  13.9431	Segment0	Comp0
 2004.4662	  -6.9731	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:"WHLC East (mm)":WeSn  \
-O >> whlc.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whlc.ps
2001.6270	-1.2	7.9
2002.4702	-0.9	17.0
2002.4730	1.0	23.2
2002.4757	-9.1	17.6
2004.4580	-0.9	21.0
2004.4607	1.0	16.5
2004.4634	1.8	10.9
2004.4662	4.6	7.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whlc.ps
2001.6270	-1.2	7.9
2002.4702	-0.9	17.0
2002.4730	1.0	23.2
2002.4757	-9.1	17.6
2004.4580	-0.9	21.0
2004.4607	1.0	16.5
2004.4634	1.8	10.9
2004.4662	4.6	7.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whlc.ps
SBreak	Segment0	Comp1
 2001.6270	  -3.9052	Segment0	Comp1
 2004.4662	   1.6217	Segment0	Comp1
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:"WHLC Up (mm)":WeSn  \
-O >> whlc.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whlc.ps
2001.6270	8.8	23.6
2002.4702	3.7	40.2
2002.4730	9.4	59.9
2002.4757	-1.3	42.0
2004.4580	-6.5	65.8
2004.4607	-1.3	44.0
2004.4634	2.7	29.1
2004.4662	-9.8	22.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whlc.ps
2001.6270	8.8	23.6
2002.4702	3.7	40.2
2002.4730	9.4	59.9
2002.4757	-1.3	42.0
2004.4580	-6.5	65.8
2004.4607	-1.3	44.0
2004.4634	2.7	29.1
2004.4662	-9.8	22.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> whlc.ps
SBreak	Segment0	Comp2
 2001.6270	   7.7229	Segment0	Comp2
 2004.4662	  -3.9114	Segment0	Comp2
EOT