tectonic-utils 0.1.2__py3-none-any.whl

tectonic_utils/seismo/second_focal_plane.py

@@ -0,0 +1,138 @@
+"""
+Utility tool for computing the second possible strike, rake, and dip in a focal mechanism,
+given the strike, rake, and dip of the first focal plane.
+Inspiration from Vince Cronin's website.
+"""
+
+# Internal information:
+# Variable "a" is the reference strike azimuth.
+# Variable "b" is the unit vector toward the reference strike azimuth in
+# a coordinate system in which the +Y axis is north, the +X axis is east, and the +Z axis is up.
+# Variable "c" is the unit dip vector.
+# Variable "d" is the vector normal to the input nodal plane, and is the slip vector for the auxilary plane.
+# Vector "d" points up if the input rake is a positive number, and down if the input rake is a negative number.
+# Variable "f" is the basis vector for the X axis of the geographic coordinate system, which points east.
+# Variable "g" is the basis vector for the Y axis of the geographic coordinate system, which points north.
+# Variable "h" is the basis vector for the Z axis of the geographic coordinate system, which points up.
+# Variable "i" is the basis vector for the X axis of the nodal plane coordinate system.
+# If the rake is a negative number, the strike vector is the X axis and the dip vector is the Y axis.
+# If the rake is a positive vector, the dip vector is the X axis and the strike vector is the Y axis.
+# Variable "j" is the basis vec  tor for the Y axis of the nodal plane coordinate system.
+# Variable "d" is the basis vector of the Z axis of the nodal plane coordinate system.
+# Variable "k" is the unit vector in the nodal plane coordinate system that is coincident with the slip vector
+# on the (input) nodal plane.
+# Variable "m" is the transformation matrix from the nodal plane coordinate system to the geographic coordinate system.
+# Variable "rakeVector" is the unit location vector along the slip direction on the nodal plane, expressed in the
+# geographic coordinate system.
+# Variable "n" is the dip vector of the auxilary plane, expressed in radian measure.
+# Variable "s" is the unit vector projection of the rake vector onto the X - Y (horizontal) plane of the
+# geographic coordinate system.
+# Explanation of the result.
+# The first element is the trend of the dip vector of the auxillary plane, in degrees.
+# The second element is the dip angle of the auxillary plane, in degrees.
+# The third element is the trend of the slip vector of the auxillary plane, in degrees.
+
+
+import numpy as np
+
+
+def vector_norm(v):
+    """Computes norm of vector."""
+    norm = np.sqrt(v[0]*v[0] + v[1]*v[1] + v[2]*v[2])
+    return norm
+
+
+def unit_vector(v):
+    """Computes unit vector associated with 3-component vector."""
+    norm = vector_norm(v)
+    unit_vector = [v[0]/norm, v[1]/norm, v[2]/norm]
+    return unit_vector
+
+
+def vector_angle(v1, v2):
+    """Returns angle between two vectors, in radians."""
+    vectorAngle = np.arccos(np.dot(v1, v2)/(vector_norm(v1)*vector_norm(v2)))  # Returns in radians
+    return vectorAngle
+
+
+def find_aux_plane(strike, dip, rake):
+    """
+    Find the second auxiliary plane in a double couple focal mechanism from the given slip plane/vector.
+
+    :param strike: strike of first focal plane, in degrees
+    :type strike: float
+    :param dip: dip of first focal plane, in degrees
+    :type dip: float
+    :param rake: rake of slip vector on first focal plane, in degrees
+    :type rake: float
+    :returns: [strike, dip, rake] of second focal plane
+    :rtype: list of 3 floats
+    """
+    a = strike  # degrees
+    dipaz = strike + 90  # degrees
+    b = np.array([np.sin(np.deg2rad(a)), np.cos(np.deg2rad(a)), 0])
+    c = np.array([np.sin(np.deg2rad(dipaz))*np.cos(np.deg2rad(dip)),
+                  np.cos(np.deg2rad(dipaz))*np.cos(np.deg2rad(dip)),
+                  -np.sin(np.deg2rad(dip))])
+
+    if rake < 0:
+        d = np.cross(b, c)
+    else:
+        d = np.cross(c, b)
+
+    f = np.array([1, 0, 0])
+    g = np.array([0, 1, 0])
+    h = np.array([0, 0, 1])
+
+    if rake < 0:
+        i = b
+        j = c
+        k = np.array([np.cos(np.deg2rad(rake)), -np.sin(np.deg2rad(rake)), 0])
+    else:
+        i = c
+        j = b
+        k = np.array([-np.sin(np.deg2rad(rake)), np.cos(np.deg2rad(rake)), 0])
+
+    m = np.array([[np.dot(f, i), np.dot(f, j), np.dot(f, d)],
+                  [np.dot(g, i), np.dot(g, j), np.dot(g, d)],
+                  [np.dot(h, i), np.dot(h, j), np.dot(h, d)]])
+
+    rakeVector = np.dot(m, k)
+    n = np.arccos(abs(rakeVector[2]))  # auxDipAngle, radians
+
+    s = unit_vector([rakeVector[0], rakeVector[1], 0])
+
+    if rakeVector[0] < 0:
+        rakeTrend = 360-vector_angle(s, g)*(180/np.pi)
+    else:
+        rakeTrend = vector_angle(s, g)*(180/np.pi)
+
+    if rake < 0:
+        if rakeTrend > 270:
+            auxStrike = rakeTrend-270
+        else:
+            auxStrike = rakeTrend+90
+    else:
+        if rakeTrend < 90:
+            auxStrike = rakeTrend+270
+        else:
+            auxStrike = rakeTrend-90
+    auxStrikeRad = auxStrike*(np.pi/180)
+    auxStrikeVector = [np.sin(auxStrikeRad), np.cos(auxStrikeRad), 0]
+
+    if auxStrike > 270:
+        auxDipAz = (auxStrike-270)
+    else:
+        auxDipAz = auxStrike+90
+    auxDipAngle = n*180/np.pi
+
+    if rake < 0:
+        auxRake = -1*vector_angle(auxStrikeVector, d)*180/np.pi
+    else:
+        auxRake = vector_angle(auxStrikeVector, d)*180/np.pi
+
+    auxStrikeAz = auxDipAz-90
+    if auxStrikeAz < 0:
+        auxStrikeAz = auxStrikeAz+360
+
+    return [auxStrikeAz, auxDipAngle, auxRake]

tectonic_utils/seismo/test/test_WC.py

@@ -0,0 +1,19 @@
+import unittest
+from .. import wells_and_coppersmith as wc
+
+
+class Tests(unittest.TestCase):
+
+    def test_simple_focal_plane(self):
+        print("Testing WC")
+        mag = 6.8
+        RLD = wc.RLD_from_M(mag, "N")
+        RW = wc.RW_from_M(mag, "N")
+        expected_RLD = 33.11311214825911  # km
+        expected_RW = 17.378008287493753  # km
+        self.assertEqual(RLD, expected_RLD, 'Wells and Coppersmith wrong')
+        self.assertEqual(RW, expected_RW, 'Wells and Coppersmith wrong')
+
+
+if __name__ == '__main__':
+    unittest.main()

tectonic_utils/seismo/test/test_second_focal_plane.py

@@ -0,0 +1,16 @@
+import unittest
+from .. import second_focal_plane
+
+
+class Tests(unittest.TestCase):
+
+    def test_simple_focal_plane(self):
+        print("Testing second focal plane.")
+        np1 = [38, 51, 44]
+        expected_outcome = [276.71179351444425, 57.32650724684216, 131.61502732561564]
+        np2 = second_focal_plane.find_aux_plane(np1[0], np1[1], np1[2])
+        self.assertEqual(np2, expected_outcome, 'Wrong focal plane')
+
+
+if __name__ == '__main__':
+    unittest.main()

tectonic_utils/seismo/wells_and_coppersmith.py

@@ -0,0 +1,167 @@
+"""
+Set of functions that calculate important scaling relationships between earthquake magnitude and other parameters.
+From Wells and Coppersmith, 1994.
+
+* FAULT_TYPE = SS, R, N, All
+* a, b = parameters in equation
+* sa, sb = parameter uncertainties
+* s = estimate uncertainty
+* r = R correlation coefficient
+
+Average displacements, Max displacements, and inverse relationships not included.
+"""
+
+import numpy as np
+import sys
+from . import moment_calculations
+
+
+def check_fault_types(fault_type):
+    """Ensure that provided fault type is within ["SS", "N", "R", "ALL"]"""
+    if fault_type == "SS" or fault_type == "N" or fault_type == "R" or fault_type == "ALL":
+        return
+    else:
+        sys.exit("Error! Fault Type should be SS, N, R, or ALL!  Exiting...")
+
+
+def rectangular_slip(length, width, magnitude, mu=30*1e9):
+    """Magnitude to slip on rectangular patch
+
+    :param length: fault length, in m
+    :type length: float
+    :param width: fault width, in m
+    :type width: float
+    :param magnitude: magnitude
+    :type magnitude: float
+    :param mu: shear modulus, in Pa, defaults to 30 GPa
+    :type mu: float, optional
+    :returns: slip, in m
+    :rtype: float
+    """
+    area = length * width  # meters^2
+    moment_nm = moment_calculations.moment_from_mw(magnitude)
+    slip = moment_nm / (area * mu)  # Moment = Area * slip * mu.
+    return slip
+
+
+def get_magnitude_from_rectangle(length, width, slip, mu=30*1e9):
+    """
+    :param length: fault length, in m
+    :type length: float
+    :param width: fault width, in m
+    :type width: float
+    :param slip: fault slip, in m
+    :type slip: float
+    :param mu: shear modulus, in Pa, defaults to 30 GPa
+    :type mu: float, optional
+    :returns: magnitude
+    :rtype: float
+    """
+    area = length * width  # meters^2
+    moment_nm = moment_calculations.moment_from_muad(mu, area, slip)  # newton-meters
+    magnitude = moment_calculations.mw_from_moment(moment_nm)
+    return magnitude
+
+
+def SLR_from_M(M, fault_type):
+    """
+    Surface rupture width from Magnitude, computed using log(SLR)=a + b*M
+
+    :param M: magnitude
+    :type M: float
+    :param fault_type: fault_type ["SS", "N", "R", "ALL"]
+    :type fault_type: string
+    :returns: surface rupture length, in km
+    :rtype: float
+    """
+    # row has the form [a, sa, b, sb, s, r, mmin, mmax]
+    check_fault_types(fault_type)
+    if fault_type == "SS":
+        [a, _, b, _, _, _, mmin, mmax] = [-3.55, 0.37, 0.74, 0.05, 0.23, 0.91, 5.6, 8.1]
+    elif fault_type == "R":
+        [a, _, b, _, _, _, mmin, mmax] = [-2.86, 0.55, 0.63, 0.08, 0.20, 0.88, 5.4, 7.4]
+    elif fault_type == "N":
+        [a, _, b, _, _, _, mmin, mmax] = [-2.01, 0.65, 0.50, 0.10, 0.21, 0.81, 5.2, 7.3]
+    else:
+        [a, _, b, _, _, _, mmin, mmax] = [-3.22, 0.27, 0.69, 0.04, 0.22, 0.89, 5.2, 8.1]
+    assert(mmin <= M <= mmax), ValueError("M"+str(M)+" outside valid bounds for Wells and Coppersmith.")
+    SLR = np.power(10, (a + b * M))
+    return SLR
+
+
+def RW_from_M(M, fault_type):
+    """
+    Downdip rupture width from Magnitude, computed using log(RW)=a + b*M
+
+    :param M: magnitude
+    :type M: float
+    :param fault_type: fault_type ["SS", "N", "R", "ALL"]
+    :type fault_type: string
+    :returns: downdip rupture width, in km
+    :rtype: float
+    """
+    # row has the form [a, sa, b, sb, s, r, mmin, mmax]
+    check_fault_types(fault_type)
+    if fault_type == "SS":
+        [a, _, b, _, _, _, mmin, mmax] = [-0.76, 0.12, 0.27, 0.02, 0.14, 0.84, 4.8, 8.1]
+    elif fault_type == "R":
+        [a, _, b, _, _, _, mmin, mmax] = [-1.61, 0.20, 0.41, 0.03, 0.15, 0.90, 4.8, 7.6]
+    elif fault_type == "N":
+        [a, _, b, _, _, _, mmin, mmax] = [-1.14, 0.28, 0.35, 0.05, 0.12, 0.86, 5.2, 7.3]
+    else:
+        [a, _, b, _, _, _, mmin, mmax] = [-1.01, 0.10, 0.32, 0.02, 0.15, 0.84, 4.8, 8.1]
+    assert(mmin <= M <= mmax), ValueError("M"+str(M)+" outside valid bounds for Wells and Coppersmith.")
+    RW = np.power(10, (a + b * M))
+    return RW
+
+
+def RLD_from_M(M, fault_type):
+    """
+    Subsurface rupture length from Magnitude, computed using log(RLD)=a + b*M
+
+    :param M: magnitude
+    :type M: float
+    :param fault_type: fault_type ["SS", "N", "R", "ALL"]
+    :type fault_type: string
+    :returns: subsurface rupture length, in km
+    :rtype: float
+    """
+    # row has the form [a, sa, b, sb, s, r, mmin, mmax]
+    check_fault_types(fault_type)
+    if fault_type == "SS":
+        [a, _, b, _, _, _, mmin, mmax] = [-2.57, 0.12, 0.62, 0.02, 0.15, 0.96, 4.8, 8.1]
+    elif fault_type == "R":
+        [a, _, b, _, _, _, mmin, mmax] = [-2.42, 0.21, 0.58, 0.03, 0.16, 0.93, 4.8, 7.6]
+    elif fault_type == "N":
+        [a, _, b, _, _, _, mmin, mmax] = [-1.88, 0.37, 0.50, 0.06, 0.17, 0.88, 5.2, 7.3]
+    else:
+        [a, _, b, _, _, _, mmin, mmax] = [-2.44, 0.11, 0.59, 0.02, 0.16, 0.94, 4.8, 8.1]
+    assert(mmin <= M <= mmax), ValueError("M"+str(M)+" outside valid bounds for Wells and Coppersmith.")
+    RLD = np.power(10, (a + b * M))
+    return RLD
+
+
+def RA_from_M(M, fault_type):
+    """
+    Rupture Area from Magnitude, computed using log(RA)=a + b*M
+
+    :param M: magnitude
+    :type M: float
+    :param fault_type: fault_type ["SS", "N", "R", "ALL"]
+    :type fault_type: string
+    :returns: rupture area, in km^2
+    :rtype: float
+    """
+    # row has the form [a, sa, b, sb, s, r, mmin, mmax]
+    check_fault_types(fault_type)
+    if fault_type == "SS":
+        [a, _, b, _, _, _, mmin, mmax] = [-3.42, 0.18, 0.90, 0.03, 0.22, 0.96, 4.8, 7.9]
+    elif fault_type == "R":
+        [a, _, b, _, _, _, mmin, mmax] = [-3.99, 0.36, 0.98, 0.06, 0.26, 0.94, 4.8, 7.6]
+    elif fault_type == "N":
+        [a, _, b, _, _, _, mmin, mmax] = [-2.87, 0.50, 0.82, 0.08, 0.22, 0.92, 5.2, 7.3]
+    else:  # ALL
+        [a, _, b, _, _, _, mmin, mmax] = [-3.49, 0.16, 0.91, 0.03, 0.24, 0.95, 4.8, 7.9]
+    assert(mmin <= M <= mmax), ValueError("M"+str(M)+" outside valid bounds for Wells and Coppersmith.")
+    RA = np.power(10, (a + b * M))
+    return RA

tectonic_utils-0.1.2.dist-info/LICENSE.md

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2020 Kathryn Materna
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

tectonic_utils-0.1.2.dist-info/METADATA

@@ -0,0 +1,82 @@
+Metadata-Version: 2.1
+Name: tectonic-utils
+Version: 0.1.2
+Summary: A small package of useful geophysics utilities
+Home-page: https://github.com/kmaterna/Tectonic_Utils
+License: MIT
+Author: Kathryn Materna
+Author-email: kmaterna@berkeley.edu
+Requires-Python: >=3.7
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.7
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Requires-Dist: netCDF4
+Requires-Dist: numpy (>=1.20,<2.0)
+Requires-Dist: scipy
+Project-URL: Repository, https://github.com/kmaterna/Tectonic_Utils
+Description-Content-Type: text/markdown
+
+# Tectonic Utils
+
+A simple library with utility functions for "everyday use" in active tectonics research. 
+There are geodesy and seismology utility functions in Python, 
+some originally from libraries in Matlab or other languages. 
+There are read-write functions to access commonly used formats. 
+Please explore and let me know if you'd like to contribute! 
+
+Functions in this library that I use regularly:
+* Haversine formula
+* Euler Pole rotations
+* Geographic to XYZ Earth-centered coordinates
+* InSAR look vectors (Fialko et al., 2001)
+* Fault vector operations (Aki and Richards, 1980)
+* Moment magnitude conversions (Hanks and Kanamori, 1979)
+* Earthquake magnitude scaling (Wells and Coppersmith, 1994)
+* Reading GMT multi-segment files into Python
+
+![cover](https://github.com/kmaterna/Tectonic_Utils/blob/master/Tectonic_Utils/cover_picture.png)
+
+## Dependencies
+* Python 3
+* numpy
+* scipy 
+* netCDF4 (pip install netCDF4)
+
+
+## Installation
+On Unix (other systems not tested yet), install from PyPI using: 
+
+```pip install Tectonic-Utils```
+
+To test that it works, import a utility in Python. I use an import structure like: 
+```
+$ python
+>>> from Tectonic_Utils.geodesy import haversine
+```
+
+## Documentation
+The preliminary sphinx html docs can be found here: https://kmaterna.github.io/Tectonic_Utils/index.html
+
+## Additional Thoughts
+I wrote this library for my own research because (**in theory**) 
+putting all my Python utility functions 
+into one unit-tested repository reduces bugs, increases code reproducibility, and allows for faster experimentation. 
+The library is open source and available on Pip in case it would help others too. 
+
+Disclaimer: Completely bug-free code can never be guaranteed. As things come up, I occasionally fix and update the package on Pip/PyPI.    
+
+If you're using this package and you'd like to contribute or you find bugs, please let me know!  
+
+## References
+
+* Aki, K., and P. G. Richards (1980), Quantitative Seismology: Theory and Methods, W. H. Freeman, New York.
+* Fialko, Y., Simons, M., & Agnew, D. (2001). The complete (3-D) surface displacement field in the epicentral area of the 1999 Mw7.1 Hector Mine earthquake, California, from space geodetic observations. Geophysical Research Letters, 28(16), 3063–3066.
+* Hanks, T., & Kanamori, H. (1979). A moment magnitude scale. Journal of Geophysical Research, 84(B5), 2348–2350.
+* Wells, D. L., & Coppersmith, K. J. (1994). New Empirical Relationships among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement. Bulletin of the Seismological Society of America, 84(4), 974–1002. 
+

tectonic_utils-0.1.2.dist-info/RECORD

@@ -0,0 +1,51 @@
+tectonic_utils/.DS_Store,sha256=kRqKlYZrl1RGaoPXapBlMI2wNLE1cqcvn4IknBT8D2Q,6148
+tectonic_utils/__init__.py,sha256=KhFOyyd1QUbOdza1-ovPLepPkHPoNBYLX0axjvz2Aw4,68
+tectonic_utils/cover_picture.png,sha256=SioKXhjLRmjUBbVKY-4saEEO3yUVDB06P50Hv8XpLJ8,1733968
+tectonic_utils/geodesy/.DS_Store,sha256=oELC_Tpc4DH6jDWVlImUeFSN8EDmW2TsXrWR0300T5Y,6148
+tectonic_utils/geodesy/.ruff_cache/.gitignore,sha256=aEiIwOuxfzdCmLZe4oB1JsBmCUxwG8x-u-HBCV9JT8E,1
+tectonic_utils/geodesy/.ruff_cache/0.1.5/15663111236935520357,sha256=iRigcC_uE9RzPVJMM_H2-2XRweGq7TPZmatXT-J4Yhk,1014
+tectonic_utils/geodesy/.ruff_cache/CACHEDIR.TAG,sha256=WVMVbX4MVkpCclExbq8m-IcOZIOuIZf5FrYw5Pk-Ma4,43
+tectonic_utils/geodesy/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/geodesy/datums.py,sha256=z_IAhNFMSz7avDt2WhxZnU5Y2SMjMrVdszO0BwYEbsM,11000
+tectonic_utils/geodesy/euler_pole.py,sha256=0KNBnmHegL0ARA6uSySEZBtmzC3hfMUYjBoNrHAjGpk,6418
+tectonic_utils/geodesy/fault_vector_functions.py,sha256=vAdIX2evtitrheyK4LWNIEx3BNskp3J0jvyYDgl40rM,12823
+tectonic_utils/geodesy/haversine.py,sha256=aRoeztLfe8zFiNKupWN33AA9VrS5uU-ZAzQVQf5Oqao,6945
+tectonic_utils/geodesy/insar_vector_functions.py,sha256=1MLy5J2W3jHjMbtyeuLL74HPLF8CW1ztIxHFaLXsWhc,11446
+tectonic_utils/geodesy/linear_elastic.py,sha256=wdTB4m_nDbUCUPNlUxjTScfAjCxGIUqhN4wRhYNeqng,6476
+tectonic_utils/geodesy/test/.DS_Store,sha256=1lFlJ5EFymdzGAUAaI30vcaaLHt3F1LwpG7xILf9jsM,6148
+tectonic_utils/geodesy/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/geodesy/test/test_conversion_functions.py,sha256=TVrs3LCMQWPxRUMmvzKAB1f3-Oc40jZk9F4E7uZq28U,3709
+tectonic_utils/geodesy/test/test_euler_poles.py,sha256=P4rNRnalbfnkCuXnyiuzHjMoqguA1k-2pMEyfch9hvE,1375
+tectonic_utils/geodesy/test/test_insar_vector_functions.py,sha256=ChHB3dwC-cHpFuSBCxOII1KcoPVNVSIRCEY2ifWsK1o,1561
+tectonic_utils/geodesy/utilities.py,sha256=538vb-4BCYVUWkwl9V5tKS9DnjJ3KQFjRMPtrKiDjbo,1163
+tectonic_utils/geodesy/xyz2llh.py,sha256=oPDS7nWl7gJe6WD7DQqvKlrffV69avdpiy0epUwcjzE,8083
+tectonic_utils/read_write/.DS_Store,sha256=SkDIyIgu4JN5DmqhlYPZoyLgbY4WmKMc6hw6usKijGk,6148
+tectonic_utils/read_write/.ruff_cache/.gitignore,sha256=aEiIwOuxfzdCmLZe4oB1JsBmCUxwG8x-u-HBCV9JT8E,1
+tectonic_utils/read_write/.ruff_cache/0.1.5/680373307893520726,sha256=M_jn9qQ1oFwhM09Msm3b_JM7m-TQnsoZ-hR3a17VfVU,177
+tectonic_utils/read_write/.ruff_cache/CACHEDIR.TAG,sha256=WVMVbX4MVkpCclExbq8m-IcOZIOuIZf5FrYw5Pk-Ma4,43
+tectonic_utils/read_write/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/read_write/general_io.py,sha256=HfXgGqePenAqLQ0WQJEmgxygKORy_9DR5TLlSak9Ezs,2036
+tectonic_utils/read_write/netcdf_read_write.py,sha256=hFlyDNJFOPaQwihzi0heyW6XNUyLqTQK3MUKwx_MS-Q,14753
+tectonic_utils/read_write/read_kml.py,sha256=G1wzrCPRljnvRja1dV6-Na6BHqjRLciST2Oh1jgqHK0,2300
+tectonic_utils/read_write/test/.DS_Store,sha256=1lFlJ5EFymdzGAUAaI30vcaaLHt3F1LwpG7xILf9jsM,6148
+tectonic_utils/read_write/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/read_write/test/example_grd.grd,sha256=tJD_Q6IqDMdi2SuN9f0Z70yKyLpUumgDtIqUM42G-48,60395
+tectonic_utils/read_write/test/test_conversion_functions.py,sha256=4Li_lFOiXDFAQ0AvN7OsfuRbJS2V_jd_eNqyKh0B7ck,1512
+tectonic_utils/read_write/test/written_example.grd,sha256=9Ium7MjpxFOvNnr00iYb78jN0n50aiGF26ynnVMrC9M,58986
+tectonic_utils/seismo/.DS_Store,sha256=l3pS6MjdhXRwqDQMYSmLH9GgZ90rVhOAgOySNGfW8Nk,6148
+tectonic_utils/seismo/.ruff_cache/.gitignore,sha256=aEiIwOuxfzdCmLZe4oB1JsBmCUxwG8x-u-HBCV9JT8E,1
+tectonic_utils/seismo/.ruff_cache/0.1.5/12911000862714636977,sha256=y0tMlRMFpMQ44MESW-KBEqvl9vPYdp2UEvfBKMXvyW4,173
+tectonic_utils/seismo/.ruff_cache/CACHEDIR.TAG,sha256=WVMVbX4MVkpCclExbq8m-IcOZIOuIZf5FrYw5Pk-Ma4,43
+tectonic_utils/seismo/MT_calculations.py,sha256=ccd-DOVr4Vhj21WfWeicetmICOVOrST8rS0ow_GSt5w,4146
+tectonic_utils/seismo/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/seismo/moment_calculations.py,sha256=ihrBYwhvAPhfT-QZ6oM-lSktjZrnto6TC5N9plSQ9po,1146
+tectonic_utils/seismo/second_focal_plane.py,sha256=8GTSkeBg6owHk5-KsE6LM0b1NhoPjtQyLXwLu1sDCJo,5436
+tectonic_utils/seismo/test/.DS_Store,sha256=1lFlJ5EFymdzGAUAaI30vcaaLHt3F1LwpG7xILf9jsM,6148
+tectonic_utils/seismo/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tectonic_utils/seismo/test/test_WC.py,sha256=d7qNTPAPDKPu2bRlpS5f5pNwRfrN7ZbxKhBGuIu-ZuQ,544
+tectonic_utils/seismo/test/test_second_focal_plane.py,sha256=bb5-wMkeLgZPXW5qpqXo4BMNIjKQdVP8ZlAnI_uw3QA,473
+tectonic_utils/seismo/wells_and_coppersmith.py,sha256=NXbUVXDgifqTro-DeOXdFe6_eOQ82sVg6_rWvg0p-g0,6003
+tectonic_utils-0.1.2.dist-info/LICENSE.md,sha256=BoWwOxP7snl4zEufDQFfftDiyZgYqGr09E7brqIs8M4,1072
+tectonic_utils-0.1.2.dist-info/METADATA,sha256=j-vlc-aRZKbDMMTOfefURKVn4gUPTAINe52eDRLFTaQ,3524
+tectonic_utils-0.1.2.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
+tectonic_utils-0.1.2.dist-info/RECORD,,

tectonic_utils-0.1.2.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: poetry-core 1.8.1
+Root-Is-Purelib: true
+Tag: py3-none-any