#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hdqt.ps
2001.6709	0.4	7.1
2005.6655	1.1	7.1
2006.6840	-1.3	4.9
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hdqt.ps
2001.6709	0.4	7.1
2005.6655	1.1	7.1
2006.6840	-1.3	4.9
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hdqt.ps
SBreak	Segment0	Comp0
 2001.6709	   0.6003	Segment0	Comp0
 2006.6840	  -0.4407	Segment0	Comp0
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hdqt.ps
2001.6709	2.5	8.1
2005.6655	-0.7	8.4
2006.6840	-1.1	6.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hdqt.ps
2001.6709	2.5	8.1
2005.6655	-0.7	8.4
2006.6840	-1.1	6.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> hdqt.ps
SBreak	Segment0	Comp1
 2001.6709	   2.3985	Segment0	Comp1
 2006.6840	  -1.2624	Segment0	Comp1
EOT


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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> hdqt.ps
2001.6709	5.5	21.9
2005.6655	-11.3	20.8
2006.6840	5.7	18.5
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> hdqt.ps
2001.6709	5.5	21.9
2005.6655	-11.3	20.8
2006.6840	5.7	18.5
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> hdqt.ps
SBreak	Segment0	Comp2
 2001.6709	   3.4068	Segment0	Comp2
 2006.6840	  -2.2505	Segment0	Comp2
EOT