FOCangles: Angles for focal planes

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

Description

Angles for focal planes

Usage

1

Arguments

m

moment tensor

Details

Used in MapNonDouble and doNonDouble

Value

vector of 6 angles, 3 for each plane

Note

Lower Hemisphere.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

MapNonDouble, doNonDouble, PTaxes, nodalLines

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

 di = dim(moments)
    number.of.events = di[1]
moment_11 = moments[,2]
moment_22 = moments[,3]
moment_33 = moments[,4]
moment_23 = moments[,5]
moment_13 = moments[,6]
moment_12 = moments[,7]


i = 1
m=matrix( c(moment_11[i],moment_12[i],moment_13[i],
       moment_12[i],moment_22[i],moment_23[i],
       moment_13[i],moment_23[i],moment_33[i]), ncol=3, byrow=TRUE)

   angles.all = FOCangles(m)
print(angles.all)

RFOC documentation built on May 2, 2019, 1:38 p.m.