Refactor density computation. Limit Cp so that loads are always outputed even if solution is not locally converged.