goSPL v2026.6.30

Highlights

• Post-processing & analysis toolkit — new gospl-catchment per-basin outflow/water-flux tool, cross-sections, regular-grid + NetCDF export, stratigraphic sections, wedge & ELA extraction, strata-volume; plus a top-level gospl console command and Python entry points.
• Diagnostic glacial-erosion model (#454) — routed discharge → ice thickness → Glen sliding velocity → abrasion / till transport / moraine deposition / meltwater / flexural loading; melt mass-balance knob.
• Spatially-resolved provenance — per-class sink attribution including continental pit/lake sinks, with post-processing fields and [prov] verbosity.
• Losing-stream evaporation — evaporation-coherent transport-limited SPL.
• Dual-lithology strata usability — npstrata validation, surfFineFrac, dedicated infinite-bedrock sentinel, pitInletBias in pit composition.
• Parallel flat-model flexure — FV biharmonic solver ('fem', #471), replacing gFlex + FFT.
• Mesh-native orographic rain — parallel solver on the DMPlex (#474); the regular grid (regdx/latitude/nm/hw) is removed.
• Parallel correctness & performance — partition-stable flow-accumulation KSP, partition-invariant pit-fill/flat-routing, marine Picard/Schur solvers, pit-label union-find, profiler + HPC scaling harness.

⚠️ Behavior changes / migration

• Open (o) boundaries are now "non-rising" (#473, f8c1f23): each step an open edge is set to min(current, min(interior neighbours)) — it follows incision/subsidence down but aggradation can no longer raise it (previously it tracked the neighbour average). If you relied on the old non-draining behavior, switch that edge to w (true no-flux wall, new in #473).
• Boundary bc strings use the NESW letter convention o/f/w/c (open/fixed/wall/cyclic). Legacy digit strings are remapped in place (0→o, 1→f); check your inputs match the intended N,E,S,W edges.
• gFlex + FFT flexure removed — 'fem' is the only flat-model flexure.
• Regular-grid orography removed — orographic rain is now mesh-native; drop regdx/latitude/nm/hw from configs.

Notable fixes

• Multiple parallel-correctness fixes: flow-accumulation KSP partition-stability, partition-invariant drainage graph & flat-routing tie-break, several np>1 deadlocks, closed-sink deposition guards.
• Analytical hillslope-diffusion benchmark reconciled to the non-rising open-edge semantics.