def isocap(CO2, pH, K1, K2, KB, KW, TB):
    """Isocapnic quotient of HDW18, Eq. 8."""
    h = 10.0 ** -pH
    return ((K1 * CO2 * h + 4 * K1 * K2 * CO2 + KW * h + h ** 3) * (KB + h) **
        2 + KB * TB * h ** 3) / (K1 * CO2 * (2 * K2 + h) * (KB + h) ** 2)
