SourceType: Moment Tensor Source Type

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

View source: R/SourceType.R

Description

Given a vector of EigenValues, extract the source type.

Usage

1

Arguments

v

vector of decreasing eigenvalues

Details

plotting for -30 to 30 degree quadrant.

Value

list:

phi

latitude angle in degrees

lam

longitude angle in degrees

Author(s)

Jonathan M. Lees<[email protected]>

References

Tape, W.,and C.Tape(2012), A geometric comparison of source-type plots for moment tensors, Geophys. J. Int., 190, 499-510.

See Also

HAMMERprojXY, TapeBase, TapePlot

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
SourceType(c(1,-1,1) )

T1 = TapeBase()

m1 = list(Mxx=1.543,  Mxy=0.786,  Myy=0.336, Mxz=-2.441,  Myz=0.353,  Mzz=0.961)

i = 1
M1=matrix( c(m1$Mxx[i],m1$Mxy[i],m1$Mxz[i],
      m1$Mxy[i],m1$Myy[i],m1$Myz[i],
       m1$Mxz[i],m1$Myz[i],m1$Mzz[i]), ncol=3, byrow=TRUE)

 E1 = eigen(M1)
           h = SourceType( sort(E1$values, decreasing=TRUE) )
           h$dip = 90-h$phi
           ##  cat(paste(h$dip, h$lam, sep=" "), sep="\n")
           h1 = HAMMERprojXY(h$dip*pi/180, h$lam*pi/180)



TapePlot(T1)
           points(h1$x, h1$y,  pch=21, bg="red" )

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