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.