#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-320/230 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f80a160:"PRH2 North (mm)":WeSn  \
> prh2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> prh2.ps
1993.0328	-311.7	3.6
1994.3003	-275.2	3.1
1995.7131	-237.9	6.2
1996.8275	-201.9	5.6
1997.8789	-164.2	6.1
1997.9172	-163.6	5.9
1998.9165	-119.1	6.6
1998.9357	-123.5	2.9
2001.0247	-58.1	2.9
2001.0275	-61.6	3.1
2009.1891	219.5	2.8
2009.1918	219.7	2.8
2009.1945	220.6	2.8
2009.1973	222.2	3.3
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> prh2.ps
1993.0328	-311.7	3.6
1994.3003	-275.2	3.1
1995.7131	-237.9	6.2
1996.8275	-201.9	5.6
1997.8789	-164.2	6.1
1997.9172	-163.6	5.9
1998.9165	-119.1	6.6
1998.9357	-123.5	2.9
2001.0247	-58.1	2.9
2001.0275	-61.6	3.1
2009.1891	219.5	2.8
2009.1918	219.7	2.8
2009.1945	220.6	2.8
2009.1973	222.2	3.3
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> prh2.ps
SBreak	Segment0	Comp0
 1993.0328	-321.1193	Segment0	Comp0
 2009.1973	 219.0267	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-130/260 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"PRH2 East (mm)":WeSn  \
-O >> prh2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> prh2.ps
1993.0328	254.2	7.3
1994.3003	220.4	5.0
1995.7131	190.6	8.3
1996.8275	162.0	7.1
1997.8789	145.0	7.5
1997.9172	145.3	7.4
1998.9165	124.3	11.2
1998.9357	121.7	4.4
2001.0247	82.5	4.1
2001.0275	76.1	4.1
2009.1891	-124.5	3.3
2009.1918	-124.1	3.3
2009.1945	-123.8	3.3
2009.1973	-125.1	4.0
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> prh2.ps
1993.0328	254.2	7.3
1994.3003	220.4	5.0
1995.7131	190.6	8.3
1996.8275	162.0	7.1
1997.8789	145.0	7.5
1997.9172	145.3	7.4
1998.9165	124.3	11.2
1998.9357	121.7	4.4
2001.0247	82.5	4.1
2001.0275	76.1	4.1
2009.1891	-124.5	3.3
2009.1918	-124.1	3.3
2009.1945	-123.8	3.3
2009.1973	-125.1	4.0
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> prh2.ps
SBreak	Segment0	Comp1
 1993.0328	 259.6591	Segment0	Comp1
 2009.1973	-122.9463	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-30/60 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"PRH2 Up (mm)":WeSn  \
-O >> prh2.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> prh2.ps
1993.0328	44.6	20.0
1994.3003	54.7	14.6
1995.7131	32.2	22.3
1996.8275	-12.3	20.6
1997.8789	-13.4	20.8
1997.9172	10.6	21.0
1998.9165	16.7	23.5
1998.9357	-1.4	11.0
2001.0247	3.2	10.9
2001.0275	14.6	11.3
2009.1891	-13.0	10.2
2009.1918	-14.4	10.2
2009.1945	-14.5	10.2
2009.1973	-21.0	11.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> prh2.ps
1993.0328	44.6	20.0
1994.3003	54.7	14.6
1995.7131	32.2	22.3
1996.8275	-12.3	20.6
1997.8789	-13.4	20.8
1997.9172	10.6	21.0
1998.9165	16.7	23.5
1998.9357	-1.4	11.0
2001.0247	3.2	10.9
2001.0275	14.6	11.3
2009.1891	-13.0	10.2
2009.1918	-14.4	10.2
2009.1945	-14.5	10.2
2009.1973	-21.0	11.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> prh2.ps
SBreak	Segment0	Comp2
 1993.0328	  30.4858	Segment0	Comp2
 2009.1973	 -16.9464	Segment0	Comp2
EOT