def test_extent_config():
eventfile = os.path.join(datadir, 'event_oklahoma_large.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
config = {'extent': {'coefficients': {'coeffs': [27.24, 250.4, 579.1]}}}
cmp_extent = (-100.61666666666666, -93.95, 32.916666666666664, 38.35)
extent = get_extent(rupture, config)
np.testing.assert_almost_equal(cmp_extent, extent)
config = {
'extent': {
'magnitude_spans': {
'span1': [0, 6, 4, 3],
'span2': [6, 10, 6, 4]
}
}
}
extent = get_extent(rupture, config)
cmp_extent = (-99.416, -95.416, 32.679, 38.679)
np.testing.assert_almost_equal(cmp_extent, extent)
config = {'extent': {'bounds': {'extent': [-100, 32, -95, 37]}}}
extent = get_extent(rupture, config)
cmp_extent = [-100, -95, 32, 37]
np.testing.assert_almost_equal(extent, cmp_extent)
def test_extent_config():
eventfile = os.path.join(datadir, 'event_oklahoma_large.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
config = {
'extent': {
'magnitude_spans': {
'span1': [0, 6, 4, 3],
'span2': [6, 10, 6, 4]
}
}
}
extent = get_extent(rupture, config)
cmp_extent = (-99.4166667, -95.4083333, 32.675, 38.6833333)
np.testing.assert_almost_equal(cmp_extent, extent)
# Test for relative_offset
config = {
'extent': {
'magnitude_spans': {
'span1': [0, 6, 4, 3],
'span2': [6, 10, 6, 4]
},
'relative_offset': [0.25, 0.5],
}
}
extent = get_extent(rupture, config)
cmp_extent = (-98.4166667, -94.4083333, 35.675, 41.6833333)
np.testing.assert_almost_equal(cmp_extent, extent)
config = {'extent': {'bounds': {'extent': [-100, 32, -95, 37]}}}
extent = get_extent(rupture, config)
cmp_extent = [-100, -95, 32, 37]
np.testing.assert_almost_equal(extent, cmp_extent)
def test_get_extent_bad_usage():
#
# Test bad usage
#
with pytest.raises(TypeError):
W, E, S, N = get_extent()
with pytest.raises(TypeError):
W, E, S, N = get_extent(34)
def test_get_extent_bad_usage():
#
# Test bad usage
#
with pytest.raises(TypeError):
W, E, S, N = get_extent()
with pytest.raises(TypeError):
W, E, S, N = get_extent(34)
def test_get_extent_aspect():
#
# Test ruptures with weird aspect ratios
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
#
# Long horizontal rupture
#
rrep = io.StringIO(
"""
100.0 30.0 0.0
110.0 30.0 0.0
110.0 30.0 5.0
100.0 30.0 5.0
100.0 30.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
assert W == 97.18333333333334
assert E == 113.51666666666667
assert S == 25.616666666666667
assert N == 34.06666666666666
#
# Long vertical rupture
#
rrep = io.StringIO(
"""
100.0 24.0 0.0
100.0 36.0 0.0
100.0 36.0 5.0
100.0 24.0 5.0
100.0 24.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
assert W == 92.53333333333333
assert E == 108.89999999999999
assert S == 21.03333333333333
assert N == 38.46666666666667
def test_get_extent_aspect():
#
# Test ruptures with weird aspect ratios
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
#
# Long horizontal rupture
#
rrep = io.StringIO(
"""
30.0 100.0 0.0
30.0 110.0 0.0
30.0 110.0 5.0
30.0 100.0 5.0
30.0 100.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
assert W == 97.18333333333334
assert E == 113.51666666666667
assert S == 25.616666666666667
assert N == 34.06666666666666
#
# Long vertical rupture
#
rrep = io.StringIO(
"""
24.0 100.0 0.0
36.0 100.0 0.0
36.0 100.0 5.0
24.0 100.0 5.0
24.0 100.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
assert W == 92.53333333333333
assert E == 108.89999999999999
assert S == 21.03333333333333
assert N == 38.46666666666667
def test_get_extent_aspect():
#
# Test ruptures with weird aspect ratios
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
#
# Long horizontal rupture
#
rrep = io.StringIO(
"""
100.0 30.0 0.0
110.0 30.0 0.0
110.0 30.0 5.0
100.0 30.0 5.0
100.0 30.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, 92.65)
np.testing.assert_allclose(E, 119.98333333333333)
np.testing.assert_allclose(S, 19.833333333333332)
np.testing.assert_allclose(N, 38.96666666666667)
#
# Long vertical rupture
#
rrep = io.StringIO(
"""
100.0 24.0 0.0
100.0 36.0 0.0
100.0 36.0 5.0
100.0 24.0 5.0
100.0 24.0 0.0
"""
)
rupture = get_rupture(origin, rrep)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, 88.48333333333333)
np.testing.assert_allclose(E, 115.08333333333333)
np.testing.assert_allclose(S, 16.1)
np.testing.assert_allclose(N, 42.583333333333336)
def test_get_extent_small_point():
#
# Do a point rupture
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_wenchuan_small.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == 102.33333333333333
assert E == 104.41666666666667
assert S == 30.083333333333332
assert N == 31.883333333333333
def test_get_extent_small_complex():
#
# Do a complex rupture
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
faultfile = os.path.join(datadir, 'Hartzell11_fault.txt')
rupture = get_rupture(origin, faultfile)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, 95.28333333333333)
np.testing.assert_allclose(E, 115.1)
np.testing.assert_allclose(S, 23.083333333333332)
np.testing.assert_allclose(N, 39.75)
def test_get_extent_small_point():
#
# Do a point rupture
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_wenchuan_small.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, 102.33333333333333)
np.testing.assert_allclose(E, 104.41666666666667)
np.testing.assert_allclose(S, 30.083333333333332)
np.testing.assert_allclose(N, 31.883333333333333)
def test_get_extent_stable_large():
#
# Do an event in a stable region
# Large magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma_large.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, -107.45)
np.testing.assert_allclose(E, -84.8)
np.testing.assert_allclose(S, 26.166666666666668)
np.testing.assert_allclose(N, 44.15)
def test_get_extent_stable_small():
#
# Do an event in a stable region
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
np.testing.assert_allclose(W, -98.5)
np.testing.assert_allclose(E, -96.28333333333333)
np.testing.assert_allclose(S, 34.766666666666666)
np.testing.assert_allclose(N, 36.583333333333336)
def test_get_extent_stable_small():
#
# Do an event in a stable region
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == -98.5
assert E == -96.28333333333333
assert S == 34.766666666666666
assert N == 36.583333333333336
def test_get_extent_stable_large():
#
# Do an event in a stable region
# Large magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma_large.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == -106.14999999999999
assert E == -86.76666666666667
assert S == 27.55
assert N == 43.03333333333333
def test_get_extent_small_complex():
#
# Do a complex rupture
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
faultfile = os.path.join(datadir, 'Hartzell11_fault.txt')
rupture = get_rupture(origin, faultfile)
W, E, S, N = get_extent(rupture)
assert W == 100.03333333333333
assert E == 108.93333333333334
assert S == 27.9
assert N == 35.55
def test_get_extent_stable_large():
#
# Do an event in a stable region
# Large magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma_large.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == -106.14999999999999
assert E == -86.76666666666667
assert S == 27.55
assert N == 43.03333333333333
def test_get_extent_small_point():
#
# Do a point rupture
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_wenchuan_small.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == 102.33333333333333
assert E == 104.41666666666667
assert S == 30.083333333333332
assert N == 31.883333333333333
def test_get_extent_stable_small():
#
# Do an event in a stable region
# Small magnitude
#
eventfile = os.path.join(datadir, 'event_oklahoma.xml')
origin = Origin.fromFile(eventfile)
rupture = get_rupture(origin)
W, E, S, N = get_extent(rupture)
assert W == -98.5
assert E == -96.28333333333333
assert S == 34.766666666666666
assert N == 36.583333333333336
def test_get_extent_small_complex():
#
# Do a complex rupture
#
eventfile = os.path.join(datadir, 'event_wenchuan.xml')
origin = Origin.fromFile(eventfile)
faultfile = os.path.join(datadir, 'Hartzell11_fault.txt')
rupture = get_rupture(origin, faultfile)
W, E, S, N = get_extent(rupture)
assert W == 100.03333333333333
assert E == 108.93333333333334
assert S == 27.9
assert N == 35.55
