R: Angles for focal planes

FOCangles {RFOC} | R Documentation

## Angles for focal planes

### Description

Angles for focal planes

### Usage

    FOCangles(m)

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

    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)

[Package RFOC version 3.4-3 Index]