function [cp,lambda,M,rho,eta] = updtProp(T,p)
  global p_i;
  global hfTP;
  global dT;
  cp = cp_cpt(T,p);
  lambda = lambda_cpt(T,p);
  [M,rho] = molarM_cpt(T,p);
  eta = eta_cpt(T,p);
%eta*=10;
if 1
cp = 1005*ones(80,1);
lambda = 0.025*ones(80,1);
M = 29*ones(80,1);
rho = p/8314.*M./T;
eta = ones(80,1)*1e-5;
endif
endfunction

if 0
T = (300:10:50000)';
p=100*100000*ones(rows(T),1);
[cp,lambda,M,rho,eta] = updtProp(T,p);
logrho = log(rho);
logT = log(T);
plot(logT,logrho)
endif
