[WIP] Backup current_linearization_point for use in composition reaction_terms

[anne-glerum]

Create a copy of the current_linearization_point at the beginning of each timestep to compute compositional reaction_terms for fields that are updated in an iterative Advection scheme.

I'd like to discuss with @gassmoeller before actually using the backup.

For all pull requests:

• I have followed the instructions for indenting my code.

For new features/models or changes of existing features:

• I have tested my new feature locally to ensure it is correct.
• I have created a testcase for the new feature/benchmark in the tests/ directory.
• I have added a changelog entry in the doc/modules/changes directory that will inform other users of my change.