find density from temperature and liquid fraction
function rho_Tf input SI.Temperature T "Absolute temperature (K)"; input Real f "Liquid mass fraction"; output SI.SpecificEnthalpy rho "Density (kg/m3)"; end rho_Tf;