ShadowCLVD: Plot CLVD focal mechanism

Description Usage Arguments Details Value Note Author(s) See Also Examples

View source: R/ShadowCLVD.R

Description

Plot non-double couple part of the focal mechanism provided in the moment tensor.

Usage

1
ShadowCLVD(m, PLOT = TRUE, col=rgb(1, .75, .75))

Arguments

m

moment tensor

PLOT

logical, TRUE means plot

col

color, either a single color, rgb, or a color palette

Details

This code is meant to be used with doNonDouble or MapNonDouble functions for plotting the non-double couple components of the moment tensor. A color palette can be provided for some details of the radiation patterns, e.g. col=rainbow(12).

Value

Side effects and image list

Note

Lower Hemisphere.

Author(s)

Jonathan M. Lees<[email protected]>

See Also

doNonDouble, MapNonDouble

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
############  moment tensor from Harvard CMT catalog
sponent = 26
ef = 1*10^(sponent)
Mrr =  2.375*ef
Mtt = -2.777*ef
Mpp = 0.403*ef
Mrt = 2.800*ef
Mrp = 1.190*ef
Mtp = -0.539*ef

############  convert to cartesian coordinates
Mzz=Mrr
Mxx= Mtt
Myy= Mpp
Mxz= Mrt
Myz= -Mrp 
Mxy= -Mtp


m=matrix( c(Mxx,Mxy,Mxz,
      Mxy,Myy,Myz,
       Mxz,Myz,Mzz), ncol=3, byrow=TRUE)

Fi=seq(from=0, by=0.1, to=361)
  ###  dev.new()
    plot(cos(Fi*pi/180.0),sin(Fi*pi/180.0),type='l', asp=1 , ann=FALSE, axes=FALSE)
  
  ShadowCLVD(m, col='red')

RFOC documentation built on Jan. 20, 2018, 9:46 a.m.