The code calculates Griggs Rig experiment local variations 

Numerical solving of 2D stream function using Staggered Grid - variable viscosity
working for Simple shear, pure shear, inclusions in PS
based on: Gerya, 2009. Introduction to numerical geodynamic modeling
this script is available at Zenodo: https://doi.org/10.5281/zenodo.6951760
input file is generated with: calc_Griggs_v1_0.m
by Sebastian Cionoiu 

License: 
Creative Commons Attribution 4.0 International

Files:

calc_Griggs_v1_0.m 
-> main calculation file, produces a matlab file that is than plotted with the other two scripts
-> having the LRs folder in the same base folder speeds up the calculation for pre-computed matrix resolutions
-> calculation time on a desktop PC: 551x815: 1.5h / 100 time-step, 300x501: 20 min / 100 time-step
-> too low resolution (< 301x501) will create problems with the shear motion in the Pt separators (need to be thicker)

LRs (folder)
-> calc_Griggs_v1_0.m requires sparse matrices, which need time to assemble. Preambled matrices can be used, if the boundary conditions of the code are not changed (before or during run). Newly assemble matrices will be saved to this folder.
-> see lines 230-244 in calc_Griggs_v1_0.m 

plot_Griggs_v1_0.m
-> main plotting script, requires output of calc_Griggs_v1_0.m 
-> used for almost all figures in Cionoiu et al. 2022 and Cionoiu et al. 2019

plotTS_Griggs_v1_0.m
-> time segment extraction and plotting script, requires output of calc_Griggs_v1_0.m 
-> used for figure 8 in Cionoiu et al. 2022


References:
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., Tajčmanová, L., Moulas, E. and Stünitz H. (under review, 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

Gerya, T. V. (2009). Introduction to numerical geodynamic modelling. Cambridge University Press.

Gerya, T. V., & Yuen, D. A. (2003). Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Physics of the Earth and Planetary Interiors, 140(4), 293–318. https://doi.org/10.1016/j.pepi.2003.09.006