#!/bin/bash ############################################################# # # # Conor McKernon - July 2003 # # # ############################################################# #### Set gmt defaults ####################################### gmtset PAPER_MEDIA a0+ PAGE_ORIENTATION landscape gmtset ANOT_FONT_SIZE 20 ANOT_OFFSET 0.3c gmtset DEGREE_FORMAT 0 gmtset GRID_PEN 0.5p/255/255/255 FRAME_PEN 3p/0/0/0 gmtset TICK_LENGTH 0.5c TICK_PEN 1p/0/0/0 touch temp.eps temp.grd temp_north.grd temp_south.grd rm temp.eps #### Conditionals to apply correct limits for grd2cpt ####### if [ "$1" == "globe" ]; then lower="-10000" upper="10000" elif [ "$1" == "relief" ]; then lower="-8000" upper="8000" elif [ "$1" == "sealand" ]; then lower="-6000" upper="3000" elif [ "$1" == "topo" ]; then lower="-7000" upper="7000" else echo "Give colour table as argument:-" echo "Either globe, relief, sealand or topo"; exit fi echo " Using colour pallete "$1 #### Sample GEBCO 1 minute topography .grd file and create ## #### corresponding .grd shading file and .cpt shading file ## echo "###### Sampling Topography grid ######" grdsample whole_world.grd -Gtemp.grd -I5m -Lx -Q -V grdinfo temp.grd grd2cpt temp.grd -C$1 -S$lower/$upper/10 -V > temp.cpt grdgradient temp.grd -Gintensity.grd -A0/270 -Ne0.6 -V #### Setup basemap for main map ############################# echo "###############################" echo "###### Plotting Main Map ######" echo "###############################" gmtset FRAME_PEN 0p/255 GRID_PEN 0p/255 TICK_PEN 0p/255 psbasemap -JKf180/109 -R0/360/-90/90 -Ba30f10g30/a15f5g15 \ -K -V -Y24.5 -X5 -U/-1.0/-22.0 \ > temp.eps gmtset FRAME_PEN 3p/0/0/0 GRID_PEN 0.25p/255/255/255 TICK_PEN 1p/0/0/0 grdimage 5min.grd \ -I5min_intensity.grd -C$1_5min_5met.cpt -R -JK -K -O -V \ >> temp.eps pscoast -JK -R -Ba30f10g30/a15f5g15 -W0.5p/0/0/0 -Di \ -A100 -N1/0.5p/0 -K -O -V \ -G150/255/100 -S100/150/255 \ >> temp.eps #### psxy calls for flinn-engdahl region boundary files ##### echo "###### Plotting F-E Regions 1 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JK -R -N -W2p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting F-E Regions 2 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JK -R -N -W1.0p/255/255/0 -K -O \ >> temp.eps done #### psxy calls for earthquake hypocentres ################## echo "###### Plotting hypocentres ######" psxy deep.dat \ -JK -R -B -Sc -K -O -G0/0/255 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy intermediate.dat \ -JK -R -B -Sc -K -O -G0/255/0 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy shallow.dat \ -JK -R -B -Sc -K -O -G255/0/0 -: -V -W0.4p/0/0/0 \ >> temp.eps #### psxy calls for plate boundary files #################### echo "###### Plotting plate boundaries 1 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JK -R -H1 -N -W3p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting plate boundaries 2 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JK -R -H1 -N -W2p/255/25/25 -K -O \ >> temp.eps done #### psxy call for euler poles ############################## echo "###### Plotting euler poles ######" psxy euler_poles.dat \ -JK -R -H1 -N -Sc0.3c -G0/0/0 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JK -R -H1 -N -Sc0.15c -G255/255/255 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JK -R -H1 -N -Sc0.6c -W2.0p/0/0/0 -K -O -: \ >> temp.eps pstext euler_poles.txt \ -JK -R -B -H1 -N -G255/200/0 -S2.0p/0/0/0 \ -K -O -V -: -D0/0.4 \ >> temp.eps ############################################################# #### Set up new basemap coordinate system for additional #### #### items (polar plots, topography scale and hypocentre #### #### scales (for depth and size)) plus data source info #### ############################################################# echo "###################################################" echo "###### Making new basemap for reference grid ######" echo "###################################################" gmtset FRAME_PEN 0p/255 GRID_PEN 0p/255 TICK_PEN 0p/255 psbasemap -Y-23.0 -X-3.0 \ -Jx1 -R0/115q/0/81 -Bf200s -K -O -V \ >> temp.eps gmtset FRAME_PEN 3p/0/0/0 GRID_PEN 0.25p/255/255/255 TICK_PEN 1p/0/0/0 #### Plot topography scale ################################## echo "###### Plotting topography scale ######" gmtset ANOT_FONT_SIZE 15 psscale -D57.5/21.0/45/1.5h -Ba1000:"Elevation":/:"m": -I \ -C$1_5min_5met.cpt -K -O -V \ >> temp.eps gmtset ANOT_FONT_SIZE 20 #### Plot title ############################################# echo "###### Plotting title ######" echo "57.7 14.5 75 0 1 CM GLOBAL SEISMICITY AND TECTONIC BOUNDARIES " |\ pstext -Jx -R -B -N -G0/0/0 -K -O -V -N \ >> temp.eps #### Plot data sources ###################################### echo "###### Plotting data sources information ######" echo "45.00 5.9 70.00 5.9 70.00 1.0 45.00 1.0" | \ psxy -Jx -R -B -A -L -N -K -O -W4p/0/0/0 -G255/255/255 \ -V \ >> temp.eps echo "45.00 5.9 70.00 5.9 70.00 1.0 45.00 1.0" | \ psxy -Jx -R -B -A -L -N -K -O -W2p/255/255/255 -V \ >> temp.eps echo "57.5 5.2 18 0 1 CM DATA SOURCES 57.5 4.5 10 0 0 CM Topography from GTOPO30 30-second grid (resampled at 5-minutes). 57.5 4.0 10 0 0 CM Bathymetry from GEBCO 1-Minute General Bathymetric Chart of the Oceans. 57.5 3.5 10 0 0 CM Earthquake hypocentres from Harvard Global CMT Catalogue (1977-2001). 57.5 3.0 10 0 0 CM Plate Boundaries from \"An updated digital model of plate boundaries\"(Bird, P., 2003, Geochemistry Geophysics Geosystems). 57.5 2.5 10 0 0 CM Euler poles from NUVEL-1A, relative to Pacific Plate (DeMets et al, 1994, Geophys. Res. Lett.). 57.5 2.0 10 0 0 CM Coastlines from World Vector Shorelines (WVS) and CIA World Data Bank II (WDBII). 57.5 1.5 10 0 0 CM Robinson Projection from WGS-84 Ellipsoid (Standard Parallels 38 N and 38 S)." | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V -N \ >> temp.eps #### Reset gmt defaults for smaller plots ################### gmtset ANOT_FONT_SIZE 15 ANOT_OFFSET 0.15c gmtset FRAME_PEN 3p/0/0/0 gmtset GRID_PEN 0.25p/255/255/255 gmtset TICK_LENGTH 0.4c TICK_PEN 1p/0/0/0 ############################################################# #### Plot north polar plot ################################## ############################################################# echo "#######################################" echo "###### Plotting north polar plot ######" echo "#######################################" psbasemap -JF0/90/30/20c -R0/360/-90/90 -B -K -O -V -X2.5 -Y2.5\ >> temp.eps grdcut 5min.grd -Gtemp_north.grd -R0/360/60/90 -V grdimage 5min_north.grd \ -I5min_intensity_north.grd -C$1_5min_5met.cpt -R -JF -K -O -V \ >> temp.eps pscoast -JF -R -Ba30g30f5/g15 -Di -A100 -W0.5p/0/0/0 \ -N1/0.5p/0 -K -O -V \ -G150/255/100 -S100/150/255 \ >> temp.eps #### psxy calls for flinn-engdahl region boundary files ##### echo "###### Plotting F-E Regions 1 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JF -R -N -W2p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting F-E Regions 2 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JF -R -N -W1.0p/255/255/0 -K -O \ >> temp.eps done #### psxy calls for earthquake hypocentres ################## echo "###### Plotting hypocentres ######" psxy deep.dat \ -JF -R -B -Sc -K -O -G0/0/255 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy intermediate.dat \ -JF -R -B -Sc -K -O -G0/255/0 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy shallow.dat \ -JF -R -B -Sc -K -O -G255/0/0 -: -V -W0.4p/0/0/0 \ >> temp.eps #### psxy calls for plate boundary files #################### echo "###### Plotting plate boundaries 1 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JF -R -H1 -W3p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting plate boundaries 2 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JF -R -H1 -W2p/255/25/25 -K -O \ >> temp.eps done #### Plot euler poles with psxy ############################# echo "###### Plotting euler poles ######" psxy euler_poles.dat \ -JF -R -H1 -Sc0.3c -G0/0/0 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JF -R -H1 -Sc0.15c -G255/255/255 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JF -R -H1 -Sc0.6c -W2.0p/0/0/0 -K -O -: \ >> temp.eps pstext euler_poles.txt \ -JF -R -B -H1 -G255/200/0 -S2.0p/0/0/0 \ -K -O -V -: -D0/0.4 \ >> temp.eps ############################################################# #### Plot south polar plot ################################## ############################################################# echo "#######################################" echo "###### Plotting south polar plot ######" echo "#######################################" psbasemap -JF0/-90/30/20c -R0/360/-90/90 -B -K -O -V -X90.0 -Y0.0 \ >> temp.eps grdcut 5min.grd -Gtemp_south.grd -R0/360/60/90 -V grdimage 5min_south.grd \ -I5min_intensity_south.grd -C$1_5min_5met.cpt -R -JF -K -O -V \ >> temp.eps pscoast -JF -R -Ba30g30f5/g15 -Di -A100 -W0.5p/0/0/0 \ -N1/0.5p/0 -K -O -V \ -G150/255/100 -S100/150/255 \ >> temp.eps #### psxy calls for flinn-engdahl region boundary files ##### echo "###### Plotting F-E Regions 2 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JF -R -N -W2p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting F-E Regions 1 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JF -R -N -W1.0p/255/255/0 -K -O \ >> temp.eps done #### psxy calls for earthquake hypocentres ################## echo "###### Plotting hypocentres ######" psxy deep.dat \ -JF -R -Sc -K -O -G0/0/255 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy intermediate.dat \ -JF -R -Sc -K -O -G0/255/0 -: -V -W0.4p/0/0/0 \ >> temp.eps psxy shallow.dat \ -JF -R -Sc -K -O -G255/0/0 -: -V -W0.4p/0/0/0 \ >> temp.eps #### psxy calls for plate boundary files #################### echo "###### Plotting plate boundaries 1 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JF -R -H1 -W3p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting plate boundaries 2 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JF -R -H1 -W2p/255/25/25 -K -O \ >> temp.eps done #### Plotting euler poles with psxy ######################### echo "###### Plotting euler poles ######" psxy euler_poles.dat \ -JF -R -H1 -Sc0.3c -G0/0/0 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JF -R -H1 -Sc0.15c -G255/255/255 -K -O -: \ >> temp.eps psxy euler_poles.dat \ -JF -R -H1 -Sc0.6c -W2.0p/0/0/0 -K -O -: \ >> temp.eps pstext euler_poles.txt \ -JF -R -B -H1 -G255/200/0 -S2.0p/0/0/0 \ -K -O -V -: -D0/0.4 \ >> temp.eps #### Reset gmt defaults for smaller plots ################### gmtset ANOT_FONT_SIZE 10 ANOT_OFFSET 0.15c gmtset FRAME_PEN 2p/0/0/0 gmtset GRID_PEN 0.25p/255/255/255 gmtset TICK_LENGTH 0.2c TICK_PEN 1p/0/0/0 ############################################################# #### Plot small plate boundaries map ######################## ############################################################# echo "###########################################" echo "###### Plotting plate boundaries map ######" echo "###########################################" pscoast -JKf180/15 -R0/360/-90/90 -Ba60f30g60/a30f30g30 -W0.25p/0/0/0 \ -Dh -G150/255/100 -S100/150/255 -A500 -K -O -V -X-68 -Y1.0 \ >> temp.eps #### psxy calls for plate boundary files #################### echo "###### Plotting plate boundaries 1 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JK -R -N -H1 -W2p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting plate boundaries 2 ######" for j in af am an ap ar as at au bh br bs bu ca cl co cr ea eu \ ft gp id jf jz ke ma mn mo ms na nb nd ni nz ok on pa \ pm ps ri sa sb sc sl so ss su sw ti to wl ya do psxy ~/PLATE_BOUNDARIES/2003_G3_PAPER/INDIVIDUAL_PLATE_BOUNDARY_FILES/${j}.dat \ -JK -R -N -H1 -W1p/255/0/0 -K -O \ >> temp.eps done ############################################################# #### Plot small Flinn-Engdahl boundaries map ################ ############################################################# echo "################################################" echo "###### Plotting Flinn-Engdahl regions map ######" echo "################################################" pscoast -JKf180/15 -R0/360/-90/90 -Ba60f30g60/a30f30g30 -W0.25p/0/0/0 \ -Dh -G150/255/100 -S100/150/255 -A500 -K -O -V -X51.0 -Y0.0 \ >> temp.eps #### psxy calls for flinn-engdahl region boundary files ##### echo "###### Plotting F-E Regions 1 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JK -R -N -W2p/0/0/0 -K -O \ >> temp.eps done echo "###### Plotting F-E Regions 2 ######" for j in 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 \ 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 \ 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 \ 46 47 48 49 50 do psxy ~/FLINN-ENGDAHL/ADJUSTED_SEISMIC_BOUNDARIES/${j}.dat \ -JK -R -N -W1p/255/255/0 -K -O \ >> temp.eps done #### Set gmt defaults ####################################### gmtset ANOT_FONT_SIZE 20 ANOT_OFFSET 0.3c gmtset FRAME_PEN 1p/0/0/0 gmtset TICK_LENGTH 0.5c TICK_PEN 1p/0/0/0 #### Another psbasemap call ################################# gmtset FRAME_PEN 0p/255 GRID_PEN 0p/255 TICK_PEN 0p/255 psbasemap -Jx1 -Y-21.5 -X-3.0 -R0/115q/0/81 -K -O -V -Bf200s \ -Y-3.5 -X-75.5 \ >> temp.eps gmtset FRAME_PEN 3p/0/0/0 GRID_PEN 0.25p/255/255/255 TICK_PEN 1p/0/0/0 #### Plotting additional map labels ######################### echo "###### Plotting plate labels ######" pstext plate_labels.txt \ -Jx -R -B -N -H1 -G255/255/255 -S2.75p/0/0/0 -K -O -V \ >> temp.eps echo "###### Plotting ocean labels ######" pstext ocean_names.txt \ -Jx -R -B -N -H1 -G255/255/255 -S1.0p/0/0/0 -K -O -V \ >> temp.eps echo "###### Plotting place labels ######" pstext place_names.txt \ -Jx -R -B -N -H1 -G200/200/200 -S1.0p/0/0/0 -K -O -V \ >> temp.eps echo "12.5 24.0 12 0 1 CM North Pole (Gnomonic Azimuthal Projection) 102.5 24.0 12 0 1 CM South Pole (Gnomonic Azimuthal Projection) 32.0 2.2 12 0 1 CM Plate Boundaries 83.0 2.2 12 0 1 CM Flinn-Engdahl Seismic Regions" | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps echo "###### Plotting labels key ######" pstext plate_labels_key.txt \ -Jx -R -B -N -H1 -G0/0/0 -K -O -V \ >> temp.eps #### Plotting key ########################################### echo "###### Plotting Key Outline ######" echo "42.50 12.5 72.50 12.5 72.50 6.5 42.50 6.5" | \ psxy -Jx -R -B -A -L -N -K -O -W4p/0/0/0 -G255/255/255 \ -V \ >> temp.eps #### Plotting plate boundary example ######################## echo "###### Plotting extra key features ######" echo "44.0 7.5 47.0 9.5 48.0 8.0" | \ psxy -Jx -R -N -L -W4p/0/0/0 -K -O \ >> temp.eps echo "44.0 7.5 47.0 9.5 48.0 8.0" | \ psxy -Jx -R -N -L -W2p/255/0/0 -K -O \ >> temp.eps echo "46.0 10.5 14 0 0 CM Plate Boundary" | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps #### Plotting FE Region boundary example #################### echo "49.5 7.5 52.5 9.5 53.5 8.0" | \ psxy -Jx -R -N -L -W4p/0/0/0 -K -O \ >> temp.eps echo "49.5 7.5 52.5 9.5 53.5 8.0" | \ psxy -Jx -R -N -L -W2p/255/255/0 -K -O \ >> temp.eps echo "52.0 11.1 14 0 0 CM Flinn-Engdahl 52.0 10.5 14 0 0 CM Seismic Region Boundary" | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps #### Plotting Euler rotation pole example ################### echo "59.0 9.0" | \ psxy -Jx -R -Sc0.3c -G0/0/0 -K -O \ >> temp.eps echo "59.0 9.0" | \ psxy -Jx -R -Sc0.15c -G255/255/255 -K -O \ >> temp.eps echo "59.0 9.0" | \ psxy -Jx -R -Sc0.6c -W2.0p/0/0/0 -K -O \ >> temp.eps echo "59.0 9.0 14 0 0 CM EU" | \ pstext -Jx -R -G255/200/0 -S2.0p/0/0/0 \ -K -O -V -D0/0.8 \ >> temp.eps echo "59.0 11.0 14 0 0 CM Euler Pole 59.0 7.9 9 0 0 CM Poles expressed in Pacific-plate reference frame. 59.0 7.5 9 0 0 CM Rotation about each pole is anticlockwise 59.0 7.1 9 0 0 CM as seen from outside the Earth." | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps #### Plotting hypocentral depth and magnitude key ########### echo "###### Plotting deep hypocentres key ######" echo "66.0 7.3 0.15 67.0 7.3 0.18 68.0 7.3 0.21 69.0 7.3 0.24 70.0 7.3 0.27 71.0 7.3 0.30" | \ psxy -Jx -R -B -Sc -N -K -O -G0/0/255 -V -W0.4p/0/0/0 \ >> temp.eps echo "###### Plotting intermediate hypocentres key ######" echo "66.0 8.1 0.15 67.0 8.1 0.18 68.0 8.1 0.21 69.0 8.1 0.24 70.0 8.1 0.27 71.0 8.1 0.30" | \ psxy -Jx -R -B -Sc -N -K -O -G0/255/0 -V -W0.4p/0/0/0 \ >> temp.eps echo "###### Plotting shallow hypocentres key ######" echo "66.0 8.9 0.15 67.0 8.9 0.18 68.0 8.9 0.21 69.0 8.9 0.24 70.0 8.9 0.27 71.0 8.9 0.30" | \ psxy -Jx -R -B -Sc -K -O -G255/0/0 -V -W0.4p/0/0/0 \ >> temp.eps echo "65.5 9.7 14 0 0 RM Depth 65.5 8.9 10 0 0 RM 0 - 70 km 65.5 8.1 10 0 0 RM 70 - 300 km 65.5 7.3 10 0 0 RM 300 - 700 km" | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps echo "68.5 11.1 14 0 0 CM Moment Magnitude (Mw) 66.0 9.5 10 45 0 LM < 5.5 67.0 9.5 10 45 0 LM 5.5 - 6.0 68.0 9.5 10 45 0 LM 6.0 - 6.5 69.0 9.5 10 45 0 LM 6.5 - 7.0 70.0 9.5 10 45 0 LM 7.0 - 7.5 71.0 9.5 10 45 0 LM > 7.5" | \ pstext -Jx -R -B -N -G0/0/0 -K -O -V \ >> temp.eps #### Plotting key title ##################################### echo "###### Plotting Key ######" echo "42.5 12.5 42.5 11.5 44.5 11.5 44.5 12.5" | \ psxy -Jx -R -B -L -K -O -G0/0/0 -V \ >> temp.eps echo "43.5 12.0 20 0 1 CM KEY" | \ pstext -Jx -R -B -N -G255/255/255 -O -V \ >> temp.eps gmtdefaults -D > .gmtdefaults ghostview -landscape -media a0 temp.eps &