A parallel 3D numerical code that can be used to model various thermo-mechanical geodynamical processes such as mantle-lithosphere interaction, or the dynamics of magmatic systems, for rocks that have visco-elasto-plastic or (compressible) poroelastic rheologies. The code uses a marker-in-cell approach with a staggered finite difference discretization and is built on top of PETSc. A range of (Galerkin) multigrid and iterative solvers are available, for both linear and non-linear rheologies, using Picard and quasi-Newton solvers. LaMEM has been tested on a massively parallel MPI-cluster with 458,752 cores and has checkpointing/restarting capabilities.
- Used by CoE: ChEESE
- Main uses: Geodynamics
- Access Links