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

I want to solve combination of momentum and stokes equations to model visco-elastic deformation. "Viscous" media should undergo viscous deformation, and "non-viscous" media should undergo elastic deformation that recovers when load is "switched off".

Here is the equations that I take from "Introduction to numerical geodynamic modeling" http://imgur.com/a/EUivt