def test_properties_force_source():
"""
Tests some properties of the force source.
"""
src = ForceSource(latitude=0.0, longitude=0.0, f_r=1.0, f_t=2.0, f_p=3.0)
np.testing.assert_allclose(src.force_tpr, [2.0, 3.0, 1.0])
np.testing.assert_allclose(src.force_rtp, [1.0, 2.0, 3.0])
def test_str_method_of_force_source():
src = ForceSource(latitude=0.0, longitude=0.0, f_r=1.0, f_t=2.0, f_p=3.0)
assert str(src) == ("Instaseis Force Source:\n"
"\tOrigin Time : 1970-01-01T00:00:00.000000Z\n"
"\tLongitude : 0.0 deg\n"
"\tLatitude : 0.0 deg\n"
"\tFr : 1.00e+00 N\n"
"\tFt : 2.00e+00 N\n"
"\tFp : 3.00e+00 N\n")
def test_sliprate_convenience_methods_force_source():
"""
Tests some convenience methods of sliprates for force sources.
"""
src = ForceSource(latitude=0.0, longitude=90.0)
src.set_sliprate_dirac(2.0, 5)
np.testing.assert_allclose(np.array([0.5, 0, 0, 0, 0]), src.sliprate)
src = ForceSource(latitude=0.0, longitude=90.0)
src.set_sliprate_lp(2.0, 5, 0.1)
np.testing.assert_allclose(np.array(
[0.023291, 0.111382, 0.211022, 0.186723, 0.045481]), src.sliprate,
rtol=1E-3)
src = ForceSource(latitude=0.0, longitude=90.0)
src.sliprate = np.ones(5)
src.dt = 0.25
src.normalize_sliprate()
np.testing.assert_allclose(np.ones(5), src.sliprate)
def test_sliprate_convenience_methods_force_source():
"""
Tests some convenience methods of sliprates for force sources.
"""
src = ForceSource(latitude=0.0, longitude=90.0)
src.set_sliprate_dirac(2.0, 5)
np.testing.assert_allclose(np.array([0.5, 0, 0, 0, 0]), src.sliprate)
src = ForceSource(latitude=0.0, longitude=90.0)
src.set_sliprate_lp(2.0, 5, 0.1)
np.testing.assert_allclose(
np.array([0.023291, 0.111382, 0.211022, 0.186723, 0.045481]),
src.sliprate,
rtol=1e-3,
)
src = ForceSource(latitude=0.0, longitude=90.0)
src.sliprate = np.ones(5)
src.dt = 0.25
src.normalize_sliprate()
np.testing.assert_allclose(np.ones(5), src.sliprate)