function [Cpf]=Cpf(T) 
  %fluid specific heat [kJ/kg*C] at T [celsius] 
 
 Tsat=[0:10:370,374.15];
  %saturation temperature in celsius

 Cpliq=[4.218,4.194,4.182,4.179,4.179,4.181,4.185,4.191,4.198,4.207,4.218,4.230,4.244,4.262,4.282,4.306,4.334,4.366,4.403,4.446,4.494,4.550,4.613,4.685,4.769,4.866,4.985,5.134,5.307,5.520,5.794,6.143,6.604,7.241,8.225,10.07,15,55,1e6];
 %liquid specific heat [kJ/kg*C] at T [celsius] 

Cpf=interp1(Tsat,Cpliq,T);
