[PL/TH2M] Add/correct residuum output

The residuum consists of heat flux, gas and liquid mass flow rates,
and nodal forces.

While the residuum for the Newton scheme is readily available,
some computations are necessary in the Picard scheme.