ShadowCLVD {RFOC}

Plot CLVD focal mechanism

Description

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

Usage

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<[e-mail]>

See Also

doNonDouble, MapNonDouble

Examples

############  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')