function cpb = cpbmean(T)
% Simpson and TenWolde (1999)
T=T+273.15;
cpb=0.1031+0.003867*T;
end