THM: freezing related updates
See individual commits for changes.
- refactor(THM): consistent Jacobian signs
- test(THM): avoid ice formation at low positive temperatures
- feat(THM) Add epsilon_0 output
- feat(THM): ice volume fraction avg output
- feat(THM): compute ice vol. fraction in initialization
Edited by Dmitri Naumov