Published August 2, 2022 | Version v1
Software Open

Local variation of mechanical parameters in (Griggs-type) viscous deformation experiments (Matlab code)

  • 1. Institute of Earth Sciences, Heidelberg University, Germany
  • 2. Institute of Geosciences, Johannes Gutenberg University Mainz, Germany


Finite difference numerical model to calculate the local mechanical variables (e.g. stress, strain rate) within a Griggs-type experimental assembly. The model solves the Stokes’ equations for slow, viscous, incompressible flow. The code uses the MATLAB® direct solver, employs a staggered grid for discretization and the marker in cell method for model advection (Gerya & Yuen, 2003; Gerya, 2009).

The code was used to produce the data and figures for:
Cionoiu, S., Moulas, E. and Tajčmanová, L. (2019) Impact of interseismic deformation on phase transformations and rock properties in subduction zones. Scientific Reports 9, 19561
Cionoiu, S., Moulas, E., Stünitz H. and Tajčmanová, L. (2022) Locally resolved stress-state in samples during experimental deformation: insights into the effect of stress on mineral reactions, Journal of Geophysical Research: Solid Earth




Files (163.7 kB)

Name Size Download all
54.4 kB Download
95.8 kB Download
11.2 kB Download
2.4 kB Preview Download

Additional details

Related works

Journal article: 10.1016/j.pepi.2003.09.006 (DOI)
Journal article: 10.1038/s41598-019-56130-6 (DOI)