band. Made sure slip direction was dialtatant. Track residual cohesion and
do not allow this to go below 0. Tried to create a routine to track
stiffness on step where cohesion does go to zero, but this gives only
linear global convergence. Updated to handle reverse loading and guarantee
that jump is positive along slip direction.
