zoeppritz-package {zoeppritz} | R Documentation |
Calculate and plot scattering matrix coefficients for plane waves at interface.
Package: | zoeppritz |
Type: | Package |
Version: | 1.0-2 |
Date: | 2007-10-15 |
License: | GPL |
Jonathan M. Lees<jonathan.lees@unc.edu> Maintainer: Jonathan M. Lees<jonathan.lees@unc.edu>
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.
##### demo(ZOEP) ############# 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"); App = pzoeppritz( "Energy" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "P", "ALL"); App = pzoeppritz( "Potential" , alpha1, alpha2, beta1, beta2, rho1 ,rho2, "P", "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");