- Extracted the Mohr-Coulomb part from the CalcStressAndFailure function in Relaxation
        new functions:
                bool Relaxation::Mohr_Coulomb_Failure();
                Real Relaxation::Mohr_Coulomb_Breaking_Probability(Particle *prtcl, int nb);

- void ExperimentSettings::Couple_Mohr_Coulomb_Failure_Criterion_With_Griffith_Failure(bool activate);
        Adds the Griffith criterion for tensile failure to the Mohr-Coulomb part.
        In effect, this prevents tensile fracturing if sigma1-sigma3 > 4 * breaking_strength

- void Relaxation::Magmatic_Intrusion(Real sigma_3_inc, Real max_viscosity, Real xmin, Real xmax, Real ymin, Real ymax, Real zmin, Real zmax);
        this function decreases the viscosity of a volume as  a function of the tensile stress

Added several fun