Skip to contents

Build a log-spaced ky grid for one scenario

Source: R/getDnKy.R
getDnKy.Rd

Derives the effective intensities PGA, Sa(1.3Ts), Sa(1.5Ts) from the scenario UHS and constructs a log-spaced ky grid covering the calibrated range of the displacement models.

Usage

getDnKy(uhs, Ts, kyN = 30L, pRef = "mean")

Arguments

uhs

data.table with columns Tn, p, Sa (one scenario)

Ts

numeric scalar, fundamental period (s)

kyN

integer, number of grid points (default 30)

pRef

character, probability level for intensity anchor (default "mean")

Value

numeric vector of ky values (log-spaced)

Details

Default bounds come from the calibration ranges of the models in the ensemble:

  • lower: 0.01g (BM17 lowest calibration point)

  • upper: max(PGA, 0.80g). The physical limit is ky < PGA (no slip below the yield acceleration); the 0.80g floor extends the grid for low-PGA scenarios up to BM17's calibration ceiling.

Examples

uhs <- data.table::data.table(Tn = 0, p = "mean", Sa = 0.5)
getDnKy(uhs, Ts = 0.3)
#>  [1] 0.01000000 0.01163118 0.01352844 0.01573517 0.01830186 0.02128722
#>  [7] 0.02475955 0.02879828 0.03349580 0.03895957 0.04531458 0.05270620
#> [13] 0.06130353 0.07130325 0.08293409 0.09646214 0.11219685 0.13049819
#> [19] 0.15178479 0.17654363 0.20534109 0.23883592 0.27779437 0.32310764
#> [25] 0.37581233 0.43711410 0.50841530 0.59134701 0.68780637 0.80000000