Minor modification, set the default od allow_dissipation_contact to true and a condition on the if statement for damping the initial unstable zone.