function [kco2,kh2o] = GibbsEngCalcs(temperature,gibbseng)

% This function loads the gibbs energy data and calculates equilibrium
% constants.

%Gibbs free energy
    gco2=interp1(gibbseng(1:end,1),gibbseng(1:end,2),temperature)*10^3;
    gh2o=interp1(gibbseng(1:end,1),gibbseng(1:end,3),temperature)*10^3;
    gco=interp1(gibbseng(1:end,1),gibbseng(1:end,4),temperature)*10^3;

%Equilibrium constants
    kco2=exp((-(gco2-gco))*10^3/(8314.46*temperature));
    kh2o=exp((-gh2o)*10^3/(8314.46*temperature));
end