Output the total mass balance residual
The residual output is changed from max residual to
total mass balance residual. In this way the output
resemble what is actually used as convergence criteria.
Output the scaled mass balance residuals
The mass balance residuals are multiplied by the porevolume divided by
time
Merge pull request #229 from totto82/output_conv

Output the total mass balance residual