version_drift4aspect_pluginsNprms

Porting the additional ASPECT plugins and FastScape-coupled parameter files from Lei et al. (2026, EPSL 120073) forward to a current ASPECT + FastScape build, and logging every change so collaborators can review the port and flag bugs.

Study: Geodynamic controls on the habitability zone of shallow coral reefs at rifted continental margins — Zhibin Lei, Shaohong Xia, Sascha Brune, Michaël Pons, Chenglong Zhang, Tao Gou, Changrong Zhang, Chaoyan Fan.

Why this repo exists

The published v1.0 plugins and .prm files were written against an older ASPECT (circa 2.3.0-pre) and an older FastScape. Running them on a current ASPECT 3.1.0-pre + FastScape build requires two kinds of fixes:

1. Plugin source — small ASPECT 3.x API changes and build-path portability.
2. Parameter files — FastScape and ASPECT parameter names/structure drifted between versions ("version drift"); the old .prm keys must be remapped.

The first commit (EPSL_120073_v1.0) is the pristine published baseline. Every commit after it is one reviewable change with a rationale in its message.

Provenance: starting point vs. current target

Component	Original (paper / v1.0)	Currently testing against
ASPECT	Djneu/aspect branch fault_analysis (~2.3.0-pre), following Neuharth et al. (2022) — https://github.com/Djneu/aspect/tree/fault_analysis	geodynamics/aspect branch main, VERSION 3.1.0-pre, commit 2269666ef (2026-06-16) — https://github.com/geodynamics/aspect
FastScape	fastscape-lem/fastscapelib-fortran commit 18f2588	Djneu/fastscapelib-fortran branch coupling_dev, commit 887e3ea (2024-10-15), fastscapelib v2.9.0dev — https://github.com/Djneu/fastscapelib-fortran
deal.II	9.2.0	9.7.0
Plugins	Written by Anne Glerum (Neuharth et al., 2022): LAB perturbation, plastic-strain noise, composition-based mesh refinement	Same plugins, ported (see commits)

Note on FastScape fork: the current build uses the Djneu fork, which already ships Named_VTK.f90 — unlike the original workflow, it does not need to be copied from ASPECT.

What's in here

additional_plugins/
  comp_plugin/        # composition tagging / mesh-refinement plugin
  rift_plugins/       # lithosphere rift (comp / temp / geom) + rift postprocessor
prmfiles/
  M1..M6, S1..S7      # the 13 published parameter files (2.x)
  M1.wide-nosp1Myr.prm# fully ported M1 prototype that runs on 3.1.0-pre
  version_drift_refs/ # cited 2.x -> 3.x parameter mapping + working 3.x reference prms
readme.txt            # original v1.0 study readme (kept verbatim)

Port status

• ✅ Plugins compile and load on ASPECT 3.1.0-pre (API + relative build paths).
• ✅ FastScape 2.x→3.x parameter drift documented with source citations in prmfiles/version_drift_refs/README_version_drift_evidence.md.
• ✅ M1.wide-nosp1Myr.prm ported and running on 3.1.0-pre.
• ⏳ The remaining 12 published .prm files still need the same parameter remapping (later rounds).

Reproducibility caveat: FastScape v2.9.0dev here vs. commit 18f2588 in the paper — parameter names map cleanly and values are preserved, but marine/erosional physics may not be bit-identical. Sanity-check against the published figures.

Build / run

See readme.txt for the original install recipe. Each .prm loads the plugins via its first line, set Additional shared libraries = ...; edit that line to point at your own plugin build/ directory before running (these machine-specific paths are intentionally not committed).