#!/bin/csh

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

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

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dnl3.ps
2002.8619	-4.1	4.7
2002.8646	-2.4	4.1
2002.8673	-4.6	4.2
2002.8701	-6.1	4.1
2002.8728	-7.1	3.8
2002.8756	-3.3	3.8
2002.8783	-4.0	3.8
2002.8810	-4.4	3.7
2002.8838	-1.0	3.8
2002.8865	-1.2	3.6
2002.8892	-1.9	3.8
2002.8920	-0.2	4.4
2002.8947	-1.5	3.9
2002.8975	-2.4	3.9
2002.9002	0.3	3.7
2002.9029	-3.7	4.1
2002.9057	-4.0	3.8
2002.9084	-4.2	4.0
2002.9112	-7.0	3.5
2002.9139	-9.1	3.6
2002.9194	-7.9	3.8
2002.9221	-2.6	3.9
2002.9248	-2.3	4.4
2002.9276	-1.4	3.8
2002.9303	-1.7	3.8
2002.9331	-2.2	3.7
2003.3683	1.1	3.8
2003.3710	1.9	4.1
2003.3737	0.2	3.5
2003.3847	2.3	3.6
2003.3874	2.2	3.4
2003.7351	2.1	3.4
2003.7379	2.9	3.4
2003.7406	4.8	3.4
2003.7434	3.0	4.9
2003.7461	3.5	3.6
2003.7488	6.6	3.6
2003.7516	10.9	3.5
2004.4607	11.4	3.4
2004.4634	5.2	3.6
2004.4662	8.4	3.7
2004.4689	5.9	3.7
2004.4716	6.5	3.6
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dnl3.ps
2002.8619	-4.1	4.7
2002.8646	-2.4	4.1
2002.8673	-4.6	4.2
2002.8701	-6.1	4.1
2002.8728	-7.1	3.8
2002.8756	-3.3	3.8
2002.8783	-4.0	3.8
2002.8810	-4.4	3.7
2002.8838	-1.0	3.8
2002.8865	-1.2	3.6
2002.8892	-1.9	3.8
2002.8920	-0.2	4.4
2002.8947	-1.5	3.9
2002.8975	-2.4	3.9
2002.9002	0.3	3.7
2002.9029	-3.7	4.1
2002.9057	-4.0	3.8
2002.9084	-4.2	4.0
2002.9112	-7.0	3.5
2002.9139	-9.1	3.6
2002.9194	-7.9	3.8
2002.9221	-2.6	3.9
2002.9248	-2.3	4.4
2002.9276	-1.4	3.8
2002.9303	-1.7	3.8
2002.9331	-2.2	3.7
2003.3683	1.1	3.8
2003.3710	1.9	4.1
2003.3737	0.2	3.5
2003.3847	2.3	3.6
2003.3874	2.2	3.4
2003.7351	2.1	3.4
2003.7379	2.9	3.4
2003.7406	4.8	3.4
2003.7434	3.0	4.9
2003.7461	3.5	3.6
2003.7488	6.6	3.6
2003.7516	10.9	3.5
2004.4607	11.4	3.4
2004.4634	5.2	3.6
2004.4662	8.4	3.7
2004.4689	5.9	3.7
2004.4716	6.5	3.6
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> dnl3.ps
SBreak	Segment0	Comp0
 2002.8619	  -3.4072	Segment0	Comp0
 2004.4716	   8.7654	Segment0	Comp0
EOT


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

psbasemap -JX6/1.5 \
-R2002.5/2004.5/-40/50 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f20a40:"DNL3 East (mm)":WeSn  \
-O >> dnl3.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dnl3.ps
2002.8619	-32.1	6.5
2002.8646	-29.0	4.5
2002.8673	-31.1	4.5
2002.8701	-31.5	4.0
2002.8728	-27.2	3.8
2002.8756	-26.2	3.8
2002.8783	-23.3	3.8
2002.8810	-29.0	3.7
2002.8838	-23.9	3.7
2002.8865	-22.3	3.4
2002.8892	-22.1	3.9
2002.8920	-24.6	4.7
2002.8947	-19.5	3.7
2002.8975	-14.8	4.1
2002.9002	-18.5	3.5
2002.9029	-15.9	3.8
2002.9057	-21.2	3.8
2002.9084	-14.0	4.3
2002.9112	-13.7	3.6
2002.9139	-15.9	3.6
2002.9194	-12.4	3.7
2002.9221	-9.1	3.8
2002.9248	-11.2	4.3
2002.9276	-11.2	3.7
2002.9303	-14.0	4.0
2002.9331	-13.2	3.8
2003.3683	12.4	3.8
2003.3710	15.8	4.3
2003.3737	13.4	3.7
2003.3847	16.2	3.8
2003.3874	12.1	3.5
2003.7351	27.3	3.5
2003.7379	29.1	3.3
2003.7406	30.4	3.5
2003.7434	34.0	5.1
2003.7461	32.2	4.0
2003.7488	31.7	3.6
2003.7516	31.0	3.5
2004.4607	35.8	3.5
2004.4634	34.1	4.3
2004.4662	40.7	3.8
2004.4689	38.9	3.7
2004.4716	40.1	3.8
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dnl3.ps
2002.8619	-32.1	6.5
2002.8646	-29.0	4.5
2002.8673	-31.1	4.5
2002.8701	-31.5	4.0
2002.8728	-27.2	3.8
2002.8756	-26.2	3.8
2002.8783	-23.3	3.8
2002.8810	-29.0	3.7
2002.8838	-23.9	3.7
2002.8865	-22.3	3.4
2002.8892	-22.1	3.9
2002.8920	-24.6	4.7
2002.8947	-19.5	3.7
2002.8975	-14.8	4.1
2002.9002	-18.5	3.5
2002.9029	-15.9	3.8
2002.9057	-21.2	3.8
2002.9084	-14.0	4.3
2002.9112	-13.7	3.6
2002.9139	-15.9	3.6
2002.9194	-12.4	3.7
2002.9221	-9.1	3.8
2002.9248	-11.2	4.3
2002.9276	-11.2	3.7
2002.9303	-14.0	4.0
2002.9331	-13.2	3.8
2003.3683	12.4	3.8
2003.3710	15.8	4.3
2003.3737	13.4	3.7
2003.3847	16.2	3.8
2003.3874	12.1	3.5
2003.7351	27.3	3.5
2003.7379	29.1	3.3
2003.7406	30.4	3.5
2003.7434	34.0	5.1
2003.7461	32.2	4.0
2003.7488	31.7	3.6
2003.7516	31.0	3.5
2004.4607	35.8	3.5
2004.4634	34.1	4.3
2004.4662	40.7	3.8
2004.4689	38.9	3.7
2004.4716	40.1	3.8
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
-K >> dnl3.ps
SBreak	Segment0	Comp1
 2002.8619	 -18.4966	Segment0	Comp1
 2004.4716	  49.9937	Segment0	Comp1
EOT


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

psbasemap -JX6/1.5 \
-R2002.5/2004.5/-30/40 \
-G236/232/211 -Y-1.75 -P -K \
-Bf0.25a1/f10a20:"DNL3 Up (mm)":WeSn  \
-O >> dnl3.ps

#  Draw error bars
#  ---------------
psxy << EOT \
-JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \
-K >> dnl3.ps
2002.8619	-8.1	16.4
2002.8646	-13.7	14.8
2002.8673	-2.8	16.2
2002.8701	-13.2	13.9
2002.8728	-5.8	13.0
2002.8756	-1.7	13.3
2002.8783	1.7	13.3
2002.8810	0.2	12.9
2002.8838	10.6	13.6
2002.8865	2.7	13.2
2002.8892	-9.0	13.4
2002.8920	-6.3	15.7
2002.8947	-1.4	14.1
2002.8975	13.0	14.3
2002.9002	-1.5	13.5
2002.9029	6.7	14.5
2002.9057	9.8	13.6
2002.9084	3.5	14.7
2002.9112	-1.4	12.9
2002.9139	-3.8	13.8
2002.9194	-11.2	14.0
2002.9221	4.1	14.3
2002.9248	4.6	16.6
2002.9276	7.1	13.6
2002.9303	16.8	14.0
2002.9331	-7.2	13.6
2003.3683	-12.2	13.9
2003.3710	-12.6	14.0
2003.3737	-6.2	13.1
2003.3847	-4.8	13.5
2003.3874	-10.3	12.7
2003.7351	4.0	13.0
2003.7379	-9.0	12.6
2003.7406	-1.5	12.7
2003.7434	11.2	16.8
2003.7461	9.8	13.7
2003.7488	25.3	13.0
2003.7516	30.9	12.4
2004.4607	14.7	12.7
2004.4634	-21.1	13.6
2004.4662	-7.1	13.4
2004.4689	-8.4	14.0
2004.4716	6.3	13.1
EOT


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


#  Draw points
#  -----------
psxy << EOT \
-JX -R -Sc0.03 -G255/0/0 -O \
-K >> dnl3.ps
2002.8619	-8.1	16.4
2002.8646	-13.7	14.8
2002.8673	-2.8	16.2
2002.8701	-13.2	13.9
2002.8728	-5.8	13.0
2002.8756	-1.7	13.3
2002.8783	1.7	13.3
2002.8810	0.2	12.9
2002.8838	10.6	13.6
2002.8865	2.7	13.2
2002.8892	-9.0	13.4
2002.8920	-6.3	15.7
2002.8947	-1.4	14.1
2002.8975	13.0	14.3
2002.9002	-1.5	13.5
2002.9029	6.7	14.5
2002.9057	9.8	13.6
2002.9084	3.5	14.7
2002.9112	-1.4	12.9
2002.9139	-3.8	13.8
2002.9194	-11.2	14.0
2002.9221	4.1	14.3
2002.9248	4.6	16.6
2002.9276	7.1	13.6
2002.9303	16.8	14.0
2002.9331	-7.2	13.6
2003.3683	-12.2	13.9
2003.3710	-12.6	14.0
2003.3737	-6.2	13.1
2003.3847	-4.8	13.5
2003.3874	-10.3	12.7
2003.7351	4.0	13.0
2003.7379	-9.0	12.6
2003.7406	-1.5	12.7
2003.7434	11.2	16.8
2003.7461	9.8	13.7
2003.7488	25.3	13.0
2003.7516	30.9	12.4
2004.4607	14.7	12.7
2004.4634	-21.1	13.6
2004.4662	-7.1	13.4
2004.4689	-8.4	14.0
2004.4716	6.3	13.1
EOT


#  Draw best fit line
#  ------------------
psxy << EOT \
-JX -R -MS -W3/102/102/255 -O \
>> dnl3.ps
SBreak	Segment0	Comp2
 2002.8619	  -0.1153	Segment0	Comp2
 2004.4716	   0.3084	Segment0	Comp2
EOT