requake 0.4.1__tar.gz → 0.6__tar.gz

{requake-0.4.1 → requake-0.6}/CHANGELOG.md

@@ -1,6 +1,35 @@
-# Requake - Repeating earthquakes search and analysis
+# Requake Changelog
 
-(c) 2021-2024 Claudio Satriano <satriano@ipgp.fr>
+Copyright (c) 2021-2024 Claudio Satriano <satriano@ipgp.fr>
+
+## v0.6 - 2024-05-04
+
+Note: you might want to run `requake update_config` to update your config file
+to the latest version.
+
+- Verb `plot_slip` renamed to `plot_cumulative`. This new verb has new options
+  to plot cumulative slip, cumulative moment, and cumulative number of events,
+  and to make the plot logarithmic.
+- New verb `print_catalog`: print the event catalog to screen
+- New verb `print_pairs`: print the event pairs to screen
+- `print_families`: also print minimum and maximum family magnitudes
+- Additional models to convert magnitude to slip. Currently supported models
+  are: Nadeau and Johnson (1998), Beeler et al. (2001), Eshelby (1957).
+  Model selection is done using the `mag_to_slip_model` config parameter.
+- New verb `update_config`: update an existing config file to the latest
+  version
+
+## v0.5 - 2024-04-23
+
+- New config options: `station_metadata_path` and `waveform_data_path` to
+  read station metadata and waveform data from files. Supports any metadata
+  format supported by ObsPy and SDS (SeisComp Data Structure) waveform
+  archives.
+- Filter catalog files on reading using the criteria in the config file
+- Improved time axes in `plot_timespans` and `plot_slip` for short time
+  intervals
+- `print_families`: autoset duration units based on the average duration
+- `plot_timespans`: default sorting changed to `family_number`
 
 ## v0.4.1 - 2024-03-11
 

{requake-0.4.1 → requake-0.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: requake
-Version: 0.4.1
+Version: 0.6
 Summary: Repeating earthquakes search and analysis
 Home-page: https://requake.seismicsource.org
 Author: Claudio Satriano
@@ -33,6 +33,8 @@ Requires-Dist: argcomplete
 Requires-Dist: tqdm
 Requires-Dist: tabulate
 
+<img src="https://cdn.jsdelivr.net/gh/SeismicSource/requake@dc9a6e6fd818e57e0c0d74f9c16f2bfc6d4d148c/imgs/Requake_logo.svg" width="400">
+
 # Requake
 
 Repeating earthquakes search and analysis.
@@ -40,6 +42,8 @@ Repeating earthquakes search and analysis.
 [![changelog-badge]][changelog-link]
 [![PyPI-badge]][PyPI-link]
 [![license-badge]][license-link]
+[![docs-badge]][docs-link]
+[![DOI-badge]][DOI-link]
 
 Copyright (c) 2021-2024 Claudio Satriano <satriano@ipgp.fr>
 
@@ -50,8 +54,8 @@ Requake is a command line tool to search and analyse repeating earthquakes.
 It can either scan an existing earthquake catalog to search for similar events,
 or perform template matching on a continuous waveform stream.
 
-Catalogs and waveforms are downloaded using standard
-[FDSN web services](https://www.fdsn.org/webservices/).
+Catalogs and waveforms can be read from local files or downloaded using
+standard [FDSN web services](https://www.fdsn.org/webservices/).
 
 Requake is written in Python and uses [ObsPy](https://obspy.org) as backend.
 
@@ -118,33 +122,34 @@ To get help, use:
 
     requake -h
 
-The different running modes are specified as "verbs" (positional arguments).
-Currently supported verbs are:
+Different commands are available:
 
     sample_config       write sample config file to current directory and exit
+    update_config       update an existing config file to the latest version
     read_catalog        read an event catalog from web services or from a file
+    print_catalog       print the event catalog to screen
     scan_catalog        scan an existing catalog for earthquake pairs
+    print_pairs         print pairs to screen
     plot_pair           plot traces for a given event pair
     build_families      build families of repeating earthquakes from a catalog
                         of pairs
     print_families      print families to screen
     plot_families       plot traces for one ore more event families
     plot_timespans      plot family timespans
-    plot_slip           plot cumulative slip for one or more families
+    plot_cumulative     cumulative plot for one or more families
     map_families        plot families on a map
     flag_family         flag a family of repeating earthquakes as valid or not
-                        valid. Note that all families are valid by default
-                        when first created
+                        valid.
     build_templates     build waveform templates for one or more event
                         families
     scan_templates      scan a continuous waveform stream using one or more
                         templates
 
-Certain running modes (e.g., `plot_pair`) require further arguments (use, e.g.,
+Certain commands (e.g., `plot_pair`) require further arguments (use, e.g.,
 `requake plot_pair -h` to get help).
 
-Requake supports command line tab completion for arguments, thanks to
-[argcomplete](https://kislyuk.github.io/argcomplete/).
+Requake supports command line tab completion for commands and arguments, thanks
+to [argcomplete](https://kislyuk.github.io/argcomplete/).
 To enable command line tab completion, add the following line to your `.bashrc`
 or `.zshrc`:
 
@@ -184,10 +189,35 @@ parallelization will be implemented.
 
 - `requake build_families` is fast™.
 
+## How to Cite
+
+If you used Requake for a scientific paper, please cite it as:
+
+> Satriano, C. (2024). Requake: Repeating earthquakes search and analysis (X.Y).
+> [doi: 10.5281/ZENODO.10832204]
+
+Please replace `X.Y` with the Requake version number you used.
+
+You can also cite the following abstract presented at the
+2016 AGU Fall Meeting:
+
+> Satriano, C., Doucet, A. & Bouin, M.-P. (2021).
+> Probing the creep rate along the Lesser Antilles Subduction Zone through repeating earthquakes.
+> In AGU Fall Meeting Abstracts
+> (Vol. 2021, pp. T25A-0167), [bibcode: 2021AGUFM.T25A0167S]
+
 <!-- Badges and project links -->
 [PyPI-badge]: http://img.shields.io/pypi/v/requake.svg
 [PyPI-link]: https://pypi.python.org/pypi/requake
 [license-badge]: https://img.shields.io/badge/license-GPLv3-green
 [license-link]: https://www.gnu.org/licenses/gpl-3.0.html
+[docs-badge]: https://readthedocs.org/projects/requake/badge/?version=latest
+[docs-link]: https://requake.readthedocs.io/en/latest/?badge=latest
 [changelog-badge]: https://img.shields.io/badge/Changelog-136CB6.svg
 [changelog-link]: https://github.com/SeismicSource/requake/blob/main/CHANGELOG.md
+[DOI-badge]: https://zenodo.org/badge/DOI/10.5281/zenodo.10832204.svg
+[DOI-link]: https://doi.org/10.5281/zenodo.10832204
+
+<!-- References -->
+[doi: 10.5281/ZENODO.10832204]: https://doi.org/10.5281/ZENODO.10832204
+[bibcode: 2021AGUFM.T25A0167S]: https://ui.adsabs.harvard.edu/abs/2021AGUFM.T25A0167S

{requake-0.4.1 → requake-0.6}/README.md

@@ -1,3 +1,5 @@
+<img src="imgs/Requake_logo.svg" width="400">
+
 # Requake
 
 Repeating earthquakes search and analysis.
@@ -5,6 +7,8 @@ Repeating earthquakes search and analysis.
 [![changelog-badge]][changelog-link]
 [![PyPI-badge]][PyPI-link]
 [![license-badge]][license-link]
+[![docs-badge]][docs-link]
+[![DOI-badge]][DOI-link]
 
 Copyright (c) 2021-2024 Claudio Satriano <satriano@ipgp.fr>
 
@@ -15,8 +19,8 @@ Requake is a command line tool to search and analyse repeating earthquakes.
 It can either scan an existing earthquake catalog to search for similar events,
 or perform template matching on a continuous waveform stream.
 
-Catalogs and waveforms are downloaded using standard
-[FDSN web services](https://www.fdsn.org/webservices/).
+Catalogs and waveforms can be read from local files or downloaded using
+standard [FDSN web services](https://www.fdsn.org/webservices/).
 
 Requake is written in Python and uses [ObsPy](https://obspy.org) as backend.
 
@@ -83,33 +87,34 @@ To get help, use:
 
     requake -h
 
-The different running modes are specified as "verbs" (positional arguments).
-Currently supported verbs are:
+Different commands are available:
 
     sample_config       write sample config file to current directory and exit
+    update_config       update an existing config file to the latest version
     read_catalog        read an event catalog from web services or from a file
+    print_catalog       print the event catalog to screen
     scan_catalog        scan an existing catalog for earthquake pairs
+    print_pairs         print pairs to screen
     plot_pair           plot traces for a given event pair
     build_families      build families of repeating earthquakes from a catalog
                         of pairs
     print_families      print families to screen
     plot_families       plot traces for one ore more event families
     plot_timespans      plot family timespans
-    plot_slip           plot cumulative slip for one or more families
+    plot_cumulative     cumulative plot for one or more families
     map_families        plot families on a map
     flag_family         flag a family of repeating earthquakes as valid or not
-                        valid. Note that all families are valid by default
-                        when first created
+                        valid.
     build_templates     build waveform templates for one or more event
                         families
     scan_templates      scan a continuous waveform stream using one or more
                         templates
 
-Certain running modes (e.g., `plot_pair`) require further arguments (use, e.g.,
+Certain commands (e.g., `plot_pair`) require further arguments (use, e.g.,
 `requake plot_pair -h` to get help).
 
-Requake supports command line tab completion for arguments, thanks to
-[argcomplete](https://kislyuk.github.io/argcomplete/).
+Requake supports command line tab completion for commands and arguments, thanks
+to [argcomplete](https://kislyuk.github.io/argcomplete/).
 To enable command line tab completion, add the following line to your `.bashrc`
 or `.zshrc`:
 
@@ -149,10 +154,35 @@ parallelization will be implemented.
 
 - `requake build_families` is fast™.
 
+## How to Cite
+
+If you used Requake for a scientific paper, please cite it as:
+
+> Satriano, C. (2024). Requake: Repeating earthquakes search and analysis (X.Y).
+> [doi: 10.5281/ZENODO.10832204]
+
+Please replace `X.Y` with the Requake version number you used.
+
+You can also cite the following abstract presented at the
+2016 AGU Fall Meeting:
+
+> Satriano, C., Doucet, A. & Bouin, M.-P. (2021).
+> Probing the creep rate along the Lesser Antilles Subduction Zone through repeating earthquakes.
+> In AGU Fall Meeting Abstracts
+> (Vol. 2021, pp. T25A-0167), [bibcode: 2021AGUFM.T25A0167S]
+
 <!-- Badges and project links -->
 [PyPI-badge]: http://img.shields.io/pypi/v/requake.svg
 [PyPI-link]: https://pypi.python.org/pypi/requake
 [license-badge]: https://img.shields.io/badge/license-GPLv3-green
 [license-link]: https://www.gnu.org/licenses/gpl-3.0.html
+[docs-badge]: https://readthedocs.org/projects/requake/badge/?version=latest
+[docs-link]: https://requake.readthedocs.io/en/latest/?badge=latest
 [changelog-badge]: https://img.shields.io/badge/Changelog-136CB6.svg
-[changelog-link]: https://github.com/SeismicSource/requake/blob/main/CHANGELOG.md
+[changelog-link]: https://github.com/SeismicSource/requake/blob/main/CHANGELOG.md
+[DOI-badge]: https://zenodo.org/badge/DOI/10.5281/zenodo.10832204.svg
+[DOI-link]: https://doi.org/10.5281/zenodo.10832204
+
+<!-- References -->
+[doi: 10.5281/ZENODO.10832204]: https://doi.org/10.5281/ZENODO.10832204
+[bibcode: 2021AGUFM.T25A0167S]: https://ui.adsabs.harvard.edu/abs/2021AGUFM.T25A0167S

{requake-0.4.1 → requake-0.6}/requake/_version.py

@@ -8,11 +8,11 @@ import json
 
 version_json = '''
 {
- "date": "2024-03-11T15:31:32+0100",
+ "date": "2024-05-04T17:23:20+0200",
  "dirty": false,
  "error": null,
- "full-revisionid": "fa843d40d788a71f300e374e489cb575396d71c6",
- "version": "0.4.1"
+ "full-revisionid": "dc9a6e6fd818e57e0c0d74f9c16f2bfc6d4d148c",
+ "version": "0.6"
 }
 '''  # END VERSION_JSON
 

{requake-0.4.1 → requake-0.6}/requake/catalog/__init__.py

@@ -10,3 +10,4 @@ Functions and data structures for managing earthquake catalogs.
     (https://www.gnu.org/licenses/gpl-3.0-standalone.html)
 """
 from .read_catalog import read_catalog  # noqa
+from .print_catalog import print_catalog  # noqa

{requake-0.4.1 → requake-0.6}/requake/catalog/catalog.py

@@ -182,6 +182,72 @@ class RequakeCatalog(list):
                 fp.write(ev.fdsn_text() + '\n')
 
 
+    def filter(
+            self, starttime=None, endtime=None,
+            minlatitude=None, maxlatitude=None,
+            minlongitude=None, maxlongitude=None,
+            mindepth=None, maxdepth=None,
+            minmagnitude=None, maxmagnitude=None):
+        """
+        Filter the catalog, based on the specified criteria.
+
+        Returns a new catalog.
+
+        :param starttime: start time
+        :type starttime: obspy.UTCDateTime
+        :param endtime: end time
+        :type endtime: obspy.UTCDateTime
+        :param minlatitude: minimum latitude
+        :type minlatitude: float
+        :param maxlatitude: maximum latitude
+        :type maxlatitude: float
+        :param minlongitude: minimum longitude
+        :type minlongitude: float
+        :param maxlongitude: maximum longitude
+        :type maxlongitude: float
+        :param mindepth: minimum depth (in km)
+        :type mindepth: float
+        :param maxdepth: maximum depth (in km)
+        :type maxdepth: float
+        :param minmagnitude: minimum magnitude
+        :type minmagnitude: float
+        :param maxmagnitude: maximum magnitude
+        :type maxmagnitude: float
+
+        :return: a new catalog
+        :rtype: RequakeCatalog
+        """
+        outcat = RequakeCatalog()
+        for ev in self:
+            if ev.orig_time is not None:
+                if starttime is not None and ev.orig_time < starttime:
+                    continue
+                if endtime is not None and ev.orig_time > endtime:
+                    continue
+            if ev.lat is not None:
+                if minlatitude is not None and ev.lat < minlatitude:
+                    continue
+                if maxlatitude is not None and ev.lat > maxlatitude:
+                    continue
+            if ev.lon is not None:
+                if minlongitude is not None and ev.lon < minlongitude:
+                    continue
+                if maxlongitude is not None and ev.lon > maxlongitude:
+                    continue
+            if ev.depth is not None:
+                if mindepth is not None and ev.depth < mindepth:
+                    continue
+                if maxdepth is not None and ev.depth > maxdepth:
+                    continue
+            if ev.mag is not None:
+                if minmagnitude is not None and ev.mag < minmagnitude:
+                    continue
+                if maxmagnitude is not None and ev.mag > maxmagnitude:
+                    continue
+            outcat.append(ev)
+        return outcat
+
+
 def read_stored_catalog(config):
     """
     Read the catalog stored in the output directory.
@@ -198,10 +264,13 @@ def read_stored_catalog(config):
     try:
         cat = RequakeCatalog()
         cat.read(config.scan_catalog_file)
+        cat.sort()
+        if not cat:
+            raise ValueError('Empty catalog')
         return cat
     except ValueError as m:
         raise ValueError(
-            f'Error reading catalog file {config.scan_catalog_file}'
+            f'Error reading catalog file {config.scan_catalog_file}: {m}'
         ) from m
     except FileNotFoundError as m:
         raise FileNotFoundError(

requake-0.6/requake/catalog/print_catalog.py

@@ -0,0 +1,53 @@
+# -*- coding: utf8 -*-
+# SPDX-License-Identifier: GPL-3.0-or-later
+"""
+Print the event catalog to screen.
+
+:copyright:
+    2021-2024 Claudio Satriano <satriano@ipgp.fr>
+:license:
+    GNU General Public License v3.0 or later
+    (https://www.gnu.org/licenses/gpl-3.0-standalone.html)
+"""
+import logging
+from .catalog import read_stored_catalog
+from ..config.generic_printer import generic_printer
+from ..config.rq_setup import rq_exit
+logger = logging.getLogger(__name__.rsplit('.', maxsplit=1)[-1])
+
+
+def print_catalog(config):
+    """
+    Print the event catalog to screen.
+
+    :param config: Configuration object.
+    :type config: config.Config
+    """
+    try:
+        catalog = read_stored_catalog(config)
+    except (ValueError, FileNotFoundError) as m:
+        logger.error(m)
+        rq_exit(1)
+
+    headers_fmt = [
+        ('evid', None),
+        ('origin time', None),
+        ('longitude', '.4f'),
+        ('latitude', '.4f'),
+        ('depth\n(km)', '.3f'),
+        ('mag\ntype', None),
+        ('mag', '.1f')
+    ]
+    rows = [
+        [
+            ev.evid,
+            ev.orig_time,
+            ev.lon,
+            ev.lat,
+            ev.depth,
+            ev.mag_type,
+            ev.mag
+        ]
+        for ev in catalog
+    ]
+    generic_printer(config, rows, headers_fmt)

{requake-0.4.1 → requake-0.6}/requake/catalog/read_catalog.py

@@ -9,6 +9,7 @@ Read an event catalog from web services or from a file.
     GNU General Public License v3.0 or later
     (https://www.gnu.org/licenses/gpl-3.0-standalone.html)
 """
+import os
 import logging
 import contextlib
 from ..catalog.catalog import RequakeCatalog
@@ -34,8 +35,14 @@ def _read_catalog_from_file(config):
     :raises ValueError: if the file format is not supported
     """
     catalog_file = config.args.catalog_file
+    # return ValueError if the file is empty
+    if os.stat(catalog_file).st_size == 0:
+        raise ValueError('Empty file')
+    # return ValueError if the file is a directory
+    if os.path.isdir(catalog_file):
+        raise ValueError('Is a directory')
     # try to read the catalog as a QuakeML file
-    with contextlib.suppress(TypeError):
+    with contextlib.suppress(TypeError, IndexError, ValueError):
         return read_catalog_from_quakeml(catalog_file)
     # try to read the catalog as a FDSN text file
     with contextlib.suppress(ValueError):
@@ -47,6 +54,45 @@ def _read_catalog_from_file(config):
     return read_catalog_from_csv(catalog_file)
 
 
+def _filter_catalog(catalog, config):
+    """
+    Filter an event catalog, based on the criteria specified
+    in the configuration.
+
+    :param catalog: Event catalog.
+    :type catalog: requake.catalog.RequakeCatalog
+    :param config: Configuration object.
+    :type config: requake.rq_setup.RequakeConfig
+    :return: Filtered event catalog.
+    :rtype: requake.catalog.RequakeCatalog
+    """
+    lon_min = config.catalog_lon_min
+    lon_max = config.catalog_lon_max
+    lat_min = config.catalog_lat_min
+    lat_max = config.catalog_lat_max
+    depth_min = config.catalog_depth_min
+    depth_max = config.catalog_depth_max
+    mag_min = config.catalog_mag_min
+    mag_max = config.catalog_mag_max
+    outcat = catalog
+    for start_time, end_time in zip(config.catalog_start_times,
+                                    config.catalog_end_times):
+        outcat = outcat.filter(
+            starttime=start_time, endtime=end_time,
+            minlatitude=lat_min, maxlatitude=lat_max,
+            minlongitude=lon_min, maxlongitude=lon_max,
+            mindepth=depth_min, maxdepth=depth_max,
+            minmagnitude=mag_min, maxmagnitude=mag_max
+        )
+    nevents_in = len(catalog)
+    nevents_out = len(outcat)
+    logger.info(
+        'Catalog filtered based on configuration file: '
+        f'{nevents_out} events out of {nevents_in}'
+    )
+    return outcat
+
+
 def read_catalog(config):
     """
     Read an event catalog from web services or from a file.
@@ -59,34 +105,42 @@ def read_catalog(config):
     :rtype: requake.catalog.RequakeCatalog
     """
     catalog = RequakeCatalog()
-    out_cat_file = config.scan_catalog_file
+    output_cat_file = config.scan_catalog_file
     nevs_read = 0
     if config.args.append:
         with contextlib.suppress(FileNotFoundError):
             catalog = RequakeCatalog()
-            catalog.read(out_cat_file)
+            catalog.read(output_cat_file)
             nevs_read = len(catalog)
-            logger.info(f'{nevs_read} events read from "{out_cat_file}"')
+            logger.info(f'{nevs_read} events read from "{output_cat_file}"')
     logger.info('Reading catalog...')
-    in_cat_file = config.args.catalog_file
-    if in_cat_file is not None:
+    input_cat_file = config.args.catalog_file
+    if input_cat_file is not None:
         try:
             catalog += _read_catalog_from_file(config)
+            # Filter catalog based on configuration
+            catalog = _filter_catalog(catalog, config)
         except FileNotFoundError:
-            logger.error(f'File "{in_cat_file}" not found')
+            logger.error(f'File "{input_cat_file}" not found')
             rq_exit(1)
         except ValueError as m:
-            logger.error(f'Error reading catalog file "{in_cat_file}": {m}')
+            logger.error(f'Error reading catalog file "{input_cat_file}": {m}')
             rq_exit(1)
     else:
         catalog += read_catalog_from_fdsnws(config)
     if not catalog:
         logger.error('No event read')
         rq_exit(1)
+    # Deduplicate catalog
+    len_before = len(catalog)
     catalog.deduplicate()
+    len_after = len(catalog)
+    nevs_dedup = len_before - len_after
+    if nevs_dedup > 0:
+        logger.info(f'{nevs_dedup} duplicate events removed')
     # Sort catalog in increasing time order
     catalog.sort()
     # Write catalog to output file
-    catalog.write(out_cat_file)
+    catalog.write(output_cat_file)
     nevs_written = len(catalog) - nevs_read
-    logger.info(f'{nevs_written} events written to "{out_cat_file}"')
+    logger.info(f'{nevs_written} events written to "{output_cat_file}"')

{requake-0.4.1 → requake-0.6}/requake/catalog/read_catalog_from_csv.py

@@ -183,7 +183,17 @@ def read_catalog_from_csv(filename):
                     year=year, month=month, day=day,
                     hour=hour, minute=minute, second=0) + seconds
             else:
-                orig_time = UTCDateTime(row[fields['orig_time']])
+                try:
+                    orig_time = UTCDateTime(row[fields['orig_time']])
+                except ValueError:
+                    # one last try: check if the time is in the format
+                    # YYYYMMDD.hhmmss.
+                    # Replace the dot with a space, pad with zeros
+                    # and try again
+                    orig_time = UTCDateTime(
+                        row[fields['orig_time']]
+                        .replace('.', ' ').ljust(15, '0')
+                    )
             row[None] = None
             ev = RequakeEvent()
             ev.orig_time = orig_time

{requake-0.4.1 → requake-0.6}/requake/config/configspec.conf

@@ -4,9 +4,19 @@
 ## FDSN station webservice URL to retrieve station metadata.
 ## Use a full URL or a ObsPy shortcut (e.g., IRIS, RESIF, GEOFON, ...)
 fdsn_station_url = string(default=http://example.com)
+## Alternatively, you can provide the path to a local file with station
+## metadata in any format supported by ObsPy (note that some formats may not
+## contain station coordinates, which are required for the scan,
+## see https://docs.obspy.org/packages/obspy.core.inventory.html).
+## If both fdsn_station_url and stationxml_file are provided, the
+## local file will be used.
+station_metadata_path = string(default=None)
 ## FDSN dataselect webservice URL to retrieve waveforms.
 ## Use a full URL or a ObsPy shortcut (e.g., IRIS, RESIF, GEOFON, ...)
 fdsn_dataselect_url = string(default=http://example.com)
+## Alternatively, you can provide the path to a local SDS waveform archive
+## (see https://docs.obspy.org/packages/autogen/obspy.clients.filesystem.sds.html)
+waveform_data_path = string(default=None)
 
 #### Catalog-based scan
 ### The following parameters are for a catalog-based scan:
@@ -118,3 +128,19 @@ distance_from_lat = float(default=None)
 ## Normalization will produce an average trace which is more representative
 ## of the family, but will be less effective in reducing noise.
 normalize_traces_before_averaging = boolean(default=False)
+
+#### Post-processing parameters
+## Estimation of fault slip (in cm) from repetater's magnitude
+## Choose between:
+##  NJ1998 : Nadeau and Johnson (1998)
+##  B2001 : Beeler et al. (2001)
+##  E1957 : Eshelby (1957), circular crack model
+## See Uchida (2019, https://doi.org/10.1186/s40645-019-0284-z) for details
+## on the models.
+mag_to_slip_model = string(default='NJ1998')
+## Static stress drop value for B2001 and E1957 models (MPa)
+static_stress_drop = float(default=10.0)
+## Rigidity value for B2001 and E1957 models (GPa)
+rigidity = float(default=30.0)
+## Strain hardening coefficient for B2001 model (MPa/cm)
+strain_hardening = float(default=0.5)

requake-0.6/requake/config/generic_printer.py

@@ -0,0 +1,66 @@
+# -*- coding: utf8 -*-
+# SPDX-License-Identifier: GPL-3.0-or-later
+"""
+Generic printer functions for Requake.
+
+:copyright:
+    2021-2024 Claudio Satriano <satriano@ipgp.fr>
+:license:
+    GNU General Public License v3.0 or later
+    (https://www.gnu.org/licenses/gpl-3.0-standalone.html)
+"""
+import sys
+import os
+import csv
+import contextlib
+import logging
+from tabulate import tabulate
+logger = logging.getLogger(__name__.rsplit('.', maxsplit=1)[-1])
+
+
+def generic_printer(config, rows, headers_fmt, print_headers=True):
+    """
+    A generic printer function for Requake.
+
+    :param config: Configuration object.
+    :type config: config.Config
+    :param rows: Rows to print.
+    :type rows: list of rows
+    :param headers_fmt: Headers and format strings.
+    :type headers_fmt: list of tuples of str
+    """
+    headers = [h[0] for h in headers_fmt]
+    floatfmt = [h[1] for h in headers_fmt]
+    tablefmt = config.args.format
+    if tablefmt == 'csv':
+        writer = csv.writer(sys.stdout)
+        if print_headers:
+            # replace newlines with spaces in headers
+            headers = [h.replace('\n', ' ') for h in headers]
+            writer.writerow(headers)
+    elif tablefmt == 'markdown':
+        headers = [h.replace('\n', '<br>') for h in headers]
+    if tablefmt == 'csv':
+        try:
+            writer.writerows(rows)
+        except BrokenPipeError:
+            # Redirect remaining output to devnull to avoid another
+            # BrokenPipeError at shutdown
+            devnull = os.open(os.devnull, os.O_WRONLY)
+            os.dup2(devnull, sys.stdout.fileno())
+    else:
+        format_dict = {
+            'simple': 'simple',
+            'markdown': 'github'
+        }
+        tablefmt = format_dict[config.args.format]
+        with contextlib.suppress(BrokenPipeError):
+            kwargs = {
+                'headers': headers,
+                'floatfmt': floatfmt,
+                'tablefmt': tablefmt
+            }
+            line = tabulate(rows, **kwargs)
+            if not print_headers:
+                line = line.rsplit('\n', maxsplit=1)[-1]
+            print(line)