Function takes 3x3 or vector of 6 independent elements of a moment tensor and plots the appropriate beachball for that earthquake
% Ryukyu Islands earthquake 4/18/2008 (from globalcmt.org)
M = [4.070 -1.830 -2.240 1.500 1.410 -1.110];
Mw = 5.0;
% Mt. Carmel aftershock 4/18/2008
M = [ -0.754 2.250 -1.490 -0.435 -0.094 0.730 ];
Mw = 4.6;
James Conder (2023). focalmech(fm, centerX, centerY, diam, varargin) (https://www.mathworks.com/matlabcentral/fileexchange/61227-focalmech-fm-centerx-centery-diam-varargin), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!
Added option for moment tensor to be in xyz coordinate system (default still rtf).
... and another.
Another minor tweak for efficiency.
Improved speed by vectorizing main loop. Runs 4-10x faster than previous version.
Added optional input for text above beachball. E.g.,
Fix a plotting issue that arose when matlab breaks a single contour into multiple pieces.
Fixed bug that sometimes reversed color pattern for events with >40% isotropic components.
Rewrote the portion to fill the beachballs that should behave for older versions as well as the current version of MatLab (Grzegorz's comment).