When updating the strain compositional fields in the strain_depended rheology module, it uses the current strain. I don't think this is conceptually the best idea, since I think it should but the accumulated strain up to this point, but it is also numerically very difficult to solve since the whole system can become extremely nonlinear.
This pull request changes that instead of the current strain-rate to the old strain-rate (previous time-step) is added to the strain fields.
This is a work in progress, and I still need to do some testing how this would work out, but I would like to get some early feedback on whether this makes sense to other people or not.
When updating the strain compositional fields in the strain_depended rheology module, it uses the current strain. I don't think this is conceptually the best idea, since I think it should but the accumulated strain up to this point, but it is also numerically very difficult to solve since the whole system can become extremely nonlinear.
This pull request changes that instead of the current strain-rate to the old strain-rate (previous time-step) is added to the strain fields.
This is a work in progress, and I still need to do some testing how this would work out, but I would like to get some early feedback on whether this makes sense to other people or not.