def _update_diagnostics(self, newstate):
pass<EOL>
This method is called each timestep after the new state is computed
        with the implicit solver. Daughter classes can implement this method to
        compute any diagnostic quantities using the new state.