function cpc = cpcmean(T)
% Simpson and TenWolde (1999)
T=T+273.15;
cpc=1.43+0.355e-3*(T)-7.32e4/(T)^2;
end