def test_read_single_ndk(self):
"""
Test reading a single event from an NDK file and comparing it to a
QuakeML file that has been manually checked to contain all the
information in the NDK file.
"""
filename = os.path.join(self.datapath, "C200604092050A.ndk")
cat = _read_ndk(filename)
reference = os.path.join(self.datapath, "C200604092050A.xml")
ref_cat = read_events(reference)
self.assertEqual(cat, ref_cat)
def test_reading_from_bytestring(self):
"""
Tests reading from a byte string.
"""
filename = os.path.join(self.datapath, "C200604092050A.ndk")
reference = os.path.join(self.datapath, "C200604092050A.xml")
ref_cat = read_events(reference)
with io.open(filename, "rb") as fh:
data = fh.read()
self.assertTrue(_is_ndk(data))
cat = _read_ndk(data)
self.assertEqual(cat, ref_cat)
def test_read_multiple_events(self):
"""
Tests the reading of multiple events in one file. The file has been
edited to test a variety of settings.
"""
filename = os.path.join(self.datapath, "multiple_events.ndk")
cat = _read_ndk(filename)
self.assertEqual(len(cat), 6)
# Test the type of moment tensor inverted for.
self.assertEqual([i.focal_mechanisms[0].moment_tensor.inversion_type
for i in cat],
["general", "zero trace", "double couple"] * 2)
# Test the type and duration of the moment rate function.
self.assertEqual(
[i.focal_mechanisms[0].moment_tensor.source_time_function.type
for i in cat],
["triangle", "box car"] * 3)
self.assertEqual(
[i.focal_mechanisms[0].moment_tensor.source_time_function.duration
for i in cat],
[2.6, 7.4, 9.0, 1.8, 2.0, 1.6])
# Test the type of depth setting.
self.assertEqual([i.preferred_origin().depth_type for i in cat],
["from moment tensor inversion", "from location",
"from modeling of broad-band P waveforms"] * 2)
# Check the solution type.
for event in cat[:3]:
self.assertIn("Standard",
event.focal_mechanisms[0].comments[0].text)
for event in cat[3:]:
self.assertIn("Quick",
event.focal_mechanisms[0].comments[0].text)
