ScatteringAtlas

Scattering Radiation Pattern Atlas: What anisotropic elastic properties can body waves resolve? Kazei & Alkhalifah, 2019 JGR - Solid Earth, https://doi.org/10.1029/2018JB016687

Reproducing images from the paper in MATLAB

cd ScatteringAtlas/scripts
mainCIJ

Figures will appear in

ScatteringAtlas/latex/Fig

Some functions can be useful as standalone

drawPatternsIJ(WT,i,j,path_pattern_save);

shows scattering of WT type (e.g. PP, PSV, SVSH etc.) on perturbation of C_ij element of stiffness matrix and saves eps to path_pattern_save.

drawPattern(WT,Cij,denFlag,patternType)

draws the radiation pattern of given type in 3D view (parmetrized by the vector g pointing towards geophone) for preliminarily assembled Cij matrix. Figure 3 is created with these finctions:

simpleRes(Cij, WT, KzMin, KzMax, dKz, phiMax, dPhi, denFlag)

computes the resolution in wavenumber(K_z)-azimuth domain (Figure 7).

resFunc(parSET)

quantifies resolution for all elastic parameters and finds the null-space for given scattering type

Main paper file

ScatteringAtlas/latex/CIJ_JGR.tex

[email protected]