Radiation pattern for SV waves
calculate the radiation patterns for SV waves
radSV(del, phiS, lam, ichi, phi)
del |
degrees, angle |
phiS |
degrees,angle |
lam |
degrees, angle |
ichi |
degrees, take off angle |
phi |
degrees, take off azimuth |
Given a focal mechanism strike-dip-rake and a given incident angle (take-off angle) and azimuth, return the SV amplitude
Amplitude of the SV wave
Jonathan M. Lees <jonathan.lees@unc.edu>
K.~Aki and P.~G. Richards.Quantitative seismology. University Science Books, Sausalito, Calif., 2nd edition, 2002.
radP, radSH, imageSV
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 = radSV(del, phiS, lam, dip, p) ### plot values image(x,y,G, asp=1)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.