of the entries of the global residual vector that do not need
 initial non-equilibrium compensation

Use const references where necessary.

Interface and time loop changes. Processes are in the
next commit.

The integration point writer are defined only on
the bulk mesh of a process. If other meshes (observation,
boundary, etc.) are given for output, the integration
point data must be output only on the main mesh of
a process.

Name time increment as everywhere else 'dt'.

Main reason for modification is the getFlux function, which for
particular processes (ComponentTransport) relies on all solutions being
available.
This modification ripples through all of the code, especially through
the boundary conditions.
Modified `applyNaturalBC`, `preTimestep`, `integrate`, and `assemble`
functions.

preTimestep now accepts a vector of solution together with
a process_id.
In postTimestep only the vector of solutions has been changed.

Pass process_id to FEM::assemble calls.