fix conservation of 2p tracer
Add a function equilibrateTracer which is called after each
update of the flow solution. It adapts the tracer mole fraction
according to the change in phase distribution such that the amount
of tracer in each element is conserved.
Fixes #1071.