%
% Check the function that compute
%       saturation vapor pressure for a given temperature
%
% Input:  T=Temperature [K
% Output: Es=Saturation vapor pressure [Pa]
%
% Author: Xin-Zhong Liang, 2012-08-25
%
function [Es, Es_c] = EsCheck (T)
Es = 611 * ( 10 .^ (7.5 * (T-273.15) ./ (T-35.85)) );
% check the result against another approximation
Es_c = satprs(T);
end