diff --git a/geomagio/TimeseriesUtility.py b/geomagio/TimeseriesUtility.py
index da7086d967b2e2c20c817b1381d204a1039dba80..aea308c52062148bfd7625e21fc8f78c3fb257a5 100644
--- a/geomagio/TimeseriesUtility.py
+++ b/geomagio/TimeseriesUtility.py
@@ -195,17 +195,36 @@ def merge_streams(*streams):
         stream with contiguous traces merged, and gaps filled with numpy.nan
     """
     merged = obspy.core.Stream()
+    # masked_trace = None
+
     # add unmasked, split traces to be merged
     for stream in streams:
         merged += mask_stream(stream)
+
+    # if trace is completely masked separate out to be added back in later
+    # for trace in merged:
+    #     if trace.data.mask.all():
+    #         if not masked_trace:
+    #             masked_trace = trace
+    #         else:
+    #             masked_trace += trace
+
     # split traces that contain gaps
     merged = merged.split()
+
     # merge data
     merged.merge(
             # 1 = do not interpolate
             interpolation_samples=1,
             # 1 = when there is overlap, use data from trace with last endtime
             method=1)
+
+    # trim masked trace to the same size as other traces and add back to merged stream
+    # if masked_trace:
+    #     masked_trace.trim(merged[0].stats.starttime,merged[0].stats.endtime)
+    #     merged += masked_trace
+
+
     # convert back to NaN filled array
     merged = unmask_stream(merged)
     return merged