def focmec(self):
if self.ui.source_tab.currentIndex() == 0:
return [
float(self.ui.m_rr.value()),
float(self.ui.m_tt.value()),
float(self.ui.m_pp.value()),
float(self.ui.m_rt.value()),
float(self.ui.m_rp.value()),
float(self.ui.m_tp.value()),
]
elif self.ui.source_tab.currentIndex() == 1:
source = Source.from_strike_dip_rake(
latitude=float(self.ui.source_latitude.value()),
longitude=float(self.ui.source_longitude.value()),
depth_in_m=float(self.source_depth) * 1000.0,
strike=float(self.ui.strike_slider.value()),
dip=float(self.ui.dip_slider.value()),
rake=float(self.ui.rake_slider.value()),
M0=1e16,
)
return [
source.m_rr,
source.m_tt,
source.m_pp,
source.m_rt,
source.m_rp,
source.m_tp,
]
def test_print_regressions():
"""
Guard against a regression for printing a source object.
"""
src = Source.from_strike_dip_rake(
latitude=27.77,
longitude=85.37,
depth_in_m=12000.0,
M0=1e21,
strike=32.0,
dip=62.0,
rake=90.0,
)
assert (str(src)) == ("Instaseis Source:\n"
"\tOrigin Time : 1970-01-01T00:00:00.000000Z\n"
"\tLongitude : 85.4 deg\n"
"\tLatitude : 27.8 deg\n"
"\tDepth : 1.2e+01 km km\n"
"\tMoment Magnitude : 7.93\n"
"\tScalar Moment : 1.00e+21 Nm\n"
"\tMrr : 8.29e+20 Nm\n"
"\tMtt : -2.33e+20 Nm\n"
"\tMpp : -5.96e+20 Nm\n"
"\tMrt : 2.96e+20 Nm\n"
"\tMrp : 4.74e+20 Nm\n"
"\tMtp : -3.73e+20 Nm\n")
def test_m0():
"""
Tests computation of scalar Moment.
"""
strike = 10.
dip = 20.
rake = 30.
m0 = 1e16
source = Source.from_strike_dip_rake(0., 0., 0., strike, dip, rake, m0)
assert source.M0 == m0
def test_M0():
"""
Tests computation of scalar Moment.
"""
strike = 10.
dip = 20.
rake = 30.
M0 = 1e16
source = Source.from_strike_dip_rake(0., 0., 0., strike, dip, rake, M0)
assert source.M0 == M0
def test_m0():
"""
Tests computation of scalar Moment.
"""
strike = 10.0
dip = 20.0
rake = 30.0
m0 = 1e16
source = Source.from_strike_dip_rake(0.0, 0.0, 0.0, strike, dip, rake, m0)
assert source.M0 == m0
def focmec(self):
if self.ui.source_tab.currentIndex() == 0:
return [float(self.ui.m_rr.value()), float(self.ui.m_tt.value()),
float(self.ui.m_pp.value()), float(self.ui.m_rt.value()),
float(self.ui.m_rp.value()), float(self.ui.m_tp.value())]
elif self.ui.source_tab.currentIndex() == 1:
source = Source.from_strike_dip_rake(
latitude=float(self.ui.source_latitude.value()),
longitude=float(self.ui.source_longitude.value()),
depth_in_m=float(self.source_depth) * 1000.0,
strike=float(self.ui.strike_slider.value()),
dip=float(self.ui.dip_slider.value()),
rake=float(self.ui.rake_slider.value()),
M0=1E16)
return [source.m_rr, source.m_tt,
source.m_pp, source.m_rt,
source.m_rp, source.m_tp]
def test_print_regressions():
"""
Guard against a regression for printing a source object.
"""
src = Source.from_strike_dip_rake(
latitude=27.77, longitude=85.37, depth_in_m=12000.0,
M0=1e+21, strike=32., dip=62., rake=90.)
assert(str(src)) == (
"Instaseis Source:\n"
"\tOrigin Time : 1970-01-01T00:00:00.000000Z\n"
"\tLongitude : 85.4 deg\n"
"\tLatitude : 27.8 deg\n"
"\tDepth : 1.2e+01 km km\n"
"\tMoment Magnitude : 7.93\n"
"\tScalar Moment : 1.00e+21 Nm\n"
"\tMrr : 8.29e+20 Nm\n"
"\tMtt : -2.33e+20 Nm\n"
"\tMpp : -5.96e+20 Nm\n"
"\tMrt : 2.96e+20 Nm\n"
"\tMrp : 4.74e+20 Nm\n"
"\tMtp : -3.73e+20 Nm\n")
