#!/bin/csh

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

psbasemap -JX6/1.5 \
-R1994/2005.5/-100/130 \
-G236/232/211 -Y5.25 -P -K \
-Bf0.25a1/f40a80:"WHAL North (mm)":WeSn  \
> whal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whal.ps
1994.3140	-95.5	5.6
1995.4831	-83.3	5.8
1996.4277	-69.4	5.8
1997.9418	-31.2	5.9
1997.9446	-31.8	5.9
1999.0288	-8.5	6.7
1999.0315	-7.9	6.6
2001.0795	30.1	6.8
2001.0822	24.2	5.2
2004.3402	102.8	3.3
2005.1891	123.2	6.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whal.ps
1994.3140	-95.5	5.6
1995.4831	-83.3	5.8
1996.4277	-69.4	5.8
1997.9418	-31.2	5.9
1997.9446	-31.8	5.9
1999.0288	-8.5	6.7
1999.0315	-7.9	6.6
2001.0795	30.1	6.8
2001.0822	24.2	5.2
2004.3402	102.8	3.3
2005.1891	123.2	6.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whal.ps
SBreak	Segment0	Comp0
 1994.3140	-107.2879	Segment0	Comp0
 2005.1891	 118.6827	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R1994/2005.5/-220/210 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f80a160:"WHAL East (mm)":WeSn  \
-O >> whal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whal.ps
1994.3140	203.9	7.3
1995.4831	133.9	6.9
1996.4277	137.2	7.1
1997.9418	61.3	7.6
1997.9446	61.3	7.6
1999.0288	24.2	9.0
1999.0315	20.6	9.9
2001.0795	-62.8	8.1
2001.0822	-62.1	7.5
2004.3402	-188.4	9.0
2005.1891	-218.5	7.7
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whal.ps
1994.3140	203.9	7.3
1995.4831	133.9	6.9
1996.4277	137.2	7.1
1997.9418	61.3	7.6
1997.9446	61.3	7.6
1999.0288	24.2	9.0
1999.0315	20.6	9.9
2001.0795	-62.8	8.1
2001.0822	-62.1	7.5
2004.3402	-188.4	9.0
2005.1891	-218.5	7.7
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> whal.ps
SBreak	Segment0	Comp1
 1994.3140	 200.2455	Segment0	Comp1
 2005.1891	-218.7611	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R1994/2005.5/-50/160 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f40a80:"WHAL Up (mm)":WeSn  \
-O >> whal.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> whal.ps
1994.3140	-13.6	22.8
1995.4831	8.9	20.4
1996.4277	154.5	21.0
1997.9418	-7.9	20.3
1997.9446	-8.0	20.5
1999.0288	-3.5	24.0
1999.0315	-27.3	23.8
2001.0795	-49.3	22.5
2001.0822	-18.2	19.7
2004.3402	-9.5	16.0
2005.1891	-8.0	21.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> whal.ps
1994.3140	-13.6	22.8
1995.4831	8.9	20.4
1996.4277	154.5	21.0
1997.9418	-7.9	20.3
1997.9446	-8.0	20.5
1999.0288	-3.5	24.0
1999.0315	-27.3	23.8
2001.0795	-49.3	22.5
2001.0822	-18.2	19.7
2004.3402	-9.5	16.0
2005.1891	-8.0	21.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> whal.ps
SBreak	Segment0	Comp2
 1994.3140	  24.3786	Segment0	Comp2
 2005.1891	 -25.5326	Segment0	Comp2
EOT