def test_international_date_line_2(self):
# from a bug affecting a calculation in New Zealand
fname = gettemp(characteric_source)
[[src]] = nrml.to_python(fname)
os.remove(fname)
maxdist = IntegrationDistance({'default': 200})
sitecol = SiteCollection([
Site(location=Point(176.919, -39.489),
vs30=760, vs30measured=True, z1pt0=100, z2pt5=5)])
srcfilter = SourceFilter(sitecol, maxdist)
sites = srcfilter.get_close_sites(src)
self.assertIsNotNone(sites)
def test_international_date_line_2(self):
# from a bug affecting a calculation in New Zealand
fname = gettemp(characteric_source)
[[src]] = nrml.to_python(fname)
os.remove(fname)
maxdist = IntegrationDistance({'default': 200})
sitecol = SiteCollection([
Site(location=Point(176.919, -39.489),
vs30=760, vs30measured=True, z1pt0=100, z2pt5=5)])
srcfilter = SourceFilter(sitecol, maxdist)
sites = srcfilter.get_close_sites(src)
self.assertIsNotNone(sites)
def test_point_sources(self):
sources = [
openquake.hazardlib.source.PointSource(
source_id='point1',
name='point1',
tectonic_region_type=const.TRT.ACTIVE_SHALLOW_CRUST,
mfd=openquake.hazardlib.mfd.EvenlyDiscretizedMFD(
min_mag=4, bin_width=1, occurrence_rates=[5]),
nodal_plane_distribution=openquake.hazardlib.pmf.PMF([
(1,
openquake.hazardlib.geo.NodalPlane(strike=0.0,
dip=90.0,
rake=0.0))
]),
hypocenter_distribution=openquake.hazardlib.pmf.PMF([(1, 10)]),
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
magnitude_scaling_relationship=openquake.hazardlib.scalerel.
PeerMSR(),
rupture_aspect_ratio=2,
temporal_occurrence_model=PoissonTOM(1.),
rupture_mesh_spacing=1.0,
location=Point(10, 10)),
openquake.hazardlib.source.PointSource(
source_id='point2',
name='point2',
tectonic_region_type=const.TRT.ACTIVE_SHALLOW_CRUST,
mfd=openquake.hazardlib.mfd.EvenlyDiscretizedMFD(
min_mag=4, bin_width=2, occurrence_rates=[5, 6, 7]),
nodal_plane_distribution=openquake.hazardlib.pmf.PMF([
(1,
openquake.hazardlib.geo.NodalPlane(strike=0,
dip=90,
rake=0.0)),
]),
hypocenter_distribution=openquake.hazardlib.pmf.PMF([(1, 10)]),
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
magnitude_scaling_relationship=openquake.hazardlib.scalerel.
PeerMSR(),
rupture_aspect_ratio=2,
temporal_occurrence_model=PoissonTOM(1.),
rupture_mesh_spacing=1.0,
location=Point(10, 11)),
]
sites = [
openquake.hazardlib.site.Site(Point(11, 10), 1, 2, 3),
openquake.hazardlib.site.Site(Point(10, 16), 2, 2, 3),
openquake.hazardlib.site.Site(Point(10, 10.6, 1), 3, 2, 3),
openquake.hazardlib.site.Site(Point(10, 10.7, -1), 4, 2, 3)
]
sitecol = openquake.hazardlib.site.SiteCollection(sites)
gsims = {const.TRT.ACTIVE_SHALLOW_CRUST: SadighEtAl1997()}
truncation_level = 1
imts = {'PGA': [0.1, 0.5, 1.3]}
s_filter = SourceFilter(sitecol, {const.TRT.ACTIVE_SHALLOW_CRUST: 30})
result = calc_hazard_curves(sources, s_filter, imts, gsims,
truncation_level)['PGA']
# there are two sources and four sites. The first source contains only
# one rupture, the second source contains three ruptures.
#
# the first source has 'maximum projection radius' of 0.707 km
# the second source has 'maximum projection radius' of 500.0 km
#
# the epicentral distances for source 1 are: [ 109.50558394,
# 667.16955987, 66.71695599, 77.83644865]
# the epicentral distances for source 2 are: [ 155.9412148 ,
# 555.97463322, 44.47797066, 33.35847799]
#
# Considering that the source site filtering distance is set to 30 km,
# for source 1, all sites have epicentral distance larger than
# 0.707 + 30 km. This means that source 1 ('point 1') is not considered
# in the calculation because too far.
# for source 2, the 1st, 3rd and 4th sites have epicentral distances
# smaller than 500.0 + 30 km. This means that source 2 ('point 2') is
# considered in the calculation for site 1, 3, and 4.
#
# JB distances for rupture 1 in source 2 are: [ 155.43860273,
# 555.26752644, 43.77086388, 32.65137121]
# JB distances for rupture 2 in source 2 are: [ 150.98882575,
# 548.90356541, 37.40690285, 26.28741018]
# JB distances for rupture 3 in source 2 are: [ 109.50545819,
# 55.97463322, 0. , 0. ]
#
# Considering that the rupture site filtering distance is set to 30 km,
# rupture 1 (magnitude 4) is not considered because too far, rupture 2
# (magnitude 6) affect only the 4th site, rupture 3 (magnitude 8)
# affect the 3rd and 4th sites.
self.assertEqual(result.shape, (4, 3)) # 4 sites, 3 levels
numpy.testing.assert_allclose(result[0], 0) # no contrib to site 1
numpy.testing.assert_allclose(result[1], 0) # no contrib to site 2
# test that depths are kept after filtering (sites 3 and 4 remain)
s_filter = SourceFilter(sitecol, {'default': 100})
numpy.testing.assert_array_equal(
s_filter.get_close_sites(sources[0]).depths, ([1, -1]))
def test_point_sources(self):
sources = [
openquake.hazardlib.source.PointSource(
source_id='point1', name='point1',
tectonic_region_type=const.TRT.ACTIVE_SHALLOW_CRUST,
mfd=openquake.hazardlib.mfd.EvenlyDiscretizedMFD(
min_mag=4, bin_width=1, occurrence_rates=[5]
),
nodal_plane_distribution=openquake.hazardlib.pmf.PMF([
(1, openquake.hazardlib.geo.NodalPlane(strike=0.0,
dip=90.0,
rake=0.0))
]),
hypocenter_distribution=openquake.hazardlib.pmf.PMF([(1, 10)]),
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
magnitude_scaling_relationship=
openquake.hazardlib.scalerel.PeerMSR(),
rupture_aspect_ratio=2,
temporal_occurrence_model=PoissonTOM(1.),
rupture_mesh_spacing=1.0,
location=Point(10, 10)
),
openquake.hazardlib.source.PointSource(
source_id='point2', name='point2',
tectonic_region_type=const.TRT.ACTIVE_SHALLOW_CRUST,
mfd=openquake.hazardlib.mfd.EvenlyDiscretizedMFD(
min_mag=4, bin_width=2, occurrence_rates=[5, 6, 7]
),
nodal_plane_distribution=openquake.hazardlib.pmf.PMF([
(1, openquake.hazardlib.geo.NodalPlane(strike=0,
dip=90,
rake=0.0)),
]),
hypocenter_distribution=openquake.hazardlib.pmf.PMF([(1, 10)]),
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
magnitude_scaling_relationship=
openquake.hazardlib.scalerel.PeerMSR(),
rupture_aspect_ratio=2,
temporal_occurrence_model=PoissonTOM(1.),
rupture_mesh_spacing=1.0,
location=Point(10, 11)
),
]
sites = [openquake.hazardlib.site.Site(Point(11, 10), 1, 2, 3),
openquake.hazardlib.site.Site(Point(10, 16), 2, 2, 3),
openquake.hazardlib.site.Site(Point(10, 10.6, 1), 3, 2, 3),
openquake.hazardlib.site.Site(Point(10, 10.7, -1), 4, 2, 3)]
sitecol = openquake.hazardlib.site.SiteCollection(sites)
gsims = {const.TRT.ACTIVE_SHALLOW_CRUST: SadighEtAl1997()}
truncation_level = 1
imts = {'PGA': [0.1, 0.5, 1.3]}
s_filter = SourceFilter(sitecol, {const.TRT.ACTIVE_SHALLOW_CRUST: 30})
result = calc_hazard_curves(
sources, s_filter, imts, gsims, truncation_level)['PGA']
# there are two sources and four sites. The first source contains only
# one rupture, the second source contains three ruptures.
#
# the first source has 'maximum projection radius' of 0.707 km
# the second source has 'maximum projection radius' of 500.0 km
#
# the epicentral distances for source 1 are: [ 109.50558394,
# 667.16955987, 66.71695599, 77.83644865]
# the epicentral distances for source 2 are: [ 155.9412148 ,
# 555.97463322, 44.47797066, 33.35847799]
#
# Considering that the source site filtering distance is set to 30 km,
# for source 1, all sites have epicentral distance larger than
# 0.707 + 30 km. This means that source 1 ('point 1') is not considered
# in the calculation because too far.
# for source 2, the 1st, 3rd and 4th sites have epicentral distances
# smaller than 500.0 + 30 km. This means that source 2 ('point 2') is
# considered in the calculation for site 1, 3, and 4.
#
# JB distances for rupture 1 in source 2 are: [ 155.43860273,
# 555.26752644, 43.77086388, 32.65137121]
# JB distances for rupture 2 in source 2 are: [ 150.98882575,
# 548.90356541, 37.40690285, 26.28741018]
# JB distances for rupture 3 in source 2 are: [ 109.50545819,
# 55.97463322, 0. , 0. ]
#
# Considering that the rupture site filtering distance is set to 30 km,
# rupture 1 (magnitude 4) is not considered because too far, rupture 2
# (magnitude 6) affect only the 4th site, rupture 3 (magnitude 8)
# affect the 3rd and 4th sites.
self.assertEqual(result.shape, (4, 3)) # 4 sites, 3 levels
numpy.testing.assert_allclose(result[0], 0) # no contrib to site 1
numpy.testing.assert_allclose(result[1], 0) # no contrib to site 2
# test that depths are kept after filtering (sites 3 and 4 remain)
s_filter = SourceFilter(sitecol, {'default': 100})
numpy.testing.assert_array_equal(
s_filter.get_close_sites(sources[0]).depths, ([1, -1]))
