radSV | R Documentation |
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.