Skip to content

Latest commit

 

History

History
40 lines (39 loc) · 1.29 KB

README.md

File metadata and controls

40 lines (39 loc) · 1.29 KB
Raw

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 Cij 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(Kz)-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]