pygmm.abrahamson_gregor_addo_2016.AbrahamsonGregorAddo2016

class pygmm.abrahamson_gregor_addo_2016.AbrahamsonGregorAddo2016(scenario, adjust_c1=None, adjust_c4=0, scale_atten=1.0)

Abrahamson, Gregor, and Addo (2016) ground motion moodel.

This model was developed for subduction regions and is commonly referred to as the BCHydro model.

Parameters:

scenario (pygmm.model.Scenario) – earthquake scenario

NAME = 'Abrahamson, Gregor, & Addo (2016)'

Long name of the model

ABBREV = 'AGA16'

Short name of the model

V_REF = 1000.0

COEFF = rec.array([( 0.   ,  865.1, -1.186, 1.18, 1.88, 10.,  4.2203, -1.35, 0.1, 0.9, 0., -0.0012,  1.0988, -1.42, 0.4, 3.12,  0.013 ,  0.98 , -0.0135, -0.4 , 0.9969, -1.  , 0.6, 0.43, 0.738),            ( 0.02 ,  865.1, -1.186, 1.18, 1.88, 10.,  4.2203, -1.35, 0.1, 0.9, 0., -0.0012,  1.0988, -1.42, 0.4, 3.12,  0.013 ,  0.98 , -0.0135, -0.4 , 0.9969, -1.  , 0.6, 0.43, 0.738),            ( 0.05 , 1053.5, -1.346, 1.18, 1.88, 10.,  4.5371, -1.4 , 0.1, 0.9, 0., -0.0012,  1.2536, -1.65, 0.4, 3.37,  0.013 ,  1.288, -0.0138, -0.4 , 1.103 , -1.18, 0.6, 0.43, 0.738),            ( 0.075, 1085.7, -1.471, 1.18, 1.88, 10.,  5.0733, -1.45, 0.1, 0.9, 0., -0.0012,  1.4175, -1.8 , 0.4, 3.37,  0.013 ,  1.483, -0.0142, -0.4 , 1.2732, -1.36, 0.6, 0.43, 0.738),            ( 0.1  , 1032.5, -1.624, 1.18, 1.88, 10.,  5.2892, -1.45, 0.1, 0.9, 0., -0.0012,  1.3997, -1.8 , 0.4, 3.33,  0.013 ,  1.613, -0.0145, -0.4 , 1.3042, -1.36, 0.6, 0.43, 0.738),            ( 0.15 ,  877.6, -1.931, 1.18, 1.88, 10.,  5.4563, -1.45, 0.1, 0.9, 0., -0.0014,  1.3582, -1.69, 0.4, 3.25,  0.013 ,  1.882, -0.0153, -0.4 , 1.26  , -1.3 , 0.6, 0.43, 0.738),            ( 0.2  ,  748.2, -2.188, 1.18, 1.88, 10.,  5.2684, -1.4 , 0.1, 0.9, 0., -0.0018,  1.1648, -1.49, 0.4, 3.03,  0.0129,  2.076, -0.0162, -0.35, 1.223 , -1.25, 0.6, 0.43, 0.738),            ( 0.25 ,  654.3, -2.381, 1.18, 1.88, 10.,  5.0594, -1.35, 0.1, 0.9, 0., -0.0023,  0.994 , -1.3 , 0.4, 2.8 ,  0.0129,  2.248, -0.0172, -0.31, 1.16  , -1.17, 0.6, 0.43, 0.738),            ( 0.3  ,  587.1, -2.518, 1.18, 1.88, 10.,  4.7945, -1.28, 0.1, 0.9, 0., -0.0027,  0.8821, -1.18, 0.4, 2.59,  0.0128,  2.348, -0.0183, -0.28, 1.05  , -1.06, 0.6, 0.43, 0.738),            ( 0.4  ,  503. , -2.657, 1.18, 1.88, 10.,  4.4644, -1.18, 0.1, 0.9, 0., -0.0035,  0.7046, -0.98, 0.4, 2.2 ,  0.0127,  2.427, -0.0206, -0.23, 0.8   , -0.78, 0.6, 0.43, 0.738),            ( 0.5  ,  456.6, -2.669, 1.18, 1.88, 10.,  4.0181, -1.08, 0.1, 0.9, 0., -0.0044,  0.5799, -0.82, 0.4, 1.92,  0.0125,  2.399, -0.0231, -0.19, 0.662 , -0.62, 0.6, 0.43, 0.738),            ( 0.6  ,  430.3, -2.599, 1.18, 1.88, 10.,  3.6055, -0.99, 0.1, 0.9, 0., -0.005 ,  0.5021, -0.7 , 0.4, 1.7 ,  0.0124,  2.273, -0.0256, -0.16, 0.58  , -0.5 , 0.6, 0.43, 0.738),            ( 0.75 ,  410.5, -2.401, 1.18, 1.88, 10.,  3.2174, -0.91, 0.1, 0.9, 0., -0.0058,  0.3687, -0.54, 0.4, 1.42,  0.012 ,  1.993, -0.0296, -0.12, 0.48  , -0.34, 0.6, 0.43, 0.738),            ( 1.   ,  400. , -1.955, 1.18, 1.88, 10.,  2.7981, -0.85, 0.1, 0.9, 0., -0.0062,  0.1746, -0.34, 0.4, 1.1 ,  0.0114,  1.47 , -0.0363, -0.07, 0.33  , -0.14, 0.6, 0.43, 0.738),            ( 1.5  ,  400. , -1.025, 1.18, 1.88, 10.,  2.0123, -0.77, 0.1, 0.9, 0., -0.0064, -0.082 , -0.05, 0.4, 0.7 ,  0.01  ,  0.408, -0.0493,  0.  , 0.31  ,  0.  , 0.6, 0.43, 0.738),            ( 2.   ,  400. , -0.299, 1.18, 1.88, 10.,  1.4128, -0.71, 0.1, 0.9, 0., -0.0064, -0.2821,  0.12, 0.4, 0.7 ,  0.0085, -0.401, -0.061 ,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 2.5  ,  400. ,  0.   , 1.18, 1.88, 10.,  0.9976, -0.67, 0.1, 0.9, 0., -0.0064, -0.4108,  0.25, 0.4, 0.7 ,  0.0069, -0.723, -0.0711,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 3.   ,  400. ,  0.   , 1.18, 1.88, 10.,  0.6443, -0.64, 0.1, 0.9, 0., -0.0064, -0.4466,  0.3 , 0.4, 0.7 ,  0.0054, -0.673, -0.0798,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 4.   ,  400. ,  0.   , 1.18, 1.88, 10.,  0.0657, -0.58, 0.1, 0.9, 0., -0.0064, -0.4344,  0.3 , 0.4, 0.7 ,  0.0027, -0.627, -0.0935,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 5.   ,  400. ,  0.   , 1.18, 1.88, 10., -0.4624, -0.54, 0.1, 0.9, 0., -0.0064, -0.4368,  0.3 , 0.4, 0.7 ,  0.0005, -0.596, -0.098 ,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 6.   ,  400. ,  0.   , 1.18, 1.88, 10., -0.9809, -0.5 , 0.1, 0.9, 0., -0.0064, -0.4586,  0.3 , 0.4, 0.7 , -0.0013, -0.566, -0.098 ,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            ( 7.5  ,  400. ,  0.   , 1.18, 1.88, 10., -1.6017, -0.46, 0.1, 0.9, 0., -0.0064, -0.4433,  0.3 , 0.4, 0.7 , -0.0033, -0.528, -0.098 ,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738),            (10.   ,  400. ,  0.   , 1.18, 1.88, 10., -2.2937, -0.4 , 0.1, 0.9, 0., -0.0064, -0.4828,  0.3 , 0.4, 0.7 , -0.006 , -0.504, -0.098 ,  0.  , 0.3   ,  0.  , 0.6, 0.43, 0.738)],           dtype=[('period', '<f8'), ('v_lin', '<f8'), ('b', '<f8'), ('n', '<f8'), ('c', '<f8'), ('c_4', '<f8'), ('t_1', '<f8'), ('t_2', '<f8'), ('t_3', '<f8'), ('t_4', '<f8'), ('t_5', '<f8'), ('t_6', '<f8'), ('t_7', '<f8'), ('t_8', '<f8'), ('t_9', '<f8'), ('t_10', '<f8'), ('t_11', '<f8'), ('t_12', '<f8'), ('t_13', '<f8'), ('t_14', '<f8'), ('t_15', '<f8'), ('t_16', '<f8'), ('phi', '<f8'), ('tau', '<f8'), ('sigma', '<f8')])

PERIODS = array([ 0.   ,  0.02 ,  0.05 ,  0.075,  0.1  ,  0.15 ,  0.2  ,  0.25 ,         0.3  ,  0.4  ,  0.5  ,  0.6  ,  0.75 ,  1.   ,  1.5  ,  2.   ,         2.5  ,  3.   ,  4.   ,  5.   ,  6.   ,  7.5  , 10.   ])

Indices of the periods

INDEX_PGA = 0

Index of the peak ground acceleration

INDICES_PSA = array([ 1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17,        18, 19, 20, 21, 22])

Indices for the spectral accelerations

LIMITS = {'dist_jb': (0.0, 300.0), 'mag': (3.0, 8.5), 'v_s30': (150.0, 1500.0)}

Limits of model applicability

PARAMS = [<pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.CategoricalParameter object>, <pygmm.model.CategoricalParameter object>]

Model parameters

__init__(scenario, adjust_c1=None, adjust_c4=0, scale_atten=1.0)

Initialize the model.

Args:

scenario (pygmm.model.Scenario): earthquake scenario.

property adjust_c1

property adjust_c4

property scale_atten

INDEX_PGD = None

Index of the peak ground displacement

INDEX_PGV = None

Index of the peak ground velocity

PGD_SCALE = 1.0

Scale factor to apply to get PGD in cm

PGV_SCALE = 1.0

Scale factor to apply to get PGV in cm/sec

interp_ln_spec_accels(periods, kind='linear')

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

• periods (array_like) – spectral periods to interpolate the response.

• kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

interp_ln_stds(periods, kind='linear')

Interpolate the logarithmic standard deviation.

Interpolate the logarithmic standard deviation (\(\sigma_{\ln}\)) of spectral acceleration at the provided damping at specified periods.

Parameters:

• periods (array_like) – spectral periods to interpolate the response.

• kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_stds – interpolated logarithmic standard deviations

Return type:

np.ndarray

interp_spec_accels(periods, kind='linear')

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

• periods (array_like) – spectral periods to interpolate the response.

• kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

property ln_std_pga: float

Peak ground accelaration log-standard deviation.

property ln_std_pgd: float

Peak ground displacement log-standard deviation.

property ln_std_pgv: float

Peak ground velocity log-standard deviation.

property ln_stds: ndarray

Pseudo-spectral accelerations log-standard deviation.

property periods: ndarray

Periods specified by the model.

property pga: float

Peak ground acceleration (PGA) computed by the model (g).

property pgd: float

Peak ground displacement (PGD) computed by the model (cm).

property pgv: float

Peak ground velocity (PGV) computed by the model (cm/sec).

property scenario

property spec_accels: ndarray

Pseudo-spectral accelerations computed by the model (g).