Switch dbdt suffix from _DDT to _DT, round at 6 decimals

geomagio/algorithm/DbDtAlgorithm.py

@@ -34,11 +34,11 @@ class DbDtAlgorithm(Algorithm):
         """
         out = Stream()
         for trace in stream:
-            dbdt = np.diff(trace.data)
+            dbdt = np.around(np.diff(trace.data), decimals=6)
             stats = Stats(trace.stats)
-            stats.channel = "{}_DDT".format(stats.channel)
+            stats.channel = "{}_DT".format(stats.channel)
             trace_out = create_empty_trace(
-                starttime=stats.starttime + self.period,
+                starttime=stats.starttime + stats.delta,
                 endtime=stats.endtime,
                 observatory=stats.station,
                 type=stats.location,

geomagio/edge/MiniSeedFactory.py

@@ -349,7 +349,7 @@ class MiniSeedFactory(TimeseriesFactory):
                 return channel
 
         if suffix is not None:
-            if suffix == "Dist" or suffix == "SQ" or suffix == "SV" or suffix == "DDT":
+            if suffix == "Dist" or suffix == "SQ" or suffix == "SV" or suffix == "DT":
                 # these suffixes modify location code, but use element channel
                 channel = element
             else:
@@ -421,7 +421,7 @@ class MiniSeedFactory(TimeseriesFactory):
                 location_suffix = "Q"
             elif suffix == "SV":
                 location_suffix = "V"
-            elif suffix == "DDT":
+            elif suffix == "DT":
                 location_suffix = "R"
             elif suffix not in ("Bin", "Volt"):
                 raise TimeseriesFactoryException(

test/algorithm_test/DbDtAlgorithm_test.py

@@ -10,7 +10,7 @@ def test_process():
     original script
     """
     # initialize DbDt object
-    dbdt = DbDtAlgorithm(inchannels=["H"], outchannels=["H_DDT"], period=60)
+    dbdt = DbDtAlgorithm(inchannels=["H"], outchannels=["H_DT"], period=60)
 
     # load boulder May 20 files from /etc/ directory
     hez_iaga2002_file = open("etc/dbdt/BOU202005vmin.min")
@@ -25,7 +25,7 @@ def test_process():
     result = dbdt.process(hez)
 
     # unpack channels from result
-    rh = result.select(channel="H_DDT")[0]
+    rh = result.select(channel="H_DT")[0]
     # unpack channels from BOU202005dbdt.min
     h = hez_dbdt.select(channel="H")[0]