pzoeppritz {zoeppritz}

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<[e-mail]>

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");