Haz_calc: The hazard calculator

In wltcwpf/RPSHA: This package implements LASSO to conduct earthquake event/ground motion map selection

The hazard calculator

Description

This function calculates hazard for the given IMs at the user-defined return periods or IM levels at the given sites using UCERF3 3.1 and 3.2 branches average source model.

Usage

Haz_calc(
  slat,
  slon,
  Vs30,
  z1 = -999,
  max_dist = 300,
  periods,
  IM_mat = NA,
  site_name = "Site1",
  output_dir
)

Arguments

slat	The latitude of the interest site
slon	The longitude of the interest site
Vs30	The associated Vs30 value (in m/s) to the interest site
z1	The associated z1 value (in km) to the interest site. -999 for unknown.
max_dist	The maximum cut-off distance for considered sources. Default is 300 km.
periods	The list of interest IMs. Use -1 for PGV, 0 for PGA, and other NGA-West2 periods for PSA
IM_mat	A matrix that specifies the IM exceedance levels for hazard curve calculation. If periods length is 1, then IM_mat is a vector; otherwise, IM_mat is a matrix with row number equals to periods length. Each row in IM_mat is the specified IM exceedance levels for each of IM in periods. The default is NA, which uses the default IM exceedance levels.
site_name	The name of the site
output_dir	The directory for the outputs

Examples

Haz_calc(slat = 36.1, slon = -118.1, Vs30 = 350,
z1 = -999, periods = c(-1, 0.01, 5), output_dir = '~/Desktop/RPSHA')

Haz_calc(slat = 36.5, slon = -119.5, Vs30 = 567,
z1 = 1.5, periods = c(0.01, 1, 5),
IM_mat = rbind(c(0.001, 0.0016), c(0.001, 0.0016), c(0.001, 0.0016)), site_name = 'Site2',
output_dir = '~/Desktop/RPSHA')

wltcwpf/RPSHA documentation built on July 6, 2023, 1:02 a.m.