gmprocess 2.4.0__tar.gz → 2.5.0__tar.gz

{gmprocess-2.4.0 → gmprocess-2.5.0}/CHANGELOG.md

@@ -1,5 +1,15 @@
 ## main
 
+## 2.5.0 / 2025-05-21
+
+- Enhancements
+  - Use location code to distinguish among instruments at a site.
+- Bugfix
+  - Bugfixes in the COSMOS fetcher query values.
+  - The `any_trace_failures` option was not working as intended. See issue #1261 and the merge request that closes it.
+- Dependencies
+  - Require usgs-strec >= 2.3.11, update tests to reflect strec fixes.
+
 ## 2.4.0 / 2025-04-04
 
 - Enhancements

{gmprocess-2.4.0 → gmprocess-2.5.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gmprocess
-Version: 2.4.0
+Version: 2.5.0
 Summary: USGS Automated Ground Motion Processing Software
 Author-email: Eric Thompson <emthompsone@usgs.gov>, Mike Hearne <mhearne@usgs.gov>, Brad Aagaard <baagaard@usgs.gov>, Bruce Worden <cbworden@contractor.usgs.gov>, John Rekoske <jrekoske@ucsd.edu>, Heather Hunsinger <hhunsinger@usgs.gov>, Gabe Ferragut <gferragut@usgs.gov>
 License: License
@@ -79,7 +79,7 @@ Requires-Dist: scipy>=1.7
 Requires-Dist: setuptools-scm>=6.3.2
 Requires-Dist: statsmodels>=0.12.2
 Requires-Dist: tqdm>4.67
-Requires-Dist: usgs-strec>=2.3.2
+Requires-Dist: usgs-strec>=2.3.11
 Requires-Dist: xlrd>=2.0
 Provides-Extra: dev
 Requires-Dist: build>=0.7.0; extra == "dev"

{gmprocess-2.4.0 → gmprocess-2.5.0}/code.json

@@ -1,4 +1,45 @@
 [
+    {
+        "name": "groundmotion-processing",
+        "description": "Read and process ground motion waveform data",
+        "version": "v2.5.0",
+        "status": "Production",
+        "organization": "U.S. Geological Survey",
+        "vcs": "git",
+        "repositoryURL": "https://code.usgs.gov/ghsc/esi/groundmotion-processing.git",
+        "homepageURL": "https://code.usgs.gov/ghsc/esi/groundmotion-processing",
+        "downloadURL": "https://code.usgs.gov/ghsc/esi/groundmotion-processing/-/archive/v2.5.0/groundmotion-processing-v2.5.0.zip",
+        "disclaimerURL": "https://code.usgs.gov/ghsc/esi/groundmotion-processing/-/raw/v2.5.0/DISCLAIMER.md",
+        "permissions": {
+            "licenses": [
+                {
+                    "name": "Public Domain, CC0-1.0",
+                    "URL": "https://code.usgs.gov/ghsc/esi/groundmotion-processing/-/raw/v2.5.0/LICENSE.md"
+                }
+            ],
+            "usageType": "openSource",
+            "exemptionText": null
+        },
+        "laborHours": 1,
+        "languages": [
+            "python"
+        ],
+        "tags": [
+            "earthquake",
+            "impact",
+            "strong motion",
+            "amplitudes",
+            "ground motion",
+            "usg-artificial-intelligence"
+        ],
+        "contact": {
+            "name": "Eric Thompson",
+            "email": "emthompson@usgs.gov"
+        },
+        "date": {
+            "metadataLastUpdated": "2025-05-21"
+        }
+    },
     {
         "name": "groundmotion-processing",
         "description": "Read and process ground motion waveform data",

{gmprocess-2.4.0 → gmprocess-2.5.0}/pyproject.toml

@@ -45,7 +45,7 @@ dependencies = [
     "setuptools-scm>=6.3.2",
     "statsmodels>=0.12.2",
     "tqdm>4.67",
-    "usgs-strec>=2.3.2",
+    "usgs-strec>=2.3.11",
     "xlrd>=2.0",
 ]
 

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/core/stationstream.py

@@ -225,7 +225,7 @@ class StationStream(Stream):
             all_codes = []
             for tr in self:
                 stats = tr.stats
-                all_codes.append(f"{stats.network}.{stats.station}.{stats.channel}")
+                all_codes.append(f"{stats.network}.{stats.station}.{stats.location}.{stats.channel}")
             if len(set(all_codes)) != len(all_codes):
                 for tr in self:
                     tr.fail("Nonunique channel code in StationStream.")
@@ -233,37 +233,43 @@ class StationStream(Stream):
     def __check_id(self):
         # Check that id is consistent, and set id if it passes the check.
         if len(self):
-            stats = self.traces[0].stats
-            if hasattr(self, "use_array") and self.use_array:
-                trace = self[0]
-                id_str = ".".join(
-                    [
-                        trace.stats.network,
-                        trace.stats.station,
-                        trace.stats.location,
-                        trace.stats.channel,
-                    ],
-                )
+            use_array = hasattr(self, "use_array") and self.use_array
+            stats = self[0].stats
+            if use_array:
+                id_str = ".".join([
+                    stats.network,
+                    stats.station,
+                    stats.location,
+                    stats.channel,
+                ])
             else:
-                id_str = f"{stats.network}.{stats.station}.{stats.channel[0:2]}"
+                id_str = ".".join([
+                    stats.network,
+                    stats.station,
+                    stats.location,
+                    stats.channel[0:2],
+                ])
 
             # Check that the id would be the same for all traces
-            for tr in self:
+            for tr in self: # Use self to avoid creating temporary object
                 stats = tr.stats
-                if hasattr(self, "use_array") and self.use_array:
-                    test_str = ".".join(
-                        [
-                            stats.network,
-                            stats.station,
-                            stats.location,
-                            stats.channel,
-                        ],
-                    )
+                if use_array:
+                    test_str = ".".join([
+                        stats.network,
+                        stats.station,
+                        stats.location,
+                        stats.channel,
+                    ])
                 else:
-                    test_str = f"{stats.network}.{stats.station}.{stats.channel[0:2]}"
+                    test_str = ".".join([
+                        stats.network,
+                        stats.station,
+                        stats.location,
+                        stats.channel[0:2],
+                    ])
 
                 if id_str != test_str:
-                    raise ValueError("Inconsistent stream ID for different traces.")
+                    raise ValueError(f"Inconsistent stream IDs {id_str} and {test_str} for different traces.")
             self.id = id_str
 
     def get_id(self):
@@ -277,13 +283,12 @@ class StationStream(Stream):
         return self.id
 
     def get_net_sta(self):
-        """Get just the network and station component of the ID."""
-        return ".".join(self.get_id().split(".")[0:2])
+        """Get just the network and station components of the ID."""
+        return ".".join(self.id.split(".")[0:2])
 
     def get_net_sta_loc(self):
-        """Get network, station, and location codes."""
-        stats = self[0].stats
-        return ".".join([stats.network, stats.station, stats.location])
+        """Get just the network, station, and location components of the ID."""
+        return ".".join(self.id.split(".")[0:3])
 
     def get_inst(self):
         """Get first two characters of the channel code."""

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/core/stationtrace.py

@@ -185,10 +185,10 @@ class StationTrace(Trace):
                 header["coordinates"] = coords
                 header["standard"] = standard
                 header["format_specific"] = format_specific
-            except BaseException:
+            except BaseException as ex:
                 raise ValueError(
                     "Failed to construct required metadata from header data."
-                )
+                ) from ex
 
         # Sometimes the channel names do not indicate which one is the
         # Z channel. If we have vertical_orientation information, then
@@ -206,6 +206,7 @@ class StationTrace(Trace):
                 if is_vertcomp:
                     pass
                 elif is_horizcomp:
+                    # Cannot call self.warning() before setting parameters and may update header (before calling constructor)
                     logging.warning(
                         "Found vertical orientation in trace %s, which seems to be a horizontal component. Assuming vertical orientation is incorrect.",
                         trace_id,
@@ -665,6 +666,12 @@ class StationTrace(Trace):
             frequency_domain (bool):
                 Apply filter in frequency domain?
         """
+        if not frequency_domain and zerophase:
+            logging.warning(
+                "Filter is only applied once in frequency domain, "
+                "even if number of passes is 2"
+            )
+
         if zerophase:
             number_of_passes = 2
         else:
@@ -688,12 +695,6 @@ class StationTrace(Trace):
                 )
 
             else:
-                if zerophase:
-                    logging.warning(
-                        "Filter is only applied once in frequency domain, "
-                        "even if number of passes is 2"
-                    )
-
                 self.__compute_fft()
                 filter = np.sqrt(1.0 + (self.signal_freq / freq) ** (2.0 * corners))
                 self.__apply_filter(filter)
@@ -727,12 +728,6 @@ class StationTrace(Trace):
                 )
 
             else:
-                if zerophase:
-                    logging.warning(
-                        "Filter is only applied once in frequency domain, "
-                        "even if number of passes is 2"
-                    )
-
                 self.__compute_fft()
                 filter = np.sqrt(1.0 + (freq / self.signal_freq) ** (2.0 * corners))
                 self.__apply_filter(filter)
@@ -768,11 +763,6 @@ class StationTrace(Trace):
                 )
 
             else:
-                if zerophase:
-                    logging.warning(
-                        "Filter is only applied once in frequency domain, "
-                        "even if number of passes is 2"
-                    )
                 self.__compute_fft()
                 filter = np.sqrt(1.0 + (freqmin / self.signal_freq) ** (2.0 * corners))
                 filter *= np.sqrt(1.0 + (self.signal_freq / freqmax) ** (2.0 * corners))
@@ -809,12 +799,6 @@ class StationTrace(Trace):
                 )
 
             else:
-                if zerophase:
-                    logging.warning(
-                        "Filter is only applied once in frequency domain, "
-                        "even if number of passes is 2"
-                    )
-
                 self.__compute_fft()
                 filter = np.sqrt(1.0 + (freq / self.signal_freq) ** (2.0 * corners))
                 filter *= np.sqrt(1.0 + (self.signal_freq / freq) ** (2.0 * corners))

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/core/streamarray.py

@@ -123,11 +123,11 @@ class StreamArray(object):
         """Copy method."""
         return copy.deepcopy(self)
 
-    def select(self, network=None, station=None, instrument=None):
+    def select(self, network=None, station=None, location=None, instrument=None):
         """Select Streams.
 
         Return a new StreamCollection with only those StationStreams that
-        match network, station, and/or instrument selection criteria.
+        match network, station, location, and/or instrument selection criteria.
 
         Based on obspy's `select` method for traces.
 
@@ -136,6 +136,8 @@ class StreamArray(object):
                 Network code.
             station (str):
                 Station code.
+            location (str):
+                Location code.
             instrument (str):
                 Instrument code; i.e., the first two characters of the
                 channel.
@@ -143,15 +145,19 @@ class StreamArray(object):
         sel = []
         for st in self:
             inst = st.get_inst()
-            net_sta = st.get_net_sta()
-            net = net_sta.split(".")[0]
-            sta = net_sta.split(".")[1]
+            net_sta_loc = st.get_net_sta_loc().split(".")
+            net = net_sta_loc[0]
+            sta = net_sta_loc[1]
+            loc = net_sta_loc[2]
             if network is not None:
                 if not fnmatch.fnmatch(net.upper(), network.upper()):
                     continue
             if station is not None:
                 if not fnmatch.fnmatch(sta.upper(), station.upper()):
                     continue
+            if location is not None:
+                if not fnmatch.fnmatch(loc.upper(), location.upper()):
+                    continue
             if instrument is not None:
                 if not fnmatch.fnmatch(inst.upper(), instrument.upper()):
                     continue

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/core/streamcollection.py

@@ -114,7 +114,7 @@ class StreamCollection(StreamArray):
                     process_level_preference,
                     format_preference,
                 )
-        self.__group_by_net_sta_inst()
+        self.__group_by_net_sta_loc_inst()
         self.validate()
 
     def validate(self):
@@ -214,6 +214,9 @@ class StreamCollection(StreamArray):
 
                 # Loop over preferred instruments
                 no_match = True
+                group_tr = self[group[0]][0]
+                group_network = group_tr.stats.network
+                group_station = group_tr.stats.station
                 for pref in preference:
                     # Is this instrument available in the group?
                     r = re.compile(pref[0:2])
@@ -226,11 +229,12 @@ class StreamCollection(StreamArray):
 
                         # Label all non-selected streams in the group as failed
                         to_fail = group_insts
-                        to_fail.remove(keep)
+                        for inst in inst_match:
+                            to_fail.remove(inst)
                         for tf in to_fail:
                             for st in self.select(
-                                network=self[group[0]][0].stats.network,
-                                station=self[group[0]][0].stats.station,
+                                network=group_network,
+                                station=group_station,
                                 instrument=tf,
                             ):
                                 for tr in st:
@@ -304,7 +308,7 @@ class StreamCollection(StreamArray):
         stream_descript = self.describe_string()
         print(stream_descript)
 
-    def __group_by_net_sta_inst(self):
+    def __group_by_net_sta_loc_inst(self):
         trace_list = []
         stream_params = gather_stream_parameters(self.streams)
         for st in self.streams:
@@ -320,6 +324,7 @@ class StreamCollection(StreamArray):
             matches = [idx1]
             network = trace1.stats["network"]
             station = trace1.stats["station"]
+            location = trace1.stats["location"]
             free_field = trace1.free_field
             # For instrument, use first two characters of the channel
             inst = trace1.stats["channel"][0:2]
@@ -328,6 +333,7 @@ class StreamCollection(StreamArray):
                     if (
                         network == trace2.stats["network"]
                         and station == trace2.stats["station"]
+                        and location == trace2.stats["location"]
                         and inst == trace2.stats["channel"][0:2]
                         and free_field == trace2.free_field
                     ):
@@ -653,6 +659,7 @@ def are_duplicates(tr1, tr2, max_dist_tolerance):
         if (
             tr1.stats.station == tr2.stats.station
             and tr1.stats.channel[:2] == tr2.stats.channel[:2]
+            and tr1.stats.location == tr2.stats.location
             and len(orientation_codes) == 1
             and distance < max_dist_tolerance
         ):

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/data/config_production.yml

@@ -260,7 +260,6 @@ fetchers:
                 - 180.0
                 - -90.0
                 - 90.0
-                url: https://fdsnws.raspberryshakedata.com
             - SCEDC:
                 bounds:
                 - -160.0
@@ -620,7 +619,7 @@ integration:
 # -----------------------------------------------------------------------------
 # This section is for options for differentiation
 differentiation:
-    # Frequency or time domain differentation?
+    # Frequency or time domain differentiation?
     frequency: True
 
 # -----------------------------------------------------------------------------

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/data/config_production_v1.yml

@@ -19,7 +19,7 @@ fetchers:
             max_distance: 800.0
             # NOTE: the GMM is evaluated with a number of simplifications, including 
             # Ztor=0 and dip=90 deg such that the 'distance' is both Rrup and Rjb. See 
-            # the dynamic_searc_parameters.py module for further details. 
+            # the dynamic_search_parameters.py module for further details. 
     KNETFetcher:
         # Enable this fetcher?
         enabled: True
@@ -255,7 +255,6 @@ fetchers:
                 - 180.0
                 - -90.0
                 - 90.0
-                url: https://fdsnws.raspberryshakedata.com
             - SCEDC:
                 bounds:
                 - -160.0
@@ -623,7 +622,7 @@ integration:
 # -----------------------------------------------------------------------------
 # This section is for options for differentiation
 differentiation:
-    # Frequency or time domain differentation?
+    # Frequency or time domain differentiation?
     frequency: True
 
 # -----------------------------------------------------------------------------

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/asdf/flatfile.py

@@ -74,12 +74,6 @@ class Flatfile(object):
         """
         if hasattr(self.workspace, "config"):
             config = self.workspace.config
-            default_config_file = constants.DATA_DIR / constants.CONFIG_FILE_PRODUCTION
-            with open(default_config_file, "r", encoding="utf-8") as f:
-                yaml = YAML()
-                yaml.preserve_quotes = True
-                default_config = yaml.load(f)
-            update_dict(self.workspace.config, default_config)
         else:
             config = get_config()
         any_trace_failures = config["check_stream"]["any_trace_failures"]
@@ -95,13 +89,16 @@ class Flatfile(object):
 
             self._station_metrics = station_metrics
             passed_traces = [tr_met for tr_met in stream_meta if tr_met["passed"]]
+            n_total = len(stream_meta)
 
             n_passed = len(passed_traces)
+            n_failed = n_total - n_passed
+
             for tr in passed_traces:
                 tr.update({"use_array": use_array})
 
             if any_trace_failures:
-                if n_passed < 3:
+                if n_failed > 0:
                     stream_passed = False
                 else:
                     stream_passed = True
@@ -124,7 +121,7 @@ class Flatfile(object):
                 row = self.get_imt_row(passed_traces, wml, imc)
                 if not row:
                     continue
-                self.imc_tables[str(imc)].append(row)  #
+                self.imc_tables[str(imc)].append(row)
 
     def get_imt_row(self, passed_traces, wml, imc):
         """Get one row of the IMT flatfile fom a stream.

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/asdf/stream_workspace.py

@@ -418,7 +418,7 @@ class StreamWorkspace(object):
         )
 
         for stream in streams:
-            logging.info("Adding waveforms for station %s", stream.get_id())
+            logging.info("Adding waveforms for %s", stream.get_id())
             # is this a raw file? Check the trace for provenance info.
             is_raw = not stream[0].provenance.ids
 

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/cosmos/cesmd_fetcher.py

@@ -1,6 +1,6 @@
 # stdlib imports
 import logging
-from datetime import datetime, timedelta
+from datetime import timedelta
 
 # third party imports
 import numpy as np
@@ -172,7 +172,9 @@ class CESMDFetcher(DataFetcher):
         self.lon = lon
         self.radius = radius
         self.station_radius = station_radius
-        self.eq_dt = eq_dt
+        self.dt = eq_dt
+        self.dmag = DMAG
+        self.ddepth = DDEPTH
         self.rawdir = rawdir
         self.depth = depth
         self.magnitude = magnitude
@@ -202,7 +204,7 @@ class CESMDFetcher(DataFetcher):
                 eqtime=self.time,
                 abandoned=False,
                 station_type=self.station_type,
-                eqtimewindow=self.eq_dt,  # seconds
+                eqtimewindow=self.dt,  # seconds
                 eqradius=self.radius,  # km
                 station_radius=self.station_radius,
             )
@@ -220,7 +222,7 @@ class CESMDFetcher(DataFetcher):
             if event["id"] in idlist:
                 continue
             idlist.append(event["id"])
-            tevent["time"] = datetime.strptime(event["time"], "%Y-%m-%d %H:%M:%S")
+            tevent["time"] = UTCDateTime(event["time"])
             tevent["lat"] = event["latitude"]
             tevent["lon"] = event["longitude"]
             tevent["depth"] = event["depth"]
@@ -259,8 +261,8 @@ class CESMDFetcher(DataFetcher):
         if not has_event:
             raise Exception("Could not find matching event.")
 
-        starttime = self.time - timedelta(seconds=self.eq_dt // 2)
-        endtime = self.time + timedelta(seconds=self.eq_dt // 2)
+        starttime = self.time - timedelta(seconds=self.dt // 2)
+        endtime = self.time + timedelta(seconds=self.dt // 2)
 
         if nstations < MAX_STATIONS:
             try:

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/cosmos/cosmos_writer.py

@@ -930,8 +930,6 @@ class CosmosWriter(object):
                         sta = trace.stats.station
                         cha = trace.stats.channel
                         loc = trace.stats.location
-                        # COSMOS files will have trace.stats.format_specific,
-                        # miniseed/stationxml will not.
                         if trace.stats.format_specific:
                             trace.stats["stage_1_sensitivity"] = (
                                 trace.stats.format_specific["stage_1_sensitivity"]
@@ -939,26 +937,6 @@ class CosmosWriter(object):
                             trace.stats["data_logger_sensitivity"] = (
                                 trace.stats.format_specific["data_logger_sensitivity"]
                             )
-                        else:
-                            # Do somthing with StationXML information
-                            tinv = inv.select(
-                                network=trace.stats.network,
-                                station=trace.stats.station,
-                                channel=trace.stats.channel,
-                                location=trace.stats.location,
-                                time=trace.stats.starttime,
-                            )
-                            stages = tinv[0][0][0].response.response_stages
-                            instrument_sensitivity, data_logger_sensitivity = 1, 1
-                            for idx, stage in enumerate(stages):
-                                if idx == 0:
-                                    instrument_sensitivity = stage.stage_gain
-                                else:
-                                    data_logger_sensitivity *= stage.stage_gain
-                            trace.stats["stage_1_sensitivity"] = instrument_sensitivity
-                            trace.stats["data_logger_sensitivity"] = (
-                                data_logger_sensitivity
-                            )
 
                         if trace.stats.standard.units_type != "acc":
                             msg = (

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/obspy/core.py

@@ -7,6 +7,7 @@ import os
 import re
 
 import numpy as np
+import scipy.constants as sp
 from obspy import read_inventory
 from obspy.core.stream import read
 
@@ -176,19 +177,20 @@ def read_obspy(filename, config=None, **kwargs):
         network = ttrace.stats.network
         station = ttrace.stats.station
         channel = ttrace.stats.channel
-        nsc = f"{network}.{station}.{channel}"
+        location = ttrace.stats.location
+        nslc = f"{network}.{station}.{channel}"
 
         if inventory:
             tinventory = inventory.select(
-                network=ttrace.stats.network,
-                station=ttrace.stats.station,
-                channel=ttrace.stats.channel,
-                location=ttrace.stats.location,
+                network=network,
+                station=station,
+                channel=channel,
+                location=location,
                 time=ttrace.stats.starttime,
             )
             if tinventory[0][0][0].sample_rate != ttrace.stats.sampling_rate:
                 raise ValueError(
-                    f"Sample rate is not consistent betwen data and metadata ({nsc})"
+                    f"Sample rate is not consistent betwen data and metadata ({nslc})"
                 )
         else:
             tinventory = None
@@ -197,7 +199,10 @@ def read_obspy(filename, config=None, **kwargs):
             data=ttrace.data, header=ttrace.stats, inventory=tinventory, config=config
         )
 
-        if not np.issubdtype(trace.data.dtype, np.integer):
+        if not (
+            np.issubdtype(trace.data.dtype, np.integer)
+            or np.allclose(trace.data % 1, 0)
+        ):
             trace.fail("Raw data type is not integer.")
 
         if ttrace.stats.location == "":
@@ -254,6 +259,26 @@ def read_obspy(filename, config=None, **kwargs):
             # Apply conversion factor if one was specified for this format
             trace.data *= float(config["read"]["sac_conversion_factor"])
 
+        if "_format" in trace.stats and trace.stats._format.lower() == "mseed":
+            # Unit conversions used by the COSMOS writer for real hdr 22 and 42
+            # This code assumes that specific stages in the instrument response
+            # correspond to the sensor sensitivity and the data logger sensitivity.
+
+            stages = tinventory[0][0][0].response.response_stages
+            data_logger_sensitivity = 1
+            for idx, stage in enumerate(stages):
+                if idx == 0:
+                    # Sensor sensitivity (row 9 column 2)
+                    trace.stats["format_specific"]["stage_1_sensitivity"] = (
+                        float(stage.stage_gain)
+                    ) * sp.g  # Units will be in (units per motion) / g
+                else:
+                    data_logger_sensitivity *= stage.stage_gain
+            # Data logger sensitivity (row 5 column 2)
+            trace.stats["format_specific"]["data_logger_sensitivity"] = (
+                1 / float(data_logger_sensitivity)
+            ) * 1e6  # It's the inverse bit-weight in microvolts
+
         traces.append(trace)
     if no_match:
         stream = StationStream(traces=traces, config=config)

{gmprocess-2.4.0 → gmprocess-2.5.0}/src/gmprocess/io/unam/core.py

@@ -218,9 +218,9 @@ def _read_header(filename):
     channels[2]["npts"] = npts_list[2]
 
     # locations
-    channels[0]["location"] = "--"
-    channels[1]["location"] = "--"
-    channels[1]["location"] = "--"
+    channels[0]["location"] = ""
+    channels[1]["location"] = ""
+    channels[2]["location"] = ""
 
     # get station coordinates
     coord1 = header_dict["COORDENADAS DE LA ESTACION"]