dTheta: Rotations in S1

In abarbour/stress: Tools for Working with Stresses and Focal-mechanism Inversions

Description

This calculates the change in orientation of the maximum principal stress relative to the fault prior to the earthquake (theta) for a given stress drop (dTau) relative to pre-mainshock levels (Tau).

Usage

dTheta(theta.deg = 30, dTau = 1, Tau = 1)

Arguments

theta.deg	numeric;
dTau	numeric;
Tau	numeric;

References

Hardebeck and Hauksson (2001): https://doi.org/10.1029/2001JB000292

Examples

thet.=seq(-90,90,by=2)
rats <- seq(-1,1,by=0.2)
All.dthet <- plyr::ldply(rats, function(r) dTheta(thet., dTau=r))
matplot(x=thet., t(All.dthet), type='l', col=c(5:1,NA,1:5), lty=c(1:5,1,5:1))

abarbour/stress documentation built on Oct. 5, 2019, 11:20 a.m.