ShearToast is a MATLAB 1D code. The name reflects on that is measures frictional heat in a shear zone.
The code handles:
-
temperature -
pressure
-
stress
-
power law breakdown
-
viscosity
-
latent heat in melt generation
-
confined pressure
-
fluid dependent mineral phase trans.
-
ductile and brittle deformation
-
normal and ‘slow’ earthquakes
-
a large range of rocks and ice
ShearToast is a simplified and thus much faster version of LiToastPhere. So far I use the code as my worktool to compare rockdata to geodynamic models at smaller and more specific scale.
I am particularly exited about how the code can handle a new hypothesis of confined pressure resulting from melt expansion. The plot above illustrates this point. It shows rising pressure as a results of rising differential stress in a shear zone. Eventually the stress (and thus the pressure) drops as an effect of shear heating. The heat causes the rocks to melt in a confined cell, which is modeled to cause an short-lived pressure spike. The heat spreads to the sidewalls to decay stress also outside the shear zone. You can follow the PTtime evolution is here plotted in a standard diagram.
Once the code is published I will be happy to distribute it, but so far you can contact me if you have interesting problems that you think we can address by it.
