This set of Matlab codes have been developed to compute the evolution of pore fluid pressure, reaction progress, stress and strain for a dehydrating antigorite layer in uniaxial deformation. The problem is written in the form of a coupled set of Ordinary Differential Equations; please refer to the manuscript for an in-depth explanation of the governing equations. For a simple summary, you can type:
publish mainq.mand have a look at the published html/pdf file.
Note: This code is not meant to be very versatile and is quite specific to the problem at hand.
Just clone the repository or copy the files into your working folder. Done !
You will need to have the package water95 in your Matlab path.
The main script mainq.m shuold run smoothly and produce the requested output. All the scripts and functions should be reasonably documented so the command help ought to produce something readable.
Please type publish mainq.m and you should get a decent html or pdf summary of what the code does.
This code has been written by Nicolas Brantut (Rock and Ice Physics and Seismological Laboratory, University College London, UK) and Ioannis Stefanou (Laboratoire Navier, Ecole des Ponts-et-Chaussees, Marne-la-Vallee, France), based on Ioannis’ complete solution originally written in Mathematica.
This repository is an inseparable companion to the following manuscript:
Brantut N., I. Stefanou and J. Sulem, Dehydration-induced instabilites at intermediate depths in subduction zones, submitted to J. Geophys. Res..