weac 3.0.1__py3-none-any.whl → 3.0.2__py3-none-any.whl

weac/__init__.py

@@ -2,4 +2,4 @@
 WEAC - Weak Layer Anticrack Nucleation Model
 """
 
-__version__ = "3.0.1"
+__version__ = "3.0.2"

weac/analysis/criteria_evaluator.py

@@ -9,7 +9,6 @@ import logging
 import time
 import warnings
 from dataclasses import dataclass
-from typing import List, Optional, Union
 
 # Third party imports
 import numpy as np
@@ -34,12 +33,12 @@ logger = logging.getLogger(__name__)
 class CoupledCriterionHistory:
     """Stores the history of the coupled criterion evaluation."""
 
-    skier_weights: List[float]
-    crack_lengths: List[float]
-    incr_energies: List[np.ndarray]
-    g_deltas: List[float]
-    dist_maxs: List[float]
-    dist_mins: List[float]
+    skier_weights: list[float]
+    crack_lengths: list[float]
+    incr_energies: list[np.ndarray]
+    g_deltas: list[float]
+    dist_maxs: list[float]
+    dist_mins: list[float]
 
 
 @dataclass
@@ -89,7 +88,7 @@ class CoupledCriterionResult:
     g_delta: float
     dist_ERR_envelope: float
     iterations: int
-    history: Optional[CoupledCriterionHistory]
+    history: CoupledCriterionHistory | None
     final_system: SystemModel
     max_dist_stress: float
     min_dist_stress: float
@@ -130,9 +129,9 @@ class FindMinimumForceResult:
         Whether the algorithm converged.
     critical_skier_weight : float
         The critical skier weight.
-    new_segments : List[Segment]
+    new_segments : list[Segment]
         The new segments.
-    old_segments : List[Segment]
+    old_segments : list[Segment]
         The old segments.
     iterations : int
         The number of iterations.
@@ -144,9 +143,9 @@ class FindMinimumForceResult:
 
     success: bool
     critical_skier_weight: float
-    new_segments: List[Segment]
-    old_segments: List[Segment]
-    iterations: Optional[int]
+    new_segments: list[Segment]
+    old_segments: list[Segment]
+    iterations: int | None
     max_dist_stress: float
     min_dist_stress: float
 
@@ -203,10 +202,10 @@ class CriteriaEvaluator:
 
     def stress_envelope(
         self,
-        sigma: Union[float, np.ndarray],
-        tau: Union[float, np.ndarray],
+        sigma: float | np.ndarray,
+        tau: float | np.ndarray,
         weak_layer: WeakLayer,
-        method: Optional[str] = None,
+        method: str | None = None,
     ) -> np.ndarray:
         """
         Evaluate the stress envelope for given stress components.
@@ -672,16 +671,20 @@ class CriteriaEvaluator:
                 UserWarning,
             )
         system_copy = copy.deepcopy(system)
+        system_copy.config.touchdown = True
+        system_copy.update_scenario(scenario_config=ScenarioConfig(phi=0.0))
+        l_BC = system.slab_touchdown.l_BC
+
         segments = [
             Segment(length=5e3, has_foundation=True, m=0.0),
-            Segment(length=5e3, has_foundation=False, m=0.0),
+            Segment(length=2 * l_BC, has_foundation=False, m=0.0),
         ]
         scenario_config = ScenarioConfig(
             system_type="vpst-" if vertical else "pst-",
-            phi=system.scenario.phi,
-            cut_length=5e3,
+            phi=0.0,  # Slab Touchdown works only for flat slab
+            cut_length=2 * l_BC,
         )
-        system_copy.config.touchdown = True
+        # system_copy.config.touchdown = True
         system_copy.update_scenario(segments=segments, scenario_config=scenario_config)
         touchdown_distance = system_copy.slab_touchdown.touchdown_distance
         analyzer = Analyzer(system_copy, printing_enabled=print_call_stats)
@@ -825,7 +828,7 @@ class CriteriaEvaluator:
         system: SystemModel,
         search_interval: tuple[float, float] | None = None,
         target: float = 1,
-    ) -> tuple[float, List[Segment]]:
+    ) -> tuple[float, list[Segment]]:
         """
         Finds the minimum crack length required to surpass the energy release rate envelope.
 
@@ -838,7 +841,7 @@ class CriteriaEvaluator:
         --------
         minimum_crack_length: float
             The minimum crack length required to surpass the energy release rate envelope [mm]
-        new_segments: List[Segment]
+        new_segments: list[Segment]
             The updated list of segments
         """
         old_segments = copy.deepcopy(system.scenario.segments)
@@ -928,7 +931,7 @@ class CriteriaEvaluator:
         self,
         system: SystemModel,
         skier_weight: float,
-    ) -> tuple[float, List[Segment]]:
+    ) -> tuple[float, list[Segment]]:
         """
         Finds the resulting anticrack length and updated segment configurations
         for a given skier weight.
@@ -944,7 +947,7 @@ class CriteriaEvaluator:
         -------
         new_crack_length: float
             The total length of the new cracked segments [mm]
-        new_segments: List[Segment]
+        new_segments: list[Segment]
             The updated list of segments
         """
         logger.info(
@@ -1086,7 +1089,7 @@ class CriteriaEvaluator:
 
     def _find_stress_envelope_crossings(
         self, system: SystemModel, weak_layer: WeakLayer
-    ) -> List[float]:
+    ) -> list[float]:
         """
         Finds the exact x-coordinates where the stress envelope is crossed.
         """

weac/analysis/plotter.py

@@ -6,7 +6,7 @@ This module provides plotting functions for visualizing the results of the WEAC
 import colorsys
 import logging
 import os
-from typing import List, Literal, Optional
+from typing import Literal
 
 # Third party imports
 import matplotlib.colors as mc
@@ -219,9 +219,9 @@ class Plotter:
 
     def _get_systems_to_plot(
         self,
-        system_model: Optional[SystemModel] = None,
-        system_models: Optional[List[SystemModel]] = None,
-    ) -> List[SystemModel]:
+        system_model: SystemModel | None = None,
+        system_models: list[SystemModel] | None = None,
+    ) -> list[SystemModel]:
         """Determine which systems to plot based on override parameters."""
         if system_model is not None and system_models is not None:
             raise ValueError(
@@ -236,7 +236,7 @@ class Plotter:
             "SystemModel or list of SystemModels"
         )
 
-    def _save_figure(self, filename: str, fig: Optional[Figure] = None):
+    def _save_figure(self, filename: str, fig: Figure | None = None):
         """Save figure with proper formatting."""
         if fig is None:
             fig = plt.gcf()
@@ -249,20 +249,20 @@ class Plotter:
 
     def plot_slab_profile(
         self,
-        weak_layers: List[WeakLayer] | WeakLayer,
-        slabs: List[Slab] | Slab,
+        weak_layers: list[WeakLayer] | WeakLayer,
+        slabs: list[Slab] | Slab,
         filename: str = "slab_profile",
-        labels: Optional[List[str] | str] = None,
-        colors: Optional[List[str]] = None,
+        labels: list[str] | str | None = None,
+        colors: list[str] | None = None,
     ):
         """
         Plot slab layer profiles for comparison.
 
         Parameters
         ----------
-        weak_layers : List[WeakLayer] | WeakLayer
+        weak_layers : list[WeakLayer] | WeakLayer
             The weak layer or layers to plot.
-        slabs : List[Slab] | Slab
+        slabs : list[Slab] | Slab
             The slab or slabs to plot.
         filename : str, optional
             Filename for saving plot
@@ -281,6 +281,9 @@ class Plotter:
         if isinstance(slabs, Slab):
             slabs = [slabs]
 
+        if len(weak_layers) != len(slabs):
+            raise ValueError("Number of weak layers must match number of slabs")
+
         if labels is None:
             labels = [f"System {i + 1}" for i in range(len(weak_layers))]
         elif isinstance(labels, str):
@@ -642,11 +645,11 @@ class Plotter:
 
     def plot_section_forces(
         self,
-        system_model: Optional[SystemModel] = None,
-        system_models: Optional[List[SystemModel]] = None,
+        system_model: SystemModel | None = None,
+        system_models: list[SystemModel] | None = None,
         filename: str = "section_forces",
-        labels: Optional[List[str]] = None,
-        colors: Optional[List[str]] = None,
+        labels: list[str] | None = None,
+        colors: list[str] | None = None,
     ):
         """
         Plot section forces (N, M, V) for comparison.
@@ -655,7 +658,7 @@ class Plotter:
         ----------
         system_model : SystemModel, optional
             Single system to plot (overrides default)
-        system_models : List[SystemModel], optional
+        system_models : list[SystemModel], optional
             Multiple systems to plot (overrides default)
         filename : str, optional
             Filename for saving plot
@@ -721,11 +724,11 @@ class Plotter:
 
     def plot_energy_release_rates(
         self,
-        system_model: Optional[SystemModel] = None,
-        system_models: Optional[List[SystemModel]] = None,
+        system_model: SystemModel | None = None,
+        system_models: list[SystemModel] | None = None,
         filename: str = "ERR",
-        labels: Optional[List[str]] = None,
-        colors: Optional[List[str]] = None,
+        labels: list[str] | None = None,
+        colors: list[str] | None = None,
     ):
         """
         Plot energy release rates (G_I, G_II) for comparison.
@@ -734,7 +737,7 @@ class Plotter:
         ----------
         system_model : SystemModel, optional
             Single system to plot (overrides default)
-        system_models : List[SystemModel], optional
+        system_models : list[SystemModel], optional
             Multiple systems to plot (overrides default)
         filename : str, optional
             Filename for saving plot
@@ -1035,7 +1038,7 @@ class Plotter:
         system_model: SystemModel,
         criteria_evaluator: CriteriaEvaluator,
         all_envelopes: bool = False,
-        filename: Optional[str] = None,
+        filename: str | None = None,
     ):
         """
         Plot stress envelope in τ-σ space.
@@ -1074,7 +1077,7 @@ class Plotter:
         weak_layer = system_model.weak_layer
 
         # Define a function to find the root for a given tau
-        def find_sigma_for_tau(tau_val, sigma_c, method: Optional[str] = None):
+        def find_sigma_for_tau(tau_val, sigma_c, method: str | None = None):
             # Target function to find the root of: envelope(sigma, tau) - 1 = 0
             def envelope_root_func(sigma_val):
                 return (

weac/components/model_input.py

@@ -14,7 +14,6 @@ field_name: type = Field(..., gt=0, description="Description")
 
 import json
 import logging
-from typing import List
 
 from pydantic import BaseModel, ConfigDict, Field, model_validator
 
@@ -49,13 +48,13 @@ class ModelInput(BaseModel):
         default_factory=lambda: WeakLayer(rho=125, h=20, E=1.0),
         description="Weak layer",
     )
-    layers: List[Layer] = Field(
-        default_factory=lambda: [Layer(rho=250, h=100)], description="List of layers"
+    layers: list[Layer] = Field(
+        default_factory=lambda: [Layer(rho=250, h=100)], description="list of layers"
     )
     scenario_config: ScenarioConfig = Field(
         default_factory=ScenarioConfig, description="Scenario configuration"
     )
-    segments: List[Segment] = Field(
+    segments: list[Segment] = Field(
         default_factory=lambda: [
             Segment(length=5000, has_foundation=True, m=100),
             Segment(length=5000, has_foundation=True, m=0),

weac/core/eigensystem.py

@@ -4,10 +4,8 @@ the eigenvalue problem for a layered beam on an elastic foundation.
 """
 
 import logging
-from typing import Optional
 
 import numpy as np
-from numpy.typing import NDArray
 
 from weac.components import WeakLayer
 from weac.constants import SHEAR_CORRECTION_FACTOR
@@ -58,17 +56,17 @@ class Eigensystem:
     kA55: float  # shear stiffness
     K0: float  # foundation stiffness
 
-    K: NDArray  # System Matrix
+    K: np.ndarray  # System Matrix
 
     # Eigenvalues and Eigenvectors
-    ewC: NDArray[np.complex128]  # shape (k): Complex Eigenvalues
-    ewR: NDArray[np.float64]  # shape (k): Real Eigenvalues
-    evC: NDArray[np.complex128]  # shape (6, k): Complex Eigenvectors
-    evR: NDArray[np.float64]  # shape (6, k): Real Eigenvectors
-    sR: NDArray[
+    ewC: np.ndarray[np.complex128]  # shape (k): Complex Eigenvalues
+    ewR: np.ndarray[np.float64]  # shape (k): Real Eigenvalues
+    evC: np.ndarray[np.complex128]  # shape (6, k): Complex Eigenvectors
+    evR: np.ndarray[np.float64]  # shape (6, k): Real Eigenvectors
+    sR: np.ndarray[
         np.float64
     ]  # shape (k): Real positive eigenvalue shifts (for numerical robustness)
-    sC: NDArray[
+    sC: np.ndarray[
         np.float64
     ]  # shape (k): Complex positive eigenvalue shifts (for numerical robustness)
 
@@ -114,8 +112,8 @@ class Eigensystem:
         self.K0 = B11**2 - A11 * D11
 
     def assemble_system_matrix(
-        self, kn: Optional[float], kt: Optional[float]
-    ) -> NDArray[np.float64]:
+        self, kn: float | None, kt: float | None
+    ) -> np.ndarray[np.float64]:
         """
         Assemble first-order ODE system matrix K.
 
@@ -171,14 +169,14 @@ class Eigensystem:
         return np.array(K, dtype=np.float64)
 
     def calc_eigenvalues_and_eigenvectors(
-        self, system_matrix: NDArray[np.float64]
+        self, system_matrix: np.ndarray[np.float64]
     ) -> tuple[
-        NDArray[np.complex128],
-        NDArray[np.float64],
-        NDArray[np.complex128],
-        NDArray[np.float64],
-        NDArray[np.float64],
-        NDArray[np.float64],
+        np.ndarray[np.complex128],
+        np.ndarray[np.float64],
+        np.ndarray[np.complex128],
+        np.ndarray[np.float64],
+        np.ndarray[np.float64],
+        np.ndarray[np.float64],
     ]:
         """
         Calculate eigenvalues and eigenvectors of the system matrix.
@@ -215,7 +213,9 @@ class Eigensystem:
         sR[ewR > 0], sC[ewC > 0] = -1, -1
         return ewC, ewR, evC, evR, sR, sC
 
-    def zh(self, x: float, length: float = 0, has_foundation: bool = True) -> NDArray:
+    def zh(
+        self, x: float, length: float = 0, has_foundation: bool = True
+    ) -> np.ndarray:
         """
         Compute bedded or free complementary solution at position x.
 
@@ -275,7 +275,7 @@ class Eigensystem:
 
     def zp(
         self, x: float, phi: float = 0, has_foundation=True, qs: float = 0
-    ) -> NDArray:
+    ) -> np.ndarray:
         """
         Compute bedded or free particular integrals at position x.
 
@@ -357,7 +357,7 @@ class Eigensystem:
 
         return zp
 
-    def get_load_vector(self, phi: float, qs: float = 0) -> NDArray:
+    def get_load_vector(self, phi: float, qs: float = 0) -> np.ndarray:
         """
         Compute system load vector q.
 

weac/core/scenario.py

@@ -3,7 +3,7 @@ This module defines the Scenario class, which encapsulates the physical setup of
 """
 
 import logging
-from typing import List, Sequence, Union
+from collections.abc import Sequence
 
 import numpy as np
 
@@ -45,7 +45,7 @@ class Scenario:
 
     # Inputs
     scenario_config: ScenarioConfig
-    segments: List[Segment]
+    segments: list[Segment]
     weak_layer: WeakLayer
     slab: Slab
 
@@ -69,7 +69,7 @@ class Scenario:
     def __init__(
         self,
         scenario_config: ScenarioConfig,
-        segments: List[Segment],
+        segments: list[Segment],
         weak_layer: WeakLayer,
         slab: Slab,
     ):
@@ -102,19 +102,19 @@ class Scenario:
         self._calc_crack_height()
 
     def get_segment_idx(
-        self, x: Union[float, Sequence[float], np.ndarray]
-    ) -> Union[int, np.ndarray]:
+        self, x: float | Sequence[float] | np.ndarray
+    ) -> int | np.ndarray:
         """
         Get the segment index for a given x-coordinate or coordinates.
 
         Parameters
         ----------
-        x: Union[float, Sequence[float], np.ndarray]
+        x: float | Sequence[float] | np.ndarray
             A single x-coordinate or a sequence of x-coordinates.
 
         Returns
         -------
-        Union[int, np.ndarray]
+        int | np.ndarray
             The segment index or an array of indices.
         """
         x_arr = np.asarray(x)

weac/core/slab.py

@@ -2,8 +2,6 @@
 This module defines the Slab class, which represents the snow slab and its properties.
 """
 
-from typing import List
-
 import numpy as np
 
 from weac.components import Layer
@@ -44,7 +42,7 @@ class Slab:  # pylint: disable=too-many-instance-attributes,too-few-public-metho
     """
 
     # Input data
-    layers: List[Layer]
+    layers: list[Layer]
 
     rhoi: np.ndarray  # densities of the layer i [t/mm^3]
     hi: np.ndarray  # thickness of the layer i [mm]
@@ -60,7 +58,7 @@ class Slab:  # pylint: disable=too-many-instance-attributes,too-few-public-metho
     z_cog: float  # z-coordinate of Center of Gravity [mm]
     qw: float  # Weight Load of the slab [N/mm]
 
-    def __init__(self, layers: List[Layer]) -> None:
+    def __init__(self, layers: list[Layer]) -> None:
         self.layers = layers
         self._calc_slab_params()
 

weac/core/slab_touchdown.py

@@ -4,7 +4,7 @@ Handling the touchdown situation in a PST.
 """
 
 import logging
-from typing import Literal, Optional
+from typing import Literal
 
 from scipy.optimize import brentq
 
@@ -60,7 +60,7 @@ class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-pub
         Type of touchdown mode
     touchdown_distance : float
         Length of the touchdown segment [mm]
-    collapsed_weak_layer_kR : Optional[float]
+    collapsed_weak_layer_kR : float | None
         Rotational spring stiffness of the collapsed weak layer segment
     """
 
@@ -78,7 +78,7 @@ class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-pub
         "A_free_hanging", "B_point_contact", "C_in_contact"
     ]  # Three types of contact with collapsed weak layer
     touchdown_distance: float
-    collapsed_weak_layer_kR: Optional[float] = None
+    collapsed_weak_layer_kR: float | None = None
 
     def __init__(self, scenario: Scenario, eigensystem: Eigensystem):
         self.scenario = scenario

weac/core/system_model.py

@@ -12,7 +12,6 @@ import copy
 import logging
 from collections.abc import Sequence
 from functools import cached_property
-from typing import List, Optional, Union
 
 import numpy as np
 
@@ -119,11 +118,11 @@ class SystemModel:
     fq: FieldQuantities
 
     scenario: Scenario
-    slab_touchdown: Optional[SlabTouchdown]
+    slab_touchdown: SlabTouchdown | None
     unknown_constants: np.ndarray
     uncracked_unknown_constants: np.ndarray
 
-    def __init__(self, model_input: ModelInput, config: Optional[Config] = None):
+    def __init__(self, model_input: ModelInput, config: Config | None = None):
         if config is None:
             config = Config()
         self.config = config
@@ -155,7 +154,7 @@ class SystemModel:
         return Eigensystem(weak_layer=self.weak_layer, slab=self.slab)
 
     @cached_property
-    def slab_touchdown(self) -> Optional[SlabTouchdown]:
+    def slab_touchdown(self) -> SlabTouchdown | None:
         """
         Solve for the slab touchdown.
         Modifies the scenario object in place by replacing the undercut segment
@@ -303,7 +302,7 @@ class SystemModel:
         self._invalidate_eigensystem()
 
     # Changes that affect the *slab*  -> rebuild everything
-    def update_layers(self, new_layers: List[Layer]):
+    def update_layers(self, new_layers: list[Layer]):
         """Update the layers."""
         slab = Slab(layers=new_layers)
         self.slab = slab
@@ -318,8 +317,8 @@ class SystemModel:
     # Changes that affect the *scenario*  -> only rebuild C constants
     def update_scenario(
         self,
-        segments: Optional[List[Segment]] = None,
-        scenario_config: Optional[ScenarioConfig] = None,
+        segments: list[Segment] | None = None,
+        scenario_config: ScenarioConfig | None = None,
     ):
         """
         Update fields on `scenario_config` (if present) or on the
@@ -365,7 +364,7 @@ class SystemModel:
 
     def z(
         self,
-        x: Union[float, Sequence[float], np.ndarray],
+        x: float | Sequence[float] | np.ndarray,
         C: np.ndarray,
         length: float,
         phi: float,
@@ -396,7 +395,7 @@ class SystemModel:
         z : ndarray
             Solution vector (6xN) at position x.
         """
-        if isinstance(x, (list, tuple, np.ndarray)):
+        if isinstance(x, (np.ndarray, Sequence)):
             z = np.concatenate(
                 [
                     np.dot(self.eigensystem.zh(xi, length, has_foundation), C)

weac/core/unknown_constants_solver.py

@@ -8,7 +8,7 @@ We utilize the pydantic library to define the system model.
 """
 
 import logging
-from typing import Literal, Optional
+from typing import Literal
 
 import numpy as np
 from numpy.linalg import LinAlgError
@@ -36,11 +36,10 @@ class UnknownConstantsSolver:
         scenario: Scenario,
         eigensystem: Eigensystem,
         system_type: SystemType,
-        touchdown_distance: Optional[float] = None,
-        touchdown_mode: Optional[
-            Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
-        ] = None,
-        collapsed_weak_layer_kR: Optional[float] = None,
+        touchdown_distance: float | None = None,
+        touchdown_mode: Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
+        | None = None,
+        collapsed_weak_layer_kR: float | None = None,
     ) -> np.ndarray:
         """
         Compute free constants *C* for system. \\
@@ -264,10 +263,9 @@ class UnknownConstantsSolver:
         has_foundation: bool,
         pos: Literal["l", "r", "m", "left", "right", "mid"],
         system_type: SystemType,
-        touchdown_mode: Optional[
-            Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
-        ] = None,
-        collapsed_weak_layer_kR: Optional[float] = None,
+        touchdown_mode: Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
+        | None = None,
+        collapsed_weak_layer_kR: float | None = None,
     ) -> np.ndarray:
         """
         Provide boundary or transmission conditions for beam segments.
@@ -369,10 +367,9 @@ class UnknownConstantsSolver:
         has_foundation: bool,
         pos: Literal["l", "r", "m", "left", "right", "mid"],
         system_type: SystemType,
-        touchdown_mode: Optional[
-            Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
-        ] = None,
-        collapsed_weak_layer_kR: Optional[float] = None,
+        touchdown_mode: Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
+        | None = None,
+        collapsed_weak_layer_kR: float | None = None,
     ):
         """
         Provide equations for free (pst) or infinite (skiers) ends.

weac/logging_config.py

@@ -4,10 +4,9 @@ Logging configuration for weak layer anticrack nucleation model.
 
 import os
 from logging.config import dictConfig
-from typing import Optional
 
 
-def setup_logging(level: Optional[str] = None) -> None:
+def setup_logging(level: str | None = None) -> None:
     """
     Initialise the global logging configuration exactly once.
     The level is taken from the env var WEAC_LOG_LEVEL (default WARNING).

weac/utils/snowpilot_parser.py

@@ -20,7 +20,6 @@ The `the column length` is the column length of the PropSawTest.
 """
 
 import logging
-from typing import List, Tuple
 
 import numpy as np
 from snowpylot import caaml_parser
@@ -47,28 +46,28 @@ class SnowPilotParser:
     def __init__(self, file_path: str):
         self.snowpit: SnowPit = caaml_parser(file_path)
 
-    def extract_layers(self) -> Tuple[List[Layer], List[str]]:
+    def extract_layers(self) -> tuple[list[Layer], list[str]]:
         """Extract layers from snowpit."""
         snowpit = self.snowpit
-        # Extract layers from snowpit: List[SnowpylotLayer]
-        sp_layers: List[SnowpylotLayer] = [
+        # Extract layers from snowpit: list[SnowpylotLayer]
+        sp_layers: list[SnowpylotLayer] = [
             layer
             for layer in snowpit.snow_profile.layers
             if layer.depth_top is not None
         ]
         sp_layers = sorted(sp_layers, key=lambda x: x.depth_top[0])  # type: ignore
 
-        # Extract density layers from snowpit: List[DensityObs]
-        sp_density_layers: List[DensityObs] = [
+        # Extract density layers from snowpit: list[DensityObs]
+        sp_density_layers: list[DensityObs] = [
             layer
             for layer in snowpit.snow_profile.density_profile
             if layer.depth_top is not None
         ]
         sp_density_layers = sorted(sp_density_layers, key=lambda x: x.depth_top[0])  # type: ignore
 
-        # Populate WEAC layers: List[Layer]
-        layers: List[Layer] = []
-        density_methods: List[str] = []
+        # Populate WEAC layers: list[Layer]
+        layers: list[Layer] = []
+        density_methods: list[str] = []
         for _i, layer in enumerate(sp_layers):
             # Parameters
             grain_type = None
@@ -213,14 +212,14 @@ class SnowPilotParser:
         self,
         layer_top_mm: float,
         layer_bottom_mm: float,
-        sp_density_layers: List[DensityObs],
+        sp_density_layers: list[DensityObs],
     ) -> float | None:
         """Find density measurements that overlap with the given layer depth range.
 
         Args:
             layer_top_mm: Top depth of layer in mm
             layer_bottom_mm: Bottom depth of layer in mm
-            sp_density_layers: List of density observations
+            sp_density_layers: list of density observations
 
         Returns:
             Average density from overlapping measurements, or None if no overlap
@@ -273,8 +272,8 @@ class SnowPilotParser:
         return None
 
     def extract_weak_layer_and_layers_above(
-        self, weak_layer_depth: float, layers: List[Layer]
-    ) -> Tuple[WeakLayer, List[Layer]]:
+        self, weak_layer_depth: float, layers: list[Layer]
+    ) -> tuple[WeakLayer, list[Layer]]:
         """Extract weak layer and layers above the weak layer for the given
         depth_top extracted from the stability test."""
         depth = 0

{weac-3.0.1.dist-info → weac-3.0.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: weac
-Version: 3.0.1
+Version: 3.0.2
 Summary: Weak layer anticrack nucleation model
 Author-email: 2phi GbR <mail@2phi.de>
 License-Expression: MIT
@@ -35,6 +35,9 @@ Requires-Dist: traitlets>=5.14.3; extra == "interactive"
 Provides-Extra: docs
 Requires-Dist: sphinx; extra == "docs"
 Requires-Dist: sphinxawesome-theme; extra == "docs"
+Provides-Extra: build
+Requires-Dist: build; extra == "build"
+Requires-Dist: twine; extra == "build"
 Provides-Extra: dev
 Requires-Dist: nbclient>=0.10.0; extra == "dev"
 Requires-Dist: nbconvert>=7.16.4; extra == "dev"
@@ -55,7 +58,6 @@ Requires-Dist: pycodestyle>=2.11.1; extra == "dev"
 Requires-Dist: black>=24.4.0; extra == "dev"
 Requires-Dist: isort>=5.13.0; extra == "dev"
 Requires-Dist: bump-my-version; extra == "dev"
-Requires-Dist: build; extra == "dev"
 Dynamic: license-file
 
 <!-- LOGO AND TITLE-->
@@ -128,12 +130,13 @@ Dynamic: license-file
 
 1. [About the project](#about-the-project)
 2. [Installation](#installation)
-3. [Usage](#usage)
-4. [Roadmap](#roadmap)
-5. [Release history](#release-history)
-6. [How to contribute](#how-to-contribute)
-7. [License](#license)
-8. [Contact](#contact)
+3. [Development Setup](#development-setup)
+4. [Usage](#usage)
+5. [Roadmap](#roadmap)
+6. [Release history](#release-history)
+7. [How to contribute](#how-to-contribute)
+8. [License](#license)
+9. [Contact](#contact)
 
 <!-- ABOUT THE PROJECT -->
 ## About the project
@@ -189,6 +192,100 @@ Needs (runtime dependencies are declared in [pyproject.toml](https://github.com/
 - [Snowpylot](https://github.com/connellymk/snowpylot) ≥ 1.1.3
 
 
+<!-- DEVELOPMENT SETUP -->
+## Development Setup
+
+This project uses [uv](https://github.com/astral-sh/uv) for fast Python package management and project handling.
+
+### Installing uv
+
+Install uv following the [official installation guide](https://github.com/astral-sh/uv#installation):
+
+```bash
+# On macOS and Linux
+curl -LsSf https://astral.sh/uv/install.sh | sh
+
+# Using pip (alternative)
+pip install uv
+```
+
+### Setting up the development environment
+
+Clone the repository and set up the development environment:
+
+```bash
+git clone https://github.com/2phi/weac
+cd weac
+
+# Install Python 3.12+ if not already available
+# uv will automatically use the version specified in .python-version
+
+# For basic setup (if only running the package):
+uv sync
+
+# For development (recommended for contributors):
+uv sync --extra dev
+
+# Activate the virtual environment
+source .venv/bin/activate  # On Windows: .venv\Scripts\activate
+```
+
+### Running tests
+
+Run the test suite using uv:
+
+```bash
+# Run all tests
+uv run python tests/run_tests.py
+
+# Or use pytest directly (if installed)
+uv run pytest
+```
+
+### Code formatting and linting
+
+This project uses [ruff](https://github.com/astral-sh/ruff) for fast Python linting and formatting:
+
+```bash
+# Format code
+uv run ruff format .
+
+# Check for linting issues
+uv run ruff check .
+
+# Fix auto-fixable linting issues
+uv run ruff check . --fix
+```
+
+### Building the package
+
+Build the package for distribution:
+
+```bash
+# Build wheel and source distribution
+uv build
+
+# Install in editable mode for development
+uv pip install -e .
+```
+
+### Additional uv commands
+
+```bash
+# Update dependencies
+uv sync --upgrade
+
+# Add a new dependency
+uv add package-name
+
+# Add a development dependency
+uv add --dev package-name
+
+# Show environment info
+uv run python --version
+uv run pip list
+```
+
 <!-- USAGE EXAMPLES -->
 ## Usage
 

weac-3.0.2.dist-info/RECORD

@@ -0,0 +1,32 @@
+weac/__init__.py,sha256=mOpFbjyMROJ4BdyLMOUvSfVmzVGLZTYXeIqvGZXGtl4,76
+weac/constants.py,sha256=BX8ifZhFciCuzzako1p-2Wh5CWWzR3cPdvyu501UtUs,1194
+weac/logging_config.py,sha256=GJn4Fs80dDRpbPkSPByCmIOLkEI01dbeuMMSv_IDgoM,1103
+weac/analysis/__init__.py,sha256=twF9OYrZiVmBtw6Lyzf5QUhLZNtqFFM2H3OuaYCog1k,486
+weac/analysis/analyzer.py,sha256=AR1yxM5trK1c-_lTEVn9e6e96IhWVWSRgHgQcaLy4Io,27089
+weac/analysis/criteria_evaluator.py,sha256=WY0FQnfzDfN0h1lI5b4LCUV0bwnvUAVDmS4mKCHGbmo,43599
+weac/analysis/plotter.py,sha256=6KzqiXWIDTIFE6qB947eXhtvEAuPp8B4P5idoSX8OwY,66120
+weac/components/__init__.py,sha256=94WIUVjPI3U-grN_7v0oi9bZB3Z8rCBcx9wJPOwJBAQ,439
+weac/components/config.py,sha256=tnOnJ0M-_knZBhdr052nDyyFFAZN0f2hQ68XRuXG6d8,869
+weac/components/criteria_config.py,sha256=f2agU7nXWURFAw8_68igiSk5aICUxwaB9u3Qasit-Q0,2909
+weac/components/layer.py,sha256=sF47006ORRmNnHW7QMk9h-M0vQQR6TL_PKOvInjdNsk,10560
+weac/components/model_input.py,sha256=oyX4p7PgaAtUXXtpvV8XO_KxwOXtG1TB0x7rXxNfCFA,3330
+weac/components/scenario_config.py,sha256=Tam-m9DQtdjmTm-lKQ8Dcjqje04ttJS2X3v_Nn7AUHQ,2563
+weac/components/segment.py,sha256=F279KcAAkRuJKWav_BZ4BanO96WZm4KXtKHinFZki7s,941
+weac/core/__init__.py,sha256=pRyCKD8XD3qXVUWtFG7N3cS91P5x5d8Jpr1hMEgxQ2U,233
+weac/core/eigensystem.py,sha256=b7KXi2weCY9IVlH_7lCTXzKSx_pLdWY-x8BPjHL5nKo,13736
+weac/core/field_quantities.py,sha256=ci7MvhJ4aYdbW6xxH8vHVgWtk5iypCYv6dZ6KjFNvt8,8964
+weac/core/scenario.py,sha256=vHHe-JDWc-8ECVfy6lFGVtJmIIKRc0J9iUEGR8zckdw,6050
+weac/core/slab.py,sha256=YwbAf7ogYDn3p0nXfj3IjVqwA6ro3jBSg2wXb_io_gQ,5125
+weac/core/slab_touchdown.py,sha256=k30FyvclH_Q3mCZNTbJkNhSenv3ja6qSNhhkmi5p04A,14009
+weac/core/system_model.py,sha256=HAZuquY6p1wf5QuUPA6VOVDIICeCxOaeXqk1dSfEzNY,15513
+weac/core/unknown_constants_solver.py,sha256=pzE7soTo7v-jdDECa21pOWRGBp8YeMU3k7egJOemfSY,17766
+weac/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+weac/utils/geldsetzer.py,sha256=DxvpqmWulBacl-mGeAmAuJcv2xqFBx92PEfNC3aeDzk,3500
+weac/utils/misc.py,sha256=lGz0IDDJ_3nvYjSkivPJ5Xscl1D_AmvQLSjaL7SUbKs,3674
+weac/utils/snow_types.py,sha256=eX9-5La6Oom7zh6pg5JZ4MZ6nLdWdc7RoUzm5e6b9w8,1483
+weac/utils/snowpilot_parser.py,sha256=drOLnUg7uzMJMO9qU9x79xZkqQfFPCCQonVpU9S0x8U,13682
+weac-3.0.2.dist-info/licenses/LICENSE,sha256=ojZPWKFHbFGDrlNOvuAKGH9WcKhpLHWZPcQ4SzhK91M,1082
+weac-3.0.2.dist-info/METADATA,sha256=p-UwzxwB89KHKa1zQIahEoGoymDvmhzRAieD3xzDPt0,25586
+weac-3.0.2.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+weac-3.0.2.dist-info/top_level.txt,sha256=8tyXUHPFU4Ba_5kPtpwvXo5l6GjJmOnODVBJFygpdeE,5
+weac-3.0.2.dist-info/RECORD,,

weac-3.0.1.dist-info/RECORD

@@ -1,32 +0,0 @@
-weac/__init__.py,sha256=YrbBczmXbMc2L2DRklNEOFxlfY0Z9FLDnpW8ZhbaIww,76
-weac/constants.py,sha256=BX8ifZhFciCuzzako1p-2Wh5CWWzR3cPdvyu501UtUs,1194
-weac/logging_config.py,sha256=BFphsxo5s_7E6rPOO2PbOfu0Wy8zHD4AUnnJdqvrn5I,1134
-weac/analysis/__init__.py,sha256=twF9OYrZiVmBtw6Lyzf5QUhLZNtqFFM2H3OuaYCog1k,486
-weac/analysis/analyzer.py,sha256=AR1yxM5trK1c-_lTEVn9e6e96IhWVWSRgHgQcaLy4Io,27089
-weac/analysis/criteria_evaluator.py,sha256=oIVL7Mq3fjyiVx_ms-miULdHE4-EYE2OXoMalM9eU1Q,43458
-weac/analysis/plotter.py,sha256=UdiPZKb8vzZAdB1tjVUNfV2INtZaqEiJ7EA8Y-dzxoM,66056
-weac/components/__init__.py,sha256=94WIUVjPI3U-grN_7v0oi9bZB3Z8rCBcx9wJPOwJBAQ,439
-weac/components/config.py,sha256=tnOnJ0M-_knZBhdr052nDyyFFAZN0f2hQ68XRuXG6d8,869
-weac/components/criteria_config.py,sha256=f2agU7nXWURFAw8_68igiSk5aICUxwaB9u3Qasit-Q0,2909
-weac/components/layer.py,sha256=sF47006ORRmNnHW7QMk9h-M0vQQR6TL_PKOvInjdNsk,10560
-weac/components/model_input.py,sha256=dp695ZMXkSXzJaD-ih9XO0QhnFRPKn8zUJY7g-bomhE,3354
-weac/components/scenario_config.py,sha256=Tam-m9DQtdjmTm-lKQ8Dcjqje04ttJS2X3v_Nn7AUHQ,2563
-weac/components/segment.py,sha256=F279KcAAkRuJKWav_BZ4BanO96WZm4KXtKHinFZki7s,941
-weac/core/__init__.py,sha256=pRyCKD8XD3qXVUWtFG7N3cS91P5x5d8Jpr1hMEgxQ2U,233
-weac/core/eigensystem.py,sha256=I4AIGD_uDVKIGv_sv9AXbIXhJ4Eol64BLB-lZF8X3jw,13735
-weac/core/field_quantities.py,sha256=ci7MvhJ4aYdbW6xxH8vHVgWtk5iypCYv6dZ6KjFNvt8,8964
-weac/core/scenario.py,sha256=f5EZaeFJWMFN_JNOrzVDv5_CeRruHn6jSjGbV9fdxhY,6076
-weac/core/slab.py,sha256=nidM3FvTwtj30LmM4OxOfPwqm5TpPzJrH5UvQR6w4no,5150
-weac/core/slab_touchdown.py,sha256=D0mVVeRktFpjywuhEiIZHOVJflslfKX3oesx1mMQGxg,14025
-weac/core/system_model.py,sha256=XMovKslSTvoNHyxeL7zTFMhAqVOTHnXD7MUPkaGSl9Q,15577
-weac/core/unknown_constants_solver.py,sha256=cDXMVKIdhA8AWZQbJnNpM19pvLllvIVc6Mw6xaZGqYw,17839
-weac/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-weac/utils/geldsetzer.py,sha256=DxvpqmWulBacl-mGeAmAuJcv2xqFBx92PEfNC3aeDzk,3500
-weac/utils/misc.py,sha256=lGz0IDDJ_3nvYjSkivPJ5Xscl1D_AmvQLSjaL7SUbKs,3674
-weac/utils/snow_types.py,sha256=eX9-5La6Oom7zh6pg5JZ4MZ6nLdWdc7RoUzm5e6b9w8,1483
-weac/utils/snowpilot_parser.py,sha256=5fKC0FKcEJvfFQfZbPW5cHqb_nyzIltweq7Bdz5-a0A,13713
-weac-3.0.1.dist-info/licenses/LICENSE,sha256=ojZPWKFHbFGDrlNOvuAKGH9WcKhpLHWZPcQ4SzhK91M,1082
-weac-3.0.1.dist-info/METADATA,sha256=-14ZLdvT2GKGX5n8XLJbd3wxPUHvj8Fi5hg85GfdJPs,23666
-weac-3.0.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-weac-3.0.1.dist-info/top_level.txt,sha256=8tyXUHPFU4Ba_5kPtpwvXo5l6GjJmOnODVBJFygpdeE,5
-weac-3.0.1.dist-info/RECORD,,