def Campbell_Thodos(T, Tb, Tc, Pc, M, dipole=None, hydroxyl=False):
Tr = T/Tc<EOL>Tbr = Tb/Tc<EOL>Pc = Pc/<NUM_LIT><EOL>s = Tbr * log(Pc)/(<NUM_LIT:1>-Tbr)<EOL>Lambda = Pc**(<NUM_LIT:1>/<NUM_LIT>)/(M**<NUM_LIT:0.5>*Tc**(<NUM_LIT:5>/<NUM_LIT>))<EOL>alpha = <NUM_LIT> - <NUM_LIT>*s<EOL>beta = <NUM_LIT>*s - <NUM_LIT> + <NUM_LIT>*Lambda**(<NUM_LIT>)<EOL>if dipole:<EOL><INDENT>theta = Pc*dipole**<NUM_LIT:2>/Tc**<NUM_LIT:2><EOL>alpha -= <NUM_LIT> * theta**<NUM_LIT><EOL>beta += <NUM_LIT> * theta**<NUM_LIT><EOL><DEDENT>if hydroxyl:<EOL><INDENT>beta = <NUM_LIT>*s - <NUM_LIT> + <NUM_LIT>*Lambda**(<NUM_LIT>) + <NUM_LIT>*theta**<NUM_LIT><EOL>alpha = (<NUM_LIT>*Tbr - <NUM_LIT> + <NUM_LIT>/Tbr**<NUM_LIT>)*Pc**<NUM_LIT><EOL><DEDENT>Zra = alpha + beta*(<NUM_LIT:1>-Tr)<EOL>Vs = R*Tc/(Pc*<NUM_LIT>)*Zra**(<NUM_LIT:1>+(<NUM_LIT:1>-Tr)**(<NUM_LIT:2>/<NUM_LIT>))<EOL>return Vs<EOL>
r'''Calculate saturation liquid density using the Campbell-Thodos [1]_
    CSP method.

    An old and uncommon estimation method.

    .. math::
        V_s = \frac{RT_c}{P_c}{Z_{RA}}^{[1+(1-T_r)^{2/7}]}

        Z_{RA} = \alpha + \beta(1-T_r)

        \alpha = 0.3883-0.0179s

        s = T_{br} \frac{\ln P_c}{(1-T_{br})}

        \beta = 0.00318s-0.0211+0.625\Lambda^{1.35}

        \Lambda = \frac{P_c^{1/3}} { M^{1/2} T_c^{5/6}}

    For polar compounds:

    .. math::
        \theta = P_c \mu^2/T_c^2

        \alpha = 0.3883 - 0.0179s - 130540\theta^{2.41}

        \beta = 0.00318s - 0.0211 + 0.625\Lambda^{1.35} + 9.74\times
        10^6 \theta^{3.38}

    Polar Combounds with hydroxyl groups (water, alcohols)

    .. math::
        \alpha = \left[0.690T_{br} -0.3342 + \frac{5.79\times 10^{-10}}
        {T_{br}^{32.75}}\right] P_c^{0.145}

        \beta = 0.00318s - 0.0211 + 0.625 \Lambda^{1.35} + 5.90\Theta^{0.835}

    Parameters
    ----------
    T : float
        Temperature of fluid [K]
    Tb : float
        Boiling temperature of the fluid [K]
    Tc : float
        Critical temperature of fluid [K]
    Pc : float
        Critical pressure of fluid [Pa]
    M : float
        Molecular weight of the fluid [g/mol]
    dipole : float, optional
        Dipole moment of the fluid [debye]
    hydroxyl : bool, optional
        Swith to use the hydroxyl variant for polar fluids

    Returns
    -------
    Vs : float
        Saturation liquid volume

    Notes
    -----
    If a dipole is provided, the polar chemical method is used.
    The paper is an excellent read.
    Pc is internally converted to atm.

    Examples
    --------
    Ammonia, from [1]_.

    >>> Campbell_Thodos(T=405.45, Tb=239.82, Tc=405.45, Pc=111.7*101325, M=17.03, dipole=1.47)
    7.347363635885525e-05

    References
    ----------
    .. [1] Campbell, Scott W., and George Thodos. "Prediction of Saturated
       Liquid Densities and Critical Volumes for Polar and Nonpolar
       Substances." Journal of Chemical & Engineering Data 30, no. 1
       (January 1, 1985): 102-11. doi:10.1021/je00039a032.