This repository archives scripts that allow to reproduce results of the submitted manuscript to Geoscientific Model Development. Automatic Dynamic Relaxation (DR) and Powell-Hestenes / Dynamic Relaxation (PHDR) solvers using Finite Difference (FD) and Face-Centered Finite Volume (FCFV) discretisations.
FD and FCFV Poisson solvers with heterogeneous coefficients using the DR method.

FD models based on the PHDR and DR methods: compressible/incompressible deformation, frictional plasticity, two-phase flow
FD and FCFV models of incompressible Stokes flow with many inclusions of variable viscosity.
- Clone/Download the repository
- Open the folder in a dedicated VScode window
- Open Julia's REPL and type to switch to package mode:
] - Install all necessary dependencies, type:
instantiate
Some script rely on ParallelStencil. Scripts that initialize ParallelStencil in 2D and scripts that used 3D cannot be run within the same Julia session because of the macro-based design of ParallelStencil. Instead, Julia needs to be restarted to account for the differrent macros in 2D and 3D.