Moment Tensor Source Type
Given a vector of EigenValues, extract the source type.
SourceType(v)
v |
vector of decreasing eigenvalues |
plotting for -30 to 30 degree quadrant.
list:
phi |
latitude angle in degrees |
lam |
longitude angle in degrees |
Jonathan M. Lees<jonathan.lees@unc.edu>
Tape, W.,and C.Tape(2012), A geometric comparison of source-type plots for moment tensors, Geophys. J. Int., 190, 499-510.
HAMMERprojXY, TapeBase, TapePlot
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" )Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.