import_fdsnws_eq: Return depths in km, not metres

shakyground2/import_fdsnws_eq.py

@@ -105,6 +105,9 @@ def fetch_origin(root, ns, preferredoriginID):
                         d[tag] = float(v)
                     except ValueError:
                         pass
+
+        # QuakeML depths are in metres.
+        d["depth"] = d["depth"] / 1000.0
     return d
 
 

tests/test_quakexml.py

@@ -24,10 +24,15 @@ class QuakeMLReadTestCase(unittest.TestCase):
     outfile = "data/testoutput.yaml"
 
     def test_read_good_file(self):
+        """
+        This GEOFON event has a focal mechanism, and a depth of 34 km.
+        Check that a value in km, not metres, is produced.
+        """
         infile = "data/gfz2021ekhv.xml"
         result = fetch_quakeml(infile)
         self.assertTrue(isinstance(result, dict))
         self.assertEqual(result["id"], "gfz2021ekhv")
+        self.assertLess(result["origin"]["depth"], 100.0)
 
     def test_read_good_eventid(self):
         """