RFOC (version 1.0-4)

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

Value

• Amplitude of the P wave

Details

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

References

K.~Aki and P.~G. Richards.{Quantitative seismology. University Science Books, Sausalito, Calif., 2nd edition, 2002.} [object Object] radP, radSV, imageP 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<="" p="">

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)

misc