Phase 3–8: remove SourceTheoryMotion/StaffordEtAl22Motion, add pygmm interop contracts

Deletes SourceTheoryMotion, StaffordEtAl22Motion, calc_stress_drop,
calc_geometric_spreading from motions.py and the sea22-site_amp data file
(now in pygmm.fourier_spectrum); adds _contracts.py with RvtMotion helpers;
adds pygmm to dev requirements; updates CHANGELOG for 2.0.0.

Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>

Author: arkottke

CHANGELOG.md

@@ -1,5 +1,25 @@
 # Changelog
 
+## 2.0.0 (unreleased)
+
+**Breaking changes**:
+
+- Removed `SourceTheoryMotion` and `StaffordEtAl22Motion` — use
+  `pygmm.fourier_spectrum.SourceTheoryModel` / `StaffordEtAl2022` instead.
+  Pipe the result through `RvtMotion.from_fas(model)` to reproduce previous behaviour.
+- Removed `calc_stress_drop` and `calc_geometric_spreading` module-level helpers.
+- Removed `event_kwds` parameter from `CompatibleRvtMotion.__init__`; `duration` is
+  now required. Use `CompatibleRvtMotion.from_response_spectrum(rs, duration)` or pass
+  `duration` explicitly.
+
+Added:
+
+- `RvtMotion.from_fas(fas, peak_calculator=None)` — create a motion from any object
+  with `.freqs`, `.fourier_amps`, `.duration` attributes.
+- `CompatibleRvtMotion.from_response_spectrum(rs, duration, **kw)` — create a
+  compatible motion from any object with `.periods`, `.spec_accels`, `.damping`.
+- `RvtMotion.calc_pga()` and `.calc_pgv()`.
+
 ## 0.8.0 (dev)
 
 - Added Stafford et al. (2022) ground motion model.

examples/run_example_1.py

@@ -14,8 +14,7 @@
 damping = 0.05
 method = "BJ84"
 
-event_keys = ["magnitude", "distance", "region"]
-event_kwds = {key: event[key] for key in event_keys}
+calc_kwds = {key: event[key] for key in ["magnitude", "distance", "region"]}
 
 for method in ["BJ84", "CLH56", "DK85", "LP99", "TM87", "V75"]:
     start = datetime.datetime.now()
@@ -24,8 +23,7 @@
         event["psa"],
         duration=event["duration"],
         osc_damping=damping,
-        event_kwds=event_kwds,
-        peak_calculator=pyrvt.tools.get_peak_calculator(method, event_kwds),
+        peak_calculator=pyrvt.tools.get_peak_calculator(method, calc_kwds),
     )
     psa_calc = crm.calc_osc_accels(target_freqs, damping)
     duration = datetime.datetime.now() - start

requirements_dev.txt

@@ -9,6 +9,7 @@ pre-commit
 pydocstyle
 pyexcel
 pyexcel-xlsx
+pygmm
 pystrata
 pytest>=3.5
 pytest-cov

src/pyrvt/_contracts.py

@@ -0,0 +1,38 @@
+"""Consumer-side data contracts for pyrvt ↔ pygmm interop.
+
+Canonical definitions live in ``pygmm.contracts``; this file is a private
+~35-LOC duplicate so pyrvt has no runtime dependency on pygmm.
+"""
+
+from __future__ import annotations
+
+from dataclasses import dataclass
+from typing import Optional
+
+import numpy as np
+import numpy.typing as npt
+
+
+@dataclass(frozen=True)
+class FourierSpectrum:
+    """One-sided Fourier amplitude spectrum of acceleration.
+
+    Field names match ``pygmm.contracts.FourierSpectrum``.
+    """
+
+    freqs: npt.NDArray[np.floating]
+    fourier_amps: npt.NDArray[np.floating]
+    duration: float
+
+
+@dataclass(frozen=True)
+class ResponseSpectrum:
+    """Pseudo-spectral acceleration response spectrum.
+
+    Field names match ``pygmm.contracts.ResponseSpectrum``.
+    """
+
+    periods: npt.NDArray[np.floating]
+    spec_accels: npt.NDArray[np.floating]
+    damping: float
+    duration: Optional[float] = None

src/pyrvt/_version.py

@@ -0,0 +1,24 @@
+# file generated by vcs-versioning
+# don't change, don't track in version control
+from __future__ import annotations
+
+__all__ = [
+    "__version__",
+    "__version_tuple__",
+    "version",
+    "version_tuple",
+    "__commit_id__",
+    "commit_id",
+]
+
+version: str
+__version__: str
+__version_tuple__: tuple[int | str, ...]
+version_tuple: tuple[int | str, ...]
+commit_id: str | None
+__commit_id__: str | None
+
+__version__ = version = '0.8.2.dev20+gdcab64e8f'
+__version_tuple__ = version_tuple = (0, 8, 2, 'dev20', 'gdcab64e8f')
+
+__commit_id__ = commit_id = 'gdcab64e8f'