Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

radP

Radiation pattern for P waves

Description

calculate the radiation patterns for P waves

Usage

radP(del, phiS, lam, ichi, phi)

Arguments

del

degrees, angle

phiS

degrees,angle

lam

degrees, angle

ichi

degrees, take off angle
phi

degrees, take off azimuth

Details

Given a focal mechanism strike-dip-rake and a given incident angle (take-off angle) and azimuth, return the P amplitude

Value

Amplitude of the P wave

Author(s)

Jonathan M. Lees <jonathan.lees@unc.edu>

References

K.~Aki and P.~G. Richards.Quantitative seismology. University Science Books, Sausalito, Calif., 2nd edition, 2002.

See Also

radP, radSV, imageP

Examples

phiS=65
del=25
lam=13
x = seq(-1, 1, 0.01)
y = x

X = matrix(rep(x, length(y)), nrow= length(x))
Y = t(X)
RAD2DEG = 180/pi
p = RAD2DEG*(pi/2 -atan2(Y, X))
p[p<0] = p[p<0] + 360

R = sqrt(X^2+Y^2)
R[R>1] = NaN
dip =RAD2DEG*2*asin(R/sqrt(2))

###  Calculate the radiation pattern
G = radP(del, phiS, lam, dip, p)

###  plot values
image(x,y,G, asp=1)
RFOC

Graphics for Spherical Distributions and Earthquake Focal Mechanisms

v3.4-6
GPL (>= 2)
Authors
Jonathan M. Lees [aut, cre], Keehoon Kim [ctb]
Initial release
2018-01-17

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.