def compute(self, ctx: np.recarray, imts, mean, sig, tau, phi):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.compute>`
for spec of input and result values.
"""
ctx.vs30 = np.full_like(ctx.vs30, 600.)
super().compute(ctx, imts, mean, sig, tau, phi)
tau_ss = 'tau'
log_phi_ss = np.log(10)
for m, imt in enumerate(imts):
_apply_adjustments(AkkarBommer2010.COEFFS,
self.COEFFS_FS_ROCK[imt], tau_ss, mean[m],
sig[m], tau[m], phi[m], ctx, ctx.rjb, imt,
log_phi_ss)
sig[m] = np.log(10**sig[m])
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
sites.vs30 = 620 * np.ones(len(sites.vs30))
mean, stddevs = super(ChiouYoungs2008SWISS01, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
log_phi_ss = 1
tau = self.get_tau(ChiouYoungs2008.COEFFS[imt], rup)
ln_y_ref = super(ChiouYoungs2008SWISS01, self).\
_get_ln_y_ref(rup, dists, ChiouYoungs2008.COEFFS[imt])
exp1 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] *
(sites.vs30.clip(-np.inf, 1130) - 360))
exp2 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] * (1130 - 360))
nl = self.get_nl(ChiouYoungs2008.COEFFS[imt], ln_y_ref, exp1, exp2)
mean, stddevs = _apply_adjustments(
ChiouYoungs2008.COEFFS, self.COEFFS_FS_ROCK[imt], 1,
mean, stddevs, sites, rup, dists.rjb, imt, stddev_types,
log_phi_ss, NL=nl, tau_value=tau)
return mean, stddevs
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
sites.vs30 = 620 * np.ones(len(sites.vs30))
mean, stddevs = super(ChiouYoungs2008SWISS01, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
log_phi_ss = 1
tau = self.get_tau(ChiouYoungs2008.COEFFS[imt], rup)
ln_y_ref = super(ChiouYoungs2008SWISS01, self).\
_get_ln_y_ref(rup, dists, ChiouYoungs2008.COEFFS[imt])
exp1 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] *
(sites.vs30.clip(-np.inf, 1130) - 360))
exp2 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] * (1130 - 360))
nl = self.get_nl(ChiouYoungs2008.COEFFS[imt], ln_y_ref, exp1, exp2)
mean, stddevs = _apply_adjustments(
ChiouYoungs2008.COEFFS, self.COEFFS_FS_ROCK[imt], 1,
mean, stddevs, sites, rup, dists.rjb, imt, stddev_types,
log_phi_ss, NL=nl, tau_value=tau)
return mean, stddevs
def compute(self, ctx, imts, mean, sig, tau, phi):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.compute>`
for spec of input and result values.
"""
ctx = copy.copy(ctx)
ctx.vs30 = 700 * np.ones(len(ctx.vs30))
super().compute(ctx, imts, mean, sig, tau, phi)
tau_ss = 'tauC'
log_phi_ss = 1.00
C = ZhaoEtAl2006AscSWISS05.COEFFS_ASC
for m, imt in enumerate(imts):
_apply_adjustments(
C, self.COEFFS_FS_ROCK[imt], tau_ss,
mean[m], sig[m], tau[m], phi[m],
ctx, ctx.rrup, imt, log_phi_ss)
def compute(self, ctx, imts, mean, sig, tau, phi):
ctx.vs30 = 620 * np.ones(len(ctx.vs30))
log_phi_ss = 1
super().compute(ctx, imts, mean, sig, tau, phi)
for m, imt in enumerate(imts):
t = get_tau(ChiouYoungs2008.COEFFS[imt], ctx)
ln_y_ref = _get_ln_y_ref(ctx, ChiouYoungs2008.COEFFS[imt])
exp1 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] *
(ctx.vs30.clip(-np.inf, 1130) - 360))
exp2 = np.exp(ChiouYoungs2008.COEFFS[imt]['phi3'] * (1130 - 360))
nl = get_nl(ChiouYoungs2008.COEFFS[imt], ln_y_ref, exp1, exp2)
_apply_adjustments(
ChiouYoungs2008.COEFFS, self.COEFFS_FS_ROCK[imt], 1,
mean[m], sig[m], tau[m], phi[m], ctx, ctx.rjb, imt,
log_phi_ss, NL=nl, tau_value=t)
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
for spec of input and result values.
"""
mean, stddevs = super(AkkarBommer2010SWISS01, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
tau_ss = 'tau'
log_phi_ss = np.log(10)
mean, stddevs = _apply_adjustments(
AkkarBommer2010.COEFFS, self.COEFFS_FS_ROCK[imt], tau_ss,
mean, stddevs, sites, rup, dists.rjb, imt, stddev_types,
log_phi_ss)
return mean, np.log(10 ** np.array(stddevs))
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
for spec of input and result values.
"""
mean, stddevs = super(ZhaoEtAl2006AscSWISS05, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
tau_ss = 'tauC'
log_phi_ss = 1.00
C = ZhaoEtAl2006AscSWISS05.COEFFS_ASC
mean, stddevs = _apply_adjustments(
C, self.COEFFS_FS_ROCK[imt], tau_ss,
mean, stddevs, sites, rup, dists.rrup, imt, stddev_types,
log_phi_ss)
return mean, stddevs
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
for spec of input and result values.
"""
sites.vs30 = 600 * np.ones(len(sites.vs30))
mean, stddevs = super(AkkarBommer2010SWISS01, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
tau_ss = 'tau'
log_phi_ss = np.log(10)
mean, stddevs = _apply_adjustments(
AkkarBommer2010.COEFFS, self.COEFFS_FS_ROCK[imt], tau_ss,
mean, stddevs, sites, rup, dists.rjb, imt, stddev_types,
log_phi_ss)
return mean, np.log(10 ** np.array(stddevs))
def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
for spec of input and result values.
"""
sites.vs30 = 700 * np.ones(len(sites.vs30))
mean, stddevs = super(ZhaoEtAl2006AscSWISS05, self).\
get_mean_and_stddevs(sites, rup, dists, imt, stddev_types)
tau_ss = 'tauC'
log_phi_ss = 1.00
C = ZhaoEtAl2006AscSWISS05.COEFFS_ASC
mean, stddevs = _apply_adjustments(C, self.COEFFS_FS_ROCK[imt], tau_ss,
mean, stddevs, sites, rup,
dists.rrup, imt, stddev_types,
log_phi_ss)
return mean, stddevs
