[sfepy] Re: Wrong results in lineral_viscoelastic.py from examples. What may cause this?

What if I want to check each step for plastic deformation? For example, I want to check stress on each step and if stress at some parts of the domain is over certain threshold I want to drastically lower viscosity to emulate fracture. Should I do it manually or there is a way to do it via sfepy?