GEOmap (version 2.5-11)

getmagsize: Earthquake Magnitude based on exponentional

Description

Estimate a size for plotting earthqukes recorded as a logarithmic scale

Usage

getmagsize(mag, minsize = 1, slope = 1, minmag = 0, maxmag = 8, style = 1)

Value

vector of sizes for plotting

Arguments

mag

magnitudes from catalog

minsize

minimum size

slope

slope for linear scale

minmag

min magnitude

maxmag

max magnitude

style

Style of plotting: 0= all the same size; 1(default): exponential scale; 2=linear scale

Author

Jonathan M. Lees<[e-mail]>

Details

The idea is to have a scale reflect the size of the earthquake. The default style (1) has a few parameters left over from old program geotouch.

Examples

mag = 0:9

x = runif(10, 1, 100)
y  = runif(10, 1, 100)

g = getmagsize(mag)

plot(c(0, 100), c(0, 100), asp=1, type='n')

points(x, y, pch=1, cex=g)