#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1993/2010/-180/140 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"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	-179.6	3.6
1994.3003	-160.6	3.1
1995.7131	-142.8	6.2
1996.8275	-122.0	5.6
1997.8789	-98.8	6.1
1997.9172	-98.7	5.9
1998.9165	-67.9	6.6
1998.9357	-72.6	2.9
2001.0247	-35.8	2.9
2001.0275	-39.5	3.1
2009.1891	129.6	2.8
2009.1918	129.8	2.8
2009.1945	130.6	2.8
2009.1973	132.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	-179.6	3.6
1994.3003	-160.6	3.1
1995.7131	-142.8	6.2
1996.8275	-122.0	5.6
1997.8789	-98.8	6.1
1997.9172	-98.7	5.9
1998.9165	-67.9	6.6
1998.9357	-72.6	2.9
2001.0247	-35.8	2.9
2001.0275	-39.5	3.1
2009.1891	129.6	2.8
2009.1918	129.8	2.8
2009.1945	130.6	2.8
2009.1973	132.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	-189.1453	Segment0	Comp0
 2009.1973	 129.0108	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1993/2010/-190/390 \
-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	381.9	7.3
1994.3003	333.4	5.0
1995.7131	287.1	8.3
1996.8275	245.6	7.1
1997.8789	216.3	7.5
1997.9172	216.2	7.4
1998.9165	183.5	11.2
1998.9357	180.7	4.4
2001.0247	117.1	4.1
2001.0275	110.7	4.1
2009.1891	-184.9	3.3
2009.1918	-184.5	3.3
2009.1945	-184.3	3.3
2009.1973	-185.6	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	381.9	7.3
1994.3003	333.4	5.0
1995.7131	287.1	8.3
1996.8275	245.6	7.1
1997.8789	216.3	7.5
1997.9172	216.2	7.4
1998.9165	183.5	11.2
1998.9357	180.7	4.4
2001.0247	117.1	4.1
2001.0275	110.7	4.1
2009.1891	-184.9	3.3
2009.1918	-184.5	3.3
2009.1945	-184.3	3.3
2009.1973	-185.6	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	 387.3830	Segment0	Comp1
 2009.1973	-183.4225	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	45.0	20.0
1994.3003	55.2	14.6
1995.7131	32.5	22.3
1996.8275	-12.0	20.6
1997.8789	-13.2	20.8
1997.9172	10.8	21.0
1998.9165	16.8	23.5
1998.9357	-1.1	11.0
2001.0247	3.3	10.9
2001.0275	14.6	11.3
2009.1891	-13.2	10.2
2009.1918	-14.7	10.2
2009.1945	-14.8	10.2
2009.1973	-21.3	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	45.0	20.0
1994.3003	55.2	14.6
1995.7131	32.5	22.3
1996.8275	-12.0	20.6
1997.8789	-13.2	20.8
1997.9172	10.8	21.0
1998.9165	16.8	23.5
1998.9357	-1.1	11.0
2001.0247	3.3	10.9
2001.0275	14.6	11.3
2009.1891	-13.2	10.2
2009.1918	-14.7	10.2
2009.1945	-14.8	10.2
2009.1973	-21.3	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.9549	Segment0	Comp2
 2009.1973	 -17.2072	Segment0	Comp2
EOT