In problem add actions_after_implicit_timestep_and_error_estimation()

New actions function called both unsteady newton solve functions.

This is needed if your actions_after_implicit_timestep() modify the
solution in a way that affects the error estimate.