function Cp = Cp_H2O_g(T)
Cp = (4.395 + -4.186e-3*T + 1.405e-5*T^2 + -1.564e-8*T^3 + 0.632e-11*T^4)*8.3145;
end