superimposed_shmax_PB {tectonicr} | R Documentation |
SHmax direction resulting from multiple plate boundaries considering distance to plate boundaries
Description
Calculates a \sigma_{Hmax}
direction at given coordinates,
sourced by multiple plate boundaries. This first-order approximation is the
circular mean of the superimposed theoretical directions, weighted by the
rotation rates of the underlying PoRs, the inverse distance to the plate
boundaries, and the type of plate boundary.
Usage
superimposed_shmax_PB(
x,
pbs,
model,
rotation_weighting = TRUE,
type_weights = c(divergent = 1, convergent = 3, transform_L = 2, transform_R = 2),
idp = 1
)
Arguments
x |
grid. An object of |
pbs |
plate boundaries. |
model |
|
rotation_weighting |
logical. |
type_weights |
named vector. |
idp |
numeric. Weighting power of inverse distance. The higher the
number, the less impact far-distant boundaries have. When set to |
Value
two-column matrix. azi
is the resultant azimuth (in degrees), R
is the resultant length.
See Also
Examples
na_grid <- sf::st_make_grid(san_andreas, what = "centers", cellsize = 1)
na_plate <- subset(plates, plateA == "na" | plateB == "na")
cpm <- subset(cpm_models, cpm_models$model == "NNR-MORVEL56")
# make divergent to ridge-push:
na_plate <- transform(na_plate, type = ifelse(na_plate$pair == "eu-na", "convergent", type))
res <- superimposed_shmax_PB(na_grid, na_plate, model = cpm, idp = 2)
head(res)
[Package tectonicr version 0.4.4 Index]