doNonDouble: Plot Non-double Couple Moment

In RFOC: Graphics for Spherical Distributions and Earthquake Focal Mechanisms

Plot Non-double Couple Moment

Description

Plot Non-double Couple Moment

Usage

doNonDouble(moments, sel = 1, col=rgb(1, .75, .75))

Arguments

moments	list of moments: seven elements. See details.
sel	integer vector, index of moments to plot
col	color, either a single color, rgb, or a color palette

Details

Plot, sequentially the moments using the CLVD (non-double couple component. The first element of the list is the integer index of the event. The next six elements are the moments in the following order, c(Mxx, Myy, Mzz, Mzy, Mxz, Mxy) .

If the data is in spherical coordinates, one must switch the sign of the Mrp and Mtp components, so: ⁠ Mrr = Mzz Mtt = Mxx Mpp = Myy Mrt = Mxz Mrp = -Myz Mtp = -Mxy ⁠

Value

Side effects

Note

If events are read in using spherical rather than cartesian coordinates need a conversion: ⁠ Mrr = Mzz Mtt = Mxx Mpp = Myy Mrt = Mxz Mrp = -Myz Mtp = -Mxy ⁠

Author(s)

Jonathan M. Lees<[e-mail]>

References

Ekstrom, G.; Nettles, M. and DziewoDski, A. The Global CMT Project 2004-2010: centroid-moment tensors for 13,017 earthquakes Physics of the Earth and Planetary Interiors, 2012.

See Also

MapNonDouble, ShadowCLVD, angles, nodalLines, PTaxes

Examples

mo = list(n=1, m1=1.035675e+017, m2=-1.985852e+016,
m3=-6.198052e+014, m4=1.177936e+017, m5=-7.600627e+016,
m6=-3.461405e+017)

moments = cbind(mo$n, mo$m1, mo$m2, mo$m3, mo$m4, mo$m5, mo$m6)

doNonDouble(moments)

RFOC documentation built on Sept. 8, 2023, 6:12 p.m.