From 98ceb6bec8ab0a7ea4a70b4f8db54aac0ac4231e Mon Sep 17 00:00:00 2001
From: kienle <kienle@passcal.nmt.edu>
Date: Mon, 7 Aug 2023 16:47:53 -0600
Subject: [PATCH] Remove unused code
---
.../model/reftek/from_rt2ms/core.py | 83 ++++---------------
1 file changed, 15 insertions(+), 68 deletions(-)
diff --git a/sohstationviewer/model/reftek/from_rt2ms/core.py b/sohstationviewer/model/reftek/from_rt2ms/core.py
index 221b1d05b..8cdcdf584 100644
--- a/sohstationviewer/model/reftek/from_rt2ms/core.py
+++ b/sohstationviewer/model/reftek/from_rt2ms/core.py
@@ -11,6 +11,9 @@ Maeva Pourpoint IRIS/PASSCAL
"""
import copy
+import os
+import sys
+
import obspy.io.reftek.core as obspy_rt130_core
import warnings
@@ -38,74 +41,6 @@ class DiscontinuousTrace(Trace):
self.times = times
-
-EH_PAYLOAD = {
- "station_name_extension": (35, 1, _decode_ascii),
- "station_name": (36, 4, _decode_ascii),
- "sampling_rate": (64, 4, float),
- "trigger_time": (72, 16, _parse_long_time),
- "first_sample_time": (
- eh_et_payload_last_field_start, eh_et_payload_last_field_size, _parse_long_time),
-}
-
-
-# class EHPacket(packet.EHPacket):
-# def __init__(self, data):
-# self._data = data
-# payload = self._data["payload"].tobytes()
-# for name, (start, length, converter) in EH_PAYLOAD.items():
-# data = payload[start:start + length]
-# if converter is not None:
-# data = converter(data)
-# setattr(self, name, data)
-#
-# def eh_et_info(self, nbr_DT_samples):
-# """
-# Compile EH and ET info to write to log file.
-# Returns list of strings.
-# Formatting of strings is based on earlier version of rt2ms.
-# """
-# info = []
-# # packet_tagline1 = ("\n\n{:s} exp {:02d} bytes {:04d} {:s} ID: {:s} "
-# # "seq {:04d}".format(self.type.decode(),
-# # self.experiment_number,
-# # self.byte_count,
-# # Packet.time_tag(self.time),
-# # self.unit_id.decode(),
-# # self.packet_sequence))
-# # info.append(packet_tagline1)
-# # if self.type.decode('ASCII') == 'EH':
-# # nbr_DT_samples = 0
-# # info.append("\nEvent Header")
-# # else:
-# # info.append("\nEvent Trailer")
-# # info.append("\n event = " + str(self.event_number))
-# # info.append("\n stream = " + str(self.data_stream_number + 1))
-# # info.append("\n format = " + self.data_format.decode('ASCII'))
-# # info.append("\n stream name = " + self.stream_name)
-# # info.append("\n sample rate = " + str(self.sampling_rate))
-# # info.append("\n trigger type = " + self.trigger_type)
-# trigger_time = Packet.time_tag(UTCDateTime(ns=self.trigger_time))
-# # info.append("\n trigger time = " + trigger_time)
-# first_sample_time = Packet.time_tag(UTCDateTime(ns=self.first_sample_time)) # noqa: E501
-# # info.append("\n first sample = " + first_sample_time)
-# # if self.last_sample_time:
-# # info.append("\n last sample = " + Packet.time_tag(UTCDateTime(ns=self.last_sample_time))) # noqa: E501
-# # info.append("\n bit weights = " + " ".join([val for val in self.channel_adjusted_nominal_bit_weights if val])) # noqa: E501
-# # info.append("\n true weights = " + " ".join([val for val in self.channel_true_bit_weights if val])) # noqa: E501
-# packet_tagline2 = ("\nDAS: {:s} EV: {:04d} DS: {:d} FST = {:s} TT = "
-# "{:s} NS: {:d} SPS: {:.1f} ETO: 0"
-# .format(self.unit_id.decode(),
-# self.event_number,
-# self.data_stream_number + 1,
-# first_sample_time,
-# trigger_time,
-# nbr_DT_samples,
-# self.sampling_rate))
-# info.append(packet_tagline2)
-# return info
-
-
class Reftek130Exception(ObsPyException):
pass
@@ -267,3 +202,15 @@ class Reftek130(obspy_rt130_core.Reftek130):
# check if endtime of trace is consistent
st += tr
return st
+
+
+data_path = '/Users/kle/Documents/SOHView data/RT130/2017149.92EB'
+for root, subdirs, files in os.walk(data_path):
+ for file_name in files:
+ if file_name.startswith('.'):
+ continue
+ file_path = os.path.join(root, file_name)
+ print(file_path)
+ a = (Reftek130.from_file(file_path)._data)
+ print(EHPacket(a[0]))
+ sys.exit()
--
GitLab