ur=10;
zr=10;
Ta=15;
zt=10;
rh=70;
zq=10;
Pa=1000.;
Ts=2;
A=hfbulktc(ur,zr,Ta,zt,rh,zq,Pa,Ts);

ce=A(11);  % Dalton Number
