#!/bin/csh # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R1992.5/2004/-80/70 \ -G230/225/200 -Y5.25 -P -K \ -Bf0.25a1/f20a40:"North (mm)":WeSn \ > cape.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> cape.ps 1992.6521 -77.0 7.1 1992.6548 -63.2 6.6 1993.6405 -69.5 7.2 1995.5817 -40.7 6.8 1996.6221 -19.5 7.3 1996.6248 -27.4 7.0 1998.6289 5.0 6.8 1998.6317 -5.8 6.5 2003.5544 68.0 6.3 2003.5572 68.2 6.4 2003.5599 66.2 6.6 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> cape.ps 1992.6521 -77.0 7.1 1992.6548 -63.2 6.6 1993.6405 -69.5 7.2 1995.5817 -40.7 6.8 1996.6221 -19.5 7.3 1996.6248 -27.4 7.0 1998.6289 5.0 6.8 1998.6317 -5.8 6.5 2003.5544 68.0 6.3 2003.5572 68.2 6.4 2003.5599 66.2 6.6 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ -K >> cape.ps SBreak Segment0 Comp0 1992.6521 -74.6616 Segment0 Comp0 2003.5599 66.1067 Segment0 Comp0 EOT # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R1992.5/2004/-60/30 \ -G230/225/200 -Y-1.75 -P -K \ -Bf0.25a1/f20a40:"East (mm)":WeSn \ -O >> cape.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> cape.ps 1992.6521 -44.9 8.0 1992.6548 -56.9 7.2 1993.6405 -17.7 7.8 1995.5817 -19.6 7.3 1996.6221 -4.6 9.8 1996.6248 -9.0 8.7 1998.6289 0.8 8.0 1998.6317 4.8 7.5 2003.5544 29.3 7.4 2003.5572 23.2 7.4 2003.5599 22.9 7.9 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> cape.ps 1992.6521 -44.9 8.0 1992.6548 -56.9 7.2 1993.6405 -17.7 7.8 1995.5817 -19.6 7.3 1996.6221 -4.6 9.8 1996.6248 -9.0 8.7 1998.6289 0.8 8.0 1998.6317 4.8 7.5 2003.5544 29.3 7.4 2003.5572 23.2 7.4 2003.5599 22.9 7.9 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ -K >> cape.ps SBreak Segment0 Comp1 1992.6521 -36.6057 Segment0 Comp1 2003.5599 27.1145 Segment0 Comp1 EOT # Draw basemap # ------------ gmtset ANOT_FONT Helvetica LABEL_FONT Helvetica psbasemap -JX6/1.5 \ -R1992.5/2004/-100/40 \ -G230/225/200 -Y-1.75 -P -K \ -Bf0.25a1/f20a40:"Up (mm)":WeSn \ -O >> cape.ps # Draw error bars # --------------- psxy << EOT \ -JX -R -G255/0/0 -Ey0.05/0.1/159/159/159 -O \ -K >> cape.ps 1992.6521 -40.7 22.9 1992.6548 -94.2 21.9 1993.6405 -0.6 23.5 1995.5817 26.9 21.5 1996.6221 20.1 23.5 1996.6248 24.8 26.4 1998.6289 -15.8 23.6 1998.6317 -4.6 22.3 2003.5544 38.7 22.8 2003.5572 14.2 22.6 2003.5599 21.0 24.1 EOT # Draw points # ----------- psxy << EOT \ -JX -R -Sc0.06 -G255/0/0 -O \ -K >> cape.ps 1992.6521 -40.7 22.9 1992.6548 -94.2 21.9 1993.6405 -0.6 23.5 1995.5817 26.9 21.5 1996.6221 20.1 23.5 1996.6248 24.8 26.4 1998.6289 -15.8 23.6 1998.6317 -4.6 22.3 2003.5544 38.7 22.8 2003.5572 14.2 22.6 2003.5599 21.0 24.1 EOT # Draw best fit line # ------------------ psxy << EOT \ -JX -R -MS -W3/102/102/255 -O \ >> cape.ps SBreak Segment0 Comp2 1992.6521 -29.7262 Segment0 Comp2 2003.5599 28.6264 Segment0 Comp2 EOT