#!/bin/csh # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R2001.5/2003.5/-20/30 \ -G230/225/200 -Y5.25 -P -K \ -Bf0.25a1/f10a20:"North (mm)":WeSn \ > dch1.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> dch1.ps 2001.6133 0.3 8.1 2001.6160 4.4 4.9 2001.6187 0.7 3.5 2001.6215 1.6 3.7 2001.6242 3.6 3.7 2001.6270 5.4 3.6 2002.4565 24.4 3.7 2002.4593 1.6 3.8 2002.4675 4.3 3.8 2002.4702 0.4 3.8 2002.4730 3.2 3.8 2002.4757 5.1 3.9 2002.4784 3.3 3.8 2002.4812 5.4 3.6 2002.4839 5.0 3.7 2002.4867 1.3 3.7 2002.4894 6.3 3.8 2002.4921 1.1 4.0 2002.4949 0.3 4.1 2002.4976 -2.1 4.0 2002.5003 2.4 3.9 2002.5031 -0.2 4.0 2002.5058 -0.9 4.0 2002.5086 0.9 3.7 2002.5113 0.2 4.0 2003.4805 -7.4 4.1 2003.4833 -7.5 3.8 2003.4860 -7.9 3.8 2003.4887 -4.1 3.8 2003.4915 -6.1 4.0 2003.4942 -11.0 4.0 2003.4969 -11.4 3.9 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> dch1.ps 2001.6133 0.3 8.1 2001.6160 4.4 4.9 2001.6187 0.7 3.5 2001.6215 1.6 3.7 2001.6242 3.6 3.7 2001.6270 5.4 3.6 2002.4565 24.4 3.7 2002.4593 1.6 3.8 2002.4675 4.3 3.8 2002.4702 0.4 3.8 2002.4730 3.2 3.8 2002.4757 5.1 3.9 2002.4784 3.3 3.8 2002.4812 5.4 3.6 2002.4839 5.0 3.7 2002.4867 1.3 3.7 2002.4894 6.3 3.8 2002.4921 1.1 4.0 2002.4949 0.3 4.1 2002.4976 -2.1 4.0 2002.5003 2.4 3.9 2002.5031 -0.2 4.0 2002.5058 -0.9 4.0 2002.5086 0.9 3.7 2002.5113 0.2 4.0 2003.4805 -7.4 4.1 2003.4833 -7.5 3.8 2003.4860 -7.9 3.8 2003.4887 -4.1 3.8 2003.4915 -6.1 4.0 2003.4942 -11.0 4.0 2003.4969 -11.4 3.9 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ -K >> dch1.ps SBreak Segment0 Comp0 2001.6133 5.4731 Segment0 Comp0 2003.4969 -5.7612 Segment0 Comp0 EOT # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R2001.5/2003.5/-10/60 \ -G230/225/200 -Y-1.75 -P -K \ -Bf0.25a1/f10a20:"East (mm)":WeSn \ -O >> dch1.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> dch1.ps 2001.6133 1.9 9.0 2001.6160 -0.8 5.1 2001.6187 -0.1 3.9 2001.6215 0.8 4.2 2001.6242 2.7 4.3 2001.6270 2.4 4.3 2002.4565 55.2 4.5 2002.4593 6.8 4.6 2002.4675 6.1 4.1 2002.4702 3.5 4.2 2002.4730 6.1 4.2 2002.4757 -3.4 4.5 2002.4784 -3.1 4.5 2002.4812 -3.7 4.2 2002.4839 2.4 4.3 2002.4867 -1.4 4.6 2002.4894 -1.7 5.1 2002.4921 -6.2 5.5 2002.4949 -2.6 6.5 2002.4976 -2.8 5.2 2002.5003 3.5 4.7 2002.5031 -1.9 5.0 2002.5058 5.8 5.0 2002.5086 -0.7 4.7 2002.5113 -1.8 5.0 2003.4805 -2.3 5.1 2003.4833 -2.7 5.1 2003.4860 -4.3 4.6 2003.4887 -1.7 5.0 2003.4915 -3.4 5.0 2003.4942 1.3 5.0 2003.4969 9.2 5.0 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> dch1.ps 2001.6133 1.9 9.0 2001.6160 -0.8 5.1 2001.6187 -0.1 3.9 2001.6215 0.8 4.2 2001.6242 2.7 4.3 2001.6270 2.4 4.3 2002.4565 55.2 4.5 2002.4593 6.8 4.6 2002.4675 6.1 4.1 2002.4702 3.5 4.2 2002.4730 6.1 4.2 2002.4757 -3.4 4.5 2002.4784 -3.1 4.5 2002.4812 -3.7 4.2 2002.4839 2.4 4.3 2002.4867 -1.4 4.6 2002.4894 -1.7 5.1 2002.4921 -6.2 5.5 2002.4949 -2.6 6.5 2002.4976 -2.8 5.2 2002.5003 3.5 4.7 2002.5031 -1.9 5.0 2002.5058 5.8 5.0 2002.5086 -0.7 4.7 2002.5113 -1.8 5.0 2003.4805 -2.3 5.1 2003.4833 -2.7 5.1 2003.4860 -4.3 4.6 2003.4887 -1.7 5.0 2003.4915 -3.4 5.0 2003.4942 1.3 5.0 2003.4969 9.2 5.0 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ -K >> dch1.ps SBreak Segment0 Comp1 2001.6133 0.7299 Segment0 Comp1 2003.4969 -0.7509 Segment0 Comp1 EOT # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R2001.5/2003.5/-30/30 \ -G230/225/200 -Y-1.75 -P -K \ -Bf0.25a1/f10a20:"Up (mm)":WeSn \ -O >> dch1.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> dch1.ps 2001.6133 14.7 23.1 2001.6160 21.3 15.1 2001.6187 12.3 11.2 2001.6215 3.7 12.0 2001.6242 8.4 12.0 2001.6270 0.8 11.3 2002.4565 17.8 12.0 2002.4593 4.9 12.7 2002.4675 14.3 12.3 2002.4702 3.7 12.3 2002.4730 1.9 12.7 2002.4757 -6.9 13.6 2002.4784 -1.6 13.0 2002.4812 2.0 11.9 2002.4839 -4.9 12.0 2002.4867 0.7 12.1 2002.4894 -1.6 12.4 2002.4921 6.9 13.5 2002.4949 3.4 13.6 2002.4976 -4.6 13.1 2002.5003 11.4 12.4 2002.5031 5.5 12.8 2002.5058 -6.6 12.8 2002.5086 -1.1 12.1 2002.5113 6.0 13.1 2003.4805 -2.4 13.2 2003.4833 -12.2 12.4 2003.4860 -0.4 12.5 2003.4887 -22.7 13.0 2003.4915 -19.9 13.8 2003.4942 -22.4 13.8 2003.4969 -18.7 13.4 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> dch1.ps 2001.6133 14.7 23.1 2001.6160 21.3 15.1 2001.6187 12.3 11.2 2001.6215 3.7 12.0 2001.6242 8.4 12.0 2001.6270 0.8 11.3 2002.4565 17.8 12.0 2002.4593 4.9 12.7 2002.4675 14.3 12.3 2002.4702 3.7 12.3 2002.4730 1.9 12.7 2002.4757 -6.9 13.6 2002.4784 -1.6 13.0 2002.4812 2.0 11.9 2002.4839 -4.9 12.0 2002.4867 0.7 12.1 2002.4894 -1.6 12.4 2002.4921 6.9 13.5 2002.4949 3.4 13.6 2002.4976 -4.6 13.1 2002.5003 11.4 12.4 2002.5031 5.5 12.8 2002.5058 -6.6 12.8 2002.5086 -1.1 12.1 2002.5113 6.0 13.1 2003.4805 -2.4 13.2 2003.4833 -12.2 12.4 2003.4860 -0.4 12.5 2003.4887 -22.7 13.0 2003.4915 -19.9 13.8 2003.4942 -22.4 13.8 2003.4969 -18.7 13.4 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ >> dch1.ps SBreak Segment0 Comp2 2001.6133 12.1402 Segment0 Comp2 2003.4969 -12.6011 Segment0 Comp2 EOT