Skip to contents

Calculate Principal Fault planes from stress vectors and friction

Usage

principal_fault(s1, s3, friction = 0.6)

Arguments

s1, s3

Principal stress vectors as spherical objects

friction

numeric. Coefficient of friction (0.6 by default)

Value

"Fault" object

Examples

s <- reduced_stress(angelier1990$TYM)
stress <- sigma2stress(s)
pr_TYM <- principal_fault(stress$axes[1, ], stress$axes[3, ])

stereoplot()
fault_plot(angelier1990$TYM, col = "grey")
fault_plot(pr_TYM, col = "red")
points(stress$axes, pch = 16)