Add fast "pre-solver" for fixed SFCC parameter settings

SFCCPreSolver constructs in initialconditon! a linear interpolant to approximate the enthalpy freeze curve

\theta(H)

. This is only valid when all freeze curve parmaeters (including the water content) can be assumed constant.