@staticmethod
def calc_racette_co2(nd, no, kd, ko, pop_avg_rdil):
    """Calculate CO2 production in mol/hr using the equation of Racette (1994)
                from dilution spaces and isotope turnover rates.
           :param nd: deuterium dilution space in mol
           :param no: oxygen dilution space in mol
           :param kd: deuterium turnover rate in 1/hr
           :param ko: oxygen turnover rate in 1/hr
           :return co2prod: co2 production rate in mol/hr
        """
    r_dil = (pop_avg_rdil + 1.034) / 2
    n = (no / 1.01 + nd / (1.01 * r_dil)) / 2
    co2_prod = n / 2.078 * (1.01 * ko - 1.01 * kd * r_dil
        ) - 0.0245 * n * 1.05 * (1.01 * ko - 1.01 * kd * r_dil)
    return co2_prod
