function [ es ] = satvaph2o( T )
% Saturation vapor pressure in Pa at temperature T in K
% Using wiki formula
% T can be a vector

Tc=T-273.15; % temperature in 'C
es=610.94.*exp(17.625.*Tc./(Tc+243.04));

end

