Learn R Programming

RFOC (version 2.0-00)

radSH: Radiation pattern for SH waves

Description

calculate the radiation patterns for SH waves

Usage

radSH(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

Value

  • Amplitude of the SH wave

Details

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

References

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

See Also

radP, radSV, imageSH

Examples

Run this code
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 = radSH(del, phiS, lam, dip, p)

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

Run the code above in your browser using DataLab