def create_oqnrml_source(self, use_defaults=False):
'''
Turns source into instance of the :class:
openquake.nrmllib.model.SimpleFaultSource
:param bool use_defaults:
If set to True, will use default values for rupture aspect ratio
and magnitude scaling relation. If False, will raise value error
if this information is missing
'''
if not isinstance(self.rake, float):
raise ValueError('Cannot render fault source - rake is missing!')
simple_geometry = models.SimpleFaultGeometry(
wkt=conv.simple_trace_to_wkt_linestring(self.fault_trace),
dip=self.dip,
upper_seismo_depth=self.upper_depth,
lower_seismo_depth=self.lower_depth)
return models.SimpleFaultSource(
self.id,
self.name,
self.trt,
simple_geometry,
conv.render_mag_scale_rel(self.mag_scale_rel, use_defaults),
conv.render_aspect_ratio(self.rupt_aspect_ratio, use_defaults),
conv.render_mfd(self.mfd),
self.rake)
def test_create_oqnmrl_simple_fault_source(self):
'''
Tests the conversion of a point source to an instance of the :class:
oqnrmllib.models.SimpleFaultSource
'''
trace = line.Line([point.Point(10., 10.), point.Point(11., 10.)])
#sflt_geom = SimpleFaultSurface.from_fault_data(trace, 0., 20., 90., 1.)
# Define a complete source
self.fault_source = mtkSimpleFaultSource('001',
'A Fault Source',
trt='Active Shallow Crust',
geometry=None,
dip=90.,
upper_depth=0.,
lower_depth=20.,
mag_scale_rel=None,
rupt_aspect_ratio=1.0,
mfd=models.TGRMFD(
a_val=3.,
b_val=1.0,
min_mag=5.0,
max_mag=8.0),
rake=0.)
self.fault_source.create_geometry(trace, 90., 0., 20., 1.0)
expected_geometry = models.SimpleFaultGeometry(
wkt='LINESTRING (10.0 10.0, 11.0 10.0)',
dip=90.,
upper_seismo_depth=0.,
lower_seismo_depth=20.)
expected_source = models.SimpleFaultSource('001',
'A Fault Source',
trt='Active Shallow Crust',
geometry=expected_geometry,
mag_scale_rel='WC1994',
rupt_aspect_ratio=1.0,
mfd=models.TGRMFD(
a_val=3.,
b_val=1.0,
min_mag=5.0,
max_mag=8.0),
rake=0.)
test_source = self.fault_source.create_oqnrml_source(use_defaults=True)
self.assertTrue(isinstance(test_source, models.SimpleFaultSource))
self.assertEqual(test_source.id, expected_source.id)
self.assertEqual(test_source.name, expected_source.name)
self.assertDictEqual(test_source.geometry.__dict__,
expected_source.geometry.__dict__)
self.assertAlmostEqual(test_source.mfd.b_val,
expected_source.mfd.b_val)
def _parse_simple(cls, src_elem):
"""
:param src_elem:
:class:`lxml.etree._Element` instance representing a source.
:returns:
Fully populated
:class:`openquake.nrmllib.models.SimpleFaultSource` object.
"""
simple = models.SimpleFaultSource()
cls._set_common_attrs(simple, src_elem)
simple_geom = cls._parse_simple_geometry(src_elem)
simple.geometry = simple_geom
simple.mfd = cls._parse_mfd(src_elem)
simple.rake = float(_xpath(src_elem, './/nrml:rake')[0].text)
return simple
def _expected_source_model(cls):
# Area:
area_geom = models.AreaGeometry(
wkt=('POLYGON((-122.5 37.5, -121.5 37.5, -121.5 38.5, -122.5 38.5,'
' -122.5 37.5))'),
upper_seismo_depth=0.0,
lower_seismo_depth=10.0,
)
area_mfd = models.IncrementalMFD(
min_mag=6.55,
bin_width=0.1,
occur_rates=[
0.0010614989, 8.8291627E-4, 7.3437777E-4, 6.108288E-4,
5.080653E-4
],
)
area_npd = [
models.NodalPlane(probability=decimal.Decimal("0.3"),
strike=0.0,
dip=90.0,
rake=0.0),
models.NodalPlane(probability=decimal.Decimal("0.7"),
strike=90.0,
dip=45.0,
rake=90.0),
]
area_hdd = [
models.HypocentralDepth(probability=decimal.Decimal("0.5"),
depth=4.0),
models.HypocentralDepth(probability=decimal.Decimal("0.5"),
depth=8.0),
]
area_src = models.AreaSource(
id='1',
name='Quito',
trt='Active Shallow Crust',
geometry=area_geom,
mag_scale_rel='PeerMSR',
rupt_aspect_ratio=1.5,
mfd=area_mfd,
nodal_plane_dist=area_npd,
hypo_depth_dist=area_hdd,
)
# Point:
point_geom = models.PointGeometry(
wkt='POINT(-122.0 38.0)',
upper_seismo_depth=0.0,
lower_seismo_depth=10.0,
)
point_mfd = models.TGRMFD(
a_val=-3.5,
b_val=1.0,
min_mag=5.0,
max_mag=6.5,
)
point_npd = [
models.NodalPlane(probability=decimal.Decimal("0.3"),
strike=0.0,
dip=90.0,
rake=0.0),
models.NodalPlane(probability=decimal.Decimal("0.7"),
strike=90.0,
dip=45.0,
rake=90.0),
]
point_hdd = [
models.HypocentralDepth(probability=decimal.Decimal("0.5"),
depth=4.0),
models.HypocentralDepth(probability=decimal.Decimal("0.5"),
depth=8.0),
]
point_src = models.PointSource(
id='2',
name='point',
trt='Stable Continental Crust',
geometry=point_geom,
mag_scale_rel='WC1994',
rupt_aspect_ratio=0.5,
mfd=point_mfd,
nodal_plane_dist=point_npd,
hypo_depth_dist=point_hdd,
)
# Simple:
simple_geom = models.SimpleFaultGeometry(
wkt='LINESTRING(-121.82290 37.73010, -122.03880 37.87710)',
dip=45.0,
upper_seismo_depth=10.0,
lower_seismo_depth=20.0,
)
simple_mfd = models.IncrementalMFD(
min_mag=5.0,
bin_width=0.1,
occur_rates=[
0.0010614989, 8.8291627E-4, 7.3437777E-4, 6.108288E-4,
5.080653E-4
],
)
simple_src = models.SimpleFaultSource(
id='3',
name='Mount Diablo Thrust',
trt='Active Shallow Crust',
geometry=simple_geom,
mag_scale_rel='WC1994',
rupt_aspect_ratio=1.5,
mfd=simple_mfd,
rake=30.0,
)
# Complex:
complex_geom = models.ComplexFaultGeometry(
top_edge_wkt=('LINESTRING(-124.704 40.363 0.5493260E+01, '
'-124.977 41.214 0.4988560E+01, '
'-125.140 42.096 0.4897340E+01)'),
bottom_edge_wkt=('LINESTRING(-123.829 40.347 0.2038490E+02, '
'-124.137 41.218 0.1741390E+02, '
'-124.252 42.115 0.1752740E+02)'),
int_edges=[
('LINESTRING(-124.704 40.363 0.5593260E+01, '
'-124.977 41.214 0.5088560E+01, '
'-125.140 42.096 0.4997340E+01)'),
('LINESTRING(-124.704 40.363 0.5693260E+01, '
'-124.977 41.214 0.5188560E+01, '
'-125.140 42.096 0.5097340E+01)'),
])
complex_mfd = models.TGRMFD(a_val=-3.5,
b_val=1.0,
min_mag=5.0,
max_mag=6.5)
complex_src = models.ComplexFaultSource(
id='4',
name='Cascadia Megathrust',
trt='Subduction Interface',
geometry=complex_geom,
mag_scale_rel='WC1994',
rupt_aspect_ratio=2.0,
mfd=complex_mfd,
rake=30.0,
)
# 3 Characteristic Sources:
char_src_simple = models.CharacteristicSource(
id='5',
name='characteristic source, simple fault',
trt='Volcanic',
mfd=models.TGRMFD(a_val=-3.5, b_val=1.0, min_mag=5.0, max_mag=6.5),
rake=30.0,
surface=simple_geom)
char_src_complex = models.CharacteristicSource(
id='6',
name='characteristic source, complex fault',
trt='Volcanic',
mfd=models.IncrementalMFD(
min_mag=5.0,
bin_width=0.1,
occur_rates=[
0.0010614989, 8.8291627E-4, 7.3437777E-4, 6.108288E-4,
5.080653E-4
],
),
rake=60.0,
surface=complex_geom)
char_src_multi = models.CharacteristicSource(
id='7',
name='characteristic source, multi surface',
trt='Volcanic',
mfd=models.TGRMFD(a_val=-3.6, b_val=1.0, min_mag=5.2, max_mag=6.4),
rake=90.0)
psurface_1 = models.PlanarSurface(
strike=0.0,
dip=90.0,
top_left=models.Point(longitude=-1.0, latitude=1.0, depth=21.0),
top_right=models.Point(longitude=1.0, latitude=1.0, depth=21.0),
bottom_left=models.Point(longitude=-1.0, latitude=-1.0,
depth=59.0),
bottom_right=models.Point(longitude=1.0, latitude=-1.0,
depth=59.0),
)
psurface_2 = models.PlanarSurface(
strike=20.0,
dip=45.0,
top_left=models.Point(longitude=1.0, latitude=1.0, depth=20.0),
top_right=models.Point(longitude=3.0, latitude=1.0, depth=20.0),
bottom_left=models.Point(longitude=1.0, latitude=-1.0, depth=80.0),
bottom_right=models.Point(longitude=3.0, latitude=-1.0,
depth=80.0),
)
char_src_multi.surface = [psurface_1, psurface_2]
source_model = models.SourceModel()
source_model.name = 'Some Source Model'
# Generator:
source_model.sources = (x for x in [
area_src, point_src, simple_src, complex_src, char_src_simple,
char_src_complex, char_src_multi
])
return source_model
