Three-dimensional Underworld/UWGeodynamics experiments used to test how slab steps and material properties affect far-field stress transmission along subduction zones. The scripts support the study “The Role of Slab Steps on Tectonic Loading Along Subduction Zones: Inferences on the Seismotectonics of the Sunda Convergent Margin.”
The model represents upper mantle, lower mantle, upper-plate crust and mantle, and long/short subducting slab segments. Parameter sweeps vary crustal cohesion, lithospheric thickness, slab-step length, and viscosity contrasts. Production configurations use large three-dimensional meshes and are intended for MPI-enabled HPC systems.
| Directory | Sweep |
|---|---|
Cohesion_test_SubductingPlateCrust/ |
subducting-plate crust cohesion |
Cohesion_test_UpperPlateCrust/ |
upper-plate crust cohesion |
Thickness_test_UpperPlateMantle/ |
upper-plate mantle thickness |
Step_length_test/ |
slab-step length |
Viscosity_test_SubductingPlateMantle:UpperMantle/ |
subducting-plate mantle / upper-mantle viscosity contrast |
Viscosity_test_UpperPlate:UpperMantle/ |
upper-plate / upper-mantle viscosity contrast |
The complete value inventory and a reproducible run-record template are in docs/experiments.md. The colon-containing viscosity directory names are retained for compatibility with the original run archive; quote these paths in shell commands.
These scripts require a compatible Underworld 2, UWGeodynamics, PETSc, MPI, NumPy, and h5py environment. Use the module or container versions provided by the target HPC facility and record their exact versions with each run.
The production scripts default to the original NCI scratch root and commonly define a 512 x 512 x 128 element mesh with 20 particles per cell. Review outputPath, resolution, process count, memory, wall time, and filesystem quota before submission. Do not launch a production configuration on a workstation.
Run the dependency-free repository checks locally:
git clone https://github.com/ZuhairQuakes/HPC-3D-Modelling.git
cd HPC-3D-Modelling
python tools/validate_repository.pyA representative MPI launch, after adapting the script and scheduler allocation, is:
mpirun -np <ranks> python "Step_length_test/uw_subduction_3d_var_dep_upfix_no_pp_25_step_length_330.py"Full validation requires a short cluster smoke test followed by the intended production allocation. Generated HDF5, XDMF, visualization, scheduler, and log files are ignored by Git.
Use GitHub's Cite this repository menu or CITATION.cff:
Zuhair, M., Gollapalli, T., Capitanio, F. A., Betts, P. G., & Graciosa, J. C. (2022). The Role of Slab Steps on Tectonic Loading Along Subduction Zones: Inferences on the Seismotectonics of the Sunda Convergent Margin. Tectonics, 41(9). https://doi.org/10.1029/2022TC007242
See CONTRIBUTING.md for experiment-review requirements. The code is available under the MIT License.