pzoeppritz {zoeppritz}R Documentation

Plot Scattering (Zoeppritz) Coefficients

Description

Calculate and plot the P and S-wave scattering amplitudes for a plane wave at an interface.

Usage

pzoeppritz(chtype = "Amplitude", alpha1, alpha2, beta1, beta2, rho1, rho2, chincw = "P", choutkind = "ALL")

Arguments

chtype character, type of output, one of: Amplitude, Potential, Energy
alpha1 P-wave Velocity of Upper Layer, km/s
alpha2 P-wave Velocity of Lower Layer, km/s
beta1 S-wave Velocity of Upper Layer, km/s
beta2 S-wave Velocity of Lower Layer, km/s
rho1 Density of Upper Layer, kg/m3
rho2 Density of Lower Layer, kg/m3
chincw Incident Wave: P, S
choutkind character, type of out put one of: P, S, ALL, NONE

Details

Front end for zoeppritz program.

Value

List output of zoeppritz call:
angle Incident angles, degrees
rmat Matrix of 4 by n reflection coefficients for each angle
rra Matrix of 4 by n real part of scattering matrix
rra Matrix of 4 by n imaginary part of scattering matrix
ang Matrix of 4 by n phase angle
incw integer, from input parameter
icoef integer, from input parameter
alphacrit critical angle

Note

This front end is easier to call because it is more verbose. Creates a plot of the coefficients versus incident angle. If coefficients are complex, they are replaced with NA and they are thus not plotted.

Author(s)

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

References

Young,~G.B., Braile, L. W. 1976. A computer program for the application of Zoeppritz's amplitude equations and Knott's energy equations, Bulletin of the Seismological Society of America, vol.66, no.6,1881-1885.

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

See Also

zoeppritz, pzoeppritz, piczoeppritz

Examples


#############  Incident wave in Low velocity layer

alpha1 = 4.98
  beta1 =  2.9
   rho1 = 2.667

  alpha2 = 8.0
   beta2 = 4.6
   rho2 = 3.38

  App =  pzoeppritz( "Amplitude" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "P", "ALL");

App =  pzoeppritz( "Amplitude" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "S", "ALL");

#############  Incident wave in high velocity layer
  alpha1 = 8.0
  beta1 =  4.6
   rho1 = 3.38

  alpha2 = 4.98
   beta2 = 2.9
   rho2 = 2.667

 App =  pzoeppritz( "Amplitude" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "P", "ALL");

App =  pzoeppritz( "Amplitude" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "S", "ALL");



[Package zoeppritz version 1.0-2 Index]