Consider adding an optional damping term to nr_post_update

[Goobley]

We can damp the Newton-Raphson correction if it's too large -- it's still self-consistent, obviously the process isn't converged till we're done with damping; like the Armijo line search.