pygmm.derras_bard_cotton_2014.DerrasBardCotton2014

class pygmm.derras_bard_cotton_2014.DerrasBardCotton2014(scenario)[source]

Derras, Bard and Cotton (2014, [Derras et al., 2014]) model.

Parameters:

scenario (pygmm.model.Scenario) – earthquake scenario

NAME = 'Derras, Bard & Cotton (2014)'

Long name of the model

ABBREV = 'DBC13'

Short name of the model

COEFF = {'b_1': [-1.271232488, 1.512611028, 0.591089009, -0.126622688, -0.415721222], 'b_2': [-0.076472745, -0.143493078, -0.141251167, -0.1402921, -0.142446238, -0.169403545, -0.168605778, -0.154859216, -0.169019709, -0.176613627, -0.170419807, -0.169012965, -0.163928953, -0.174008354, -0.194478386, -0.200801472, -0.203141294, -0.212134795, -0.214788445, -0.214182427, -0.213932639, -0.212008293, -0.227037484, -0.246738859, -0.24795058, -0.2518272, -0.212879565, -0.215512278, -0.232808829, -0.242614348, -0.230743172, -0.219657843, -0.210397231, -0.197795308, -0.160530216, -0.140822986, -0.107838684, -0.094014862, -0.073895337, -0.05799372, -0.042284114, -0.043692208, -0.027928475, -0.005401146, 0.01643085, 0.056436981, 0.076415498, 0.080005285, 0.07403565, 0.09299877, 0.098090915, 0.094946234, 0.087002098, 0.080633932, 0.091510963, 0.09495708, 0.095220431, 0.086504663, 0.080948314, 0.069604068, 0.059710513, 0.051176846, 0.053684364, 0.056169376], 'log10_std': {'between': [0.149, 0.155, 0.155, 0.157, 0.16, 0.162, 0.163, 0.165, 0.168, 0.17, 0.17, 0.169, 0.17, 0.169, 0.168, 0.168, 0.167, 0.165, 0.164, 0.164, 0.164, 0.163, 0.161, 0.161, 0.163, 0.164, 0.165, 0.166, 0.165, 0.165, 0.164, 0.164, 0.164, 0.164, 0.165, 0.166, 0.166, 0.166, 0.167, 0.166, 0.165, 0.166, 0.167, 0.168, 0.17, 0.172, 0.174, 0.177, 0.179, 0.18, 0.181, 0.183, 0.185, 0.184, 0.186, 0.186, 0.187, 0.188, 0.189, 0.19, 0.19, 0.189, 0.188, 0.188], 'total': [0.298, 0.309, 0.31, 0.313, 0.319, 0.323, 0.325, 0.328, 0.335, 0.338, 0.338, 0.337, 0.338, 0.337, 0.335, 0.334, 0.333, 0.33, 0.328, 0.328, 0.326, 0.325, 0.322, 0.322, 0.326, 0.328, 0.329, 0.33, 0.33, 0.329, 0.327, 0.328, 0.328, 0.327, 0.33, 0.331, 0.332, 0.332, 0.332, 0.331, 0.33, 0.331, 0.333, 0.335, 0.339, 0.343, 0.346, 0.353, 0.356, 0.359, 0.362, 0.365, 0.368, 0.368, 0.37, 0.37, 0.373, 0.375, 0.377, 0.378, 0.378, 0.376, 0.375, 0.375], 'within': [0.258, 0.267, 0.268, 0.27, 0.276, 0.279, 0.281, 0.284, 0.29, 0.293, 0.292, 0.292, 0.293, 0.292, 0.29, 0.289, 0.288, 0.285, 0.284, 0.284, 0.282, 0.281, 0.279, 0.279, 0.282, 0.284, 0.285, 0.286, 0.285, 0.284, 0.283, 0.283, 0.284, 0.283, 0.285, 0.286, 0.287, 0.287, 0.287, 0.287, 0.285, 0.286, 0.288, 0.29, 0.293, 0.297, 0.299, 0.305, 0.308, 0.31, 0.313, 0.316, 0.319, 0.319, 0.32, 0.32, 0.323, 0.324, 0.326, 0.327, 0.327, 0.326, 0.325, 0.324]}, 'min_max': {'depth_hyp': [0.0, 25.0], 'log10_dist_jb': [-1.0, 2.738], 'log10_resp': [[-3.849485002, -0.060923911], [-2.979303657, 0.98101835], [-2.985196725, 0.99145166], [-2.986059277, 1.007712442], [-2.984175815, 1.08949524], [-2.982552276, 1.156662115], [-2.956637722, 1.157095399], [-2.945980456, 1.260733107], [-2.936307635, 1.410814288], [-2.932216636, 1.457727792], [-2.930531138, 1.444991459], [-2.928542596, 1.453062036], [-2.919667313, 1.385037364], [-2.911697131, 1.471769235], [-2.8705569, 1.489542018], [-2.85933557, 1.490888494], [-2.873058765, 1.486617287], [-2.850059029, 1.50374993], [-2.834673005, 1.535394915], [-2.781647452, 1.409036446], [-2.758712763, 1.261319049], [-2.767214299, 1.232880087], [-2.753912978, 1.267300559], [-2.685249329, 1.269455682], [-2.699180576, 1.248427106], [-2.657061356, 1.26168743], [-2.743779615, 1.245493493], [-2.719549792, 1.212558872], [-2.667871945, 1.206737632], [-2.619043181, 1.183978702], [-2.667081411, 1.188206731], [-2.715295253, 1.187534293], [-2.744937774, 1.195017056], [-2.791616858, 1.162867247], [-2.908517633, 1.149294497], [-2.952789244, 1.226266131], [-3.002621528, 1.174170755], [-3.090715515, 1.128224702], [-3.186100033, 1.052990942], [-3.242837714, 1.035664147], [-3.309818391, 1.006674652], [-3.347867863, 1.022366733], [-3.376363292, 0.993946952], [-3.411540737, 0.96655149], [-3.499349533, 0.88712674], [-3.59867812, 0.786452684], [-3.699971186, 0.773370902], [-3.787935025, 0.749136058], [-3.858485145, 0.780103015], [-3.919735723, 0.735773313], [-3.986529186, 0.64299905], [-4.037034701, 0.610951381], [-4.084673593, 0.59822423], [-4.123668528, 0.587317187], [-4.238777883, 0.494881369], [-4.343911322, 0.437831896], [-4.425698673, 0.419551103], [-4.470029056, 0.427769901], [-4.531743129, 0.455067137], [-4.600329725, 0.466260482], [-4.644941317, 0.475656195], [-4.694669918, 0.47363989], [-4.750844723, 0.467695635], [-4.79929973, 0.444726741]], 'log10_v_s30': [1.9638, 3.2035], 'mag': [3.6, 7.6], 'mechanism': [1.0, 4.0]}, 'period': [-1, 0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.075, 0.1, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4], 'w_1': [[2.647891635, -1.07021796, 0.174087758, 0.092191287, -0.013763679], [-1.908675436, -0.535017369, -0.705141623, 0.167611583, -0.02661049], [0.203542143, 1.780557636, -0.080494591, 0.013596356, 0.061508209], [-0.692737498, 0.441505232, 0.775579973, -0.031717733, -0.16306571], [0.016162821, 0.218141339, -1.606099447, -0.041636256, 0.026057983]], 'w_2': [[-0.510826776, 0.070554779, 0.220914175, 0.168815839, 0.170928164], [-0.54101415, 0.254251327, 0.109777617, 0.075959095, -0.020347572], [-0.539773537, 0.254301257, 0.107974037, 0.074881998, -0.021585429], [-0.539375931, 0.248942044, 0.099446671, 0.081987137, -0.017662836], [-0.537208413, 0.240441278, 0.083145677, 0.087421772, -0.018053459], [-0.520634258, 0.229422592, 0.064081989, 0.087723565, -0.026990211], [-0.536143529, 0.237932202, 0.052958747, 0.079850429, -0.040505524], [-0.554628064, 0.271755601, 0.050356632, 0.055594074, -0.076045335], [-0.532266517, 0.291446897, 0.064617059, 0.055270794, -0.09645035], [-0.521512879, 0.29658671, 0.069576025, 0.056273712, -0.099062435], [-0.5198772, 0.301189055, 0.074326652, 0.054844753, -0.10296585], [-0.514898774, 0.30375476, 0.081772384, 0.055068191, -0.102310668], [-0.516215894, 0.313439795, 0.091938911, 0.052380539, -0.102390682], [-0.4979795, 0.317018727, 0.102152195, 0.049835913, -0.09816729], [-0.493176711, 0.327007863, 0.11418295, 0.049209182, -0.094872584], [-0.485747878, 0.331253599, 0.127205614, 0.046418473, -0.091736341], [-0.478191252, 0.335381467, 0.137325066, 0.044229522, -0.088439399], [-0.472577837, 0.335355736, 0.144140467, 0.044729064, -0.083059358], [-0.460675639, 0.331629622, 0.148101401, 0.05134718, -0.074474238], [-0.455401286, 0.345362019, 0.168485452, 0.057154914, -0.06449186], [-0.45629014, 0.359907531, 0.19163738, 0.058118391, -0.055195467], [-0.448942858, 0.352037712, 0.206011412, 0.061961861, -0.044676346], [-0.437033339, 0.340535007, 0.214707559, 0.064792403, -0.031257631], [-0.449973961, 0.343550607, 0.23202436, 0.066267603, -0.021772794], [-0.435914287, 0.327845458, 0.237413306, 0.068819597, -0.011293778], [-0.42465784, 0.321408865, 0.247731596, 0.077612939, 0.002352304], [-0.415678041, 0.314753909, 0.255990159, 0.081865016, 0.015181998], [-0.415314256, 0.31326253, 0.266736228, 0.090005027, 0.027379773], [-0.414839332, 0.313509014, 0.276335858, 0.101142858, 0.0395314], [-0.418038972, 0.304164889, 0.279426675, 0.109541013, 0.055573293], [-0.412223333, 0.294185023, 0.28018961, 0.10824389, 0.065720829], [-0.405732831, 0.287156033, 0.283138956, 0.108857627, 0.074589338], [-0.404244295, 0.282654061, 0.289679834, 0.110693756, 0.082047467], [-0.398349555, 0.275134246, 0.293426945, 0.11309985, 0.091428891], [-0.382933503, 0.253450948, 0.297505456, 0.122089801, 0.115887837], [-0.373349488, 0.240256546, 0.303362656, 0.136961296, 0.139021525], [-0.373729759, 0.228657553, 0.312794291, 0.14701375, 0.160696896], [-0.359574788, 0.213839065, 0.31311436, 0.151168709, 0.17805053], [-0.345784509, 0.197187139, 0.311285387, 0.158730867, 0.193066308], [-0.341037761, 0.18153626, 0.315005438, 0.166201307, 0.206400745], [-0.336228384, 0.171930075, 0.320637744, 0.170489, 0.214822499], [-0.33108553, 0.161317428, 0.322442037, 0.17349243, 0.222366883], [-0.32867094, 0.149886402, 0.323389852, 0.176411536, 0.229013406], [-0.323017173, 0.137567256, 0.32307006, 0.180835711, 0.230941347], [-0.318080249, 0.124328971, 0.331785162, 0.188006751, 0.238824085], [-0.314220923, 0.115845087, 0.342820043, 0.196252312, 0.246883255], [-0.310585172, 0.101115508, 0.345417663, 0.202331199, 0.251746491], [-0.30818647, 0.084023406, 0.341089481, 0.200767673, 0.250394494], [-0.306686047, 0.07470381, 0.3420609, 0.202142239, 0.245542794], [-0.313876883, 0.071033507, 0.355192942, 0.204450416, 0.244608387], [-0.314688542, 0.061555107, 0.356189983, 0.205736257, 0.242747231], [-0.315962674, 0.049823613, 0.354471767, 0.201790944, 0.236930639], [-0.313927933, 0.042350961, 0.353760318, 0.199163146, 0.233971037], [-0.311733532, 0.038001598, 0.35386736, 0.200105369, 0.233066656], [-0.320665749, 0.017504398, 0.355023584, 0.200108312, 0.232755627], [-0.321798843, 0.003176592, 0.352301716, 0.204816304, 0.233308677], [-0.322225722, -0.004659656, 0.353407995, 0.207387911, 0.229388839], [-0.322062703, -0.015891113, 0.347245916, 0.209002388, 0.228579043], [-0.323754946, -0.022829634, 0.345746481, 0.208617193, 0.229073638], [-0.326085151, -0.025158168, 0.346081206, 0.206000081, 0.228036961], [-0.326322457, -0.028077633, 0.344256668, 0.202713665, 0.222630727], [-0.328791942, -0.033029489, 0.343646922, 0.200033424, 0.219921594], [-0.330966628, -0.03828204, 0.343568322, 0.201457449, 0.219779751], [-0.331991914, -0.041837648, 0.344006313, 0.203514758, 0.220175879]]}
GRAVITY = 9.80665
PERIODS = array([-1.   ,  0.   ,  0.01 ,  0.02 ,  0.03 ,  0.04 ,  0.05 ,  0.075,         0.1  ,  0.11 ,  0.12 ,  0.13 ,  0.14 ,  0.15 ,  0.16 ,  0.17 ,         0.18 ,  0.19 ,  0.2  ,  0.22 ,  0.24 ,  0.26 ,  0.28 ,  0.3  ,         0.32 ,  0.34 ,  0.36 ,  0.38 ,  0.4  ,  0.42 ,  0.44 ,  0.46 ,         0.48 ,  0.5  ,  0.55 ,  0.6  ,  0.65 ,  0.7  ,  0.75 ,  0.8  ,         0.85 ,  0.9  ,  0.95 ,  1.   ,  1.1  ,  1.2  ,  1.3  ,  1.4  ,         1.5  ,  1.6  ,  1.7  ,  1.8  ,  1.9  ,  2.   ,  2.2  ,  2.4  ,         2.6  ,  2.8  ,  3.   ,  3.2  ,  3.4  ,  3.6  ,  3.8  ,  4.   ])

Indices of the periods

INDICES_PSA = array([ 2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17, 18,        19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,        36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52,        53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63])

Indices for the spectral accelerations

INDEX_PGA = 1

Index of the peak ground acceleration

INDEX_PGV = 0

Index of the peak ground velocity

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

Model parameters

__init__(scenario)[source]

Initialize the model.

Parameters:

scenario (Scenario)

INDEX_PGD = None

Index of the peak ground displacement

LIMITS = {}

Limits of model applicability

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).