function [Kf]=Kf(T) 
  %fluid thermal conductivity [W/m*C] at T [celsius] 
 
 Tsat=[0:10:370,374.15];
  %saturation temperature in celsius

 Kliq=[0.569,0.587,0.603,0.618,0.631,0.643,0.653,0.662,0.670,0.676,0.681,0.684,0.687,0.688,0.688,0.687,0.684,0.681,0.677,0.671,0.664,0.657,0.648,0.639,0.628,0.616,0.603,0.589,0.574,0.558,0.541,0.523,0.503,0.482,0.460,0.434,0.397,0.340,0.240];

Kf=interp1(Tsat,Kliq,T);
