def eli_hanley_dense(T, MW, Tc, Vc, Zc, omega, Cvm, Vm):
Cs = [<NUM_LIT>, -<NUM_LIT>, <NUM_LIT>, -<NUM_LIT>,<EOL><NUM_LIT>, -<NUM_LIT>, <NUM_LIT>, -<NUM_LIT>,<EOL><NUM_LIT>]<EOL>Tr = T/Tc<EOL>if Tr > <NUM_LIT:2>:<EOL><INDENT>Tr = <NUM_LIT:2><EOL><DEDENT>Vr = Vm/Vc<EOL>if Vr > <NUM_LIT:2>:<EOL><INDENT>Vr = <NUM_LIT:2><EOL><DEDENT>theta = <NUM_LIT:1> + (omega - <NUM_LIT>)*(<NUM_LIT> - <NUM_LIT>*log(Tr) + (<NUM_LIT> - <NUM_LIT>/Tr)*(Vr-<NUM_LIT:0.5>))<EOL>psi = (<NUM_LIT:1> + (omega-<NUM_LIT>)*(<NUM_LIT>*(Vr-<NUM_LIT>) - <NUM_LIT>*(Vr-<NUM_LIT>)*log(Tr)))*<NUM_LIT>/Zc<EOL>f = Tc/<NUM_LIT>*theta<EOL>h = Vc/<NUM_LIT>*psi<EOL>T0 = T/f<EOL>rho0 = <NUM_LIT>/(Vm*<NUM_LIT>)*h  <EOL>eta0 = <NUM_LIT>*sum([Cs[i]*T0**((i+<NUM_LIT:1>-<NUM_LIT:4>)/<NUM_LIT>) for i in range(len(Cs))])<EOL>k1 = <NUM_LIT>*eta0<EOL>b1 = -<NUM_LIT><EOL>b2 = <NUM_LIT><EOL>b3 = <NUM_LIT><EOL>b4 = <NUM_LIT><EOL>k2 = (b1 + b2*(b3 - log(T0/b4))**<NUM_LIT:2>)/<NUM_LIT>*rho0<EOL>a1 = -<NUM_LIT><EOL>a2 = <NUM_LIT><EOL>a3 = <NUM_LIT><EOL>a4 = -<NUM_LIT><EOL>a5 = <NUM_LIT><EOL>a6 = <NUM_LIT><EOL>a7 = -<NUM_LIT><EOL>k3 = exp(a1 + a2/T0)*(exp((a3 + a4/T0**<NUM_LIT>)*rho0**<NUM_LIT:0.1> + (rho0/<NUM_LIT> - <NUM_LIT:1>)*rho0**<NUM_LIT:0.5>*(a5 + a6/T0 + a7/T0**<NUM_LIT:2>)) - <NUM_LIT:1>)/<NUM_LIT><EOL>if T/Tc > <NUM_LIT:2>:<EOL><INDENT>dtheta = <NUM_LIT:0><EOL><DEDENT>else:<EOL><INDENT>dtheta = (omega - <NUM_LIT>)*(-<NUM_LIT>/T + (Vr-<NUM_LIT:0.5>)*<NUM_LIT>*Tc/T**<NUM_LIT:2>)<EOL><DEDENT>dfdT = Tc/<NUM_LIT>*dtheta<EOL>X = ((<NUM_LIT:1> - T/f*dfdT)*<NUM_LIT>/Zc)**<NUM_LIT><EOL>H = (<NUM_LIT>/MW)**<NUM_LIT:0.5>*f**<NUM_LIT:0.5>/h**(<NUM_LIT:2>/<NUM_LIT>)<EOL>ks = (k1*X + k2 + k3)*H<EOL>theta = <NUM_LIT:1> + (omega - <NUM_LIT>)*(<NUM_LIT> - <NUM_LIT>*log(Tr) - <NUM_LIT>/Tr)<EOL>psi = (<NUM_LIT:1> + (omega-<NUM_LIT>)*(<NUM_LIT> - <NUM_LIT>*log(Tr)))*<NUM_LIT>/Zc<EOL>f = Tc/<NUM_LIT>*theta<EOL>h = Vc/<NUM_LIT>*psi<EOL>T0 = T/f<EOL>eta0 = <NUM_LIT>*sum([Cs[i]*T0**((i+<NUM_LIT:1>-<NUM_LIT:4>)/<NUM_LIT>) for i in range(len(Cs))])<EOL>H = (<NUM_LIT>/MW)**<NUM_LIT:0.5>*f**<NUM_LIT:0.5>/h**(<NUM_LIT:2>/<NUM_LIT>)<EOL>etas = eta0*H*MW/<NUM_LIT><EOL>k = ks + etas/(MW/<NUM_LIT>)*<NUM_LIT>*(Cvm-<NUM_LIT:3>*R/<NUM_LIT>)<EOL>return k<EOL>
r'''Estimates the thermal conductivity of a gas at high pressure as a
    function of temperature using the reference fluid method of Eli and
    Hanley [1]_ as shown in [2]_.

    .. math::
        Tr = min(Tr, 2)

        Vr = min(Vr, 2)

        f = \frac{T_c}{190.4}\theta

        h = \frac{V_c}{9.92E-5}\psi

        T_0 = T/f

        \rho_0 = \frac{16.04}{V}h

        \theta = 1 + (\omega-0.011)\left(0.09057 - 0.86276\ln Tr + \left(
        0.31664 - \frac{0.46568}{Tr}\right) (V_r - 0.5)\right)

        \psi = [1 + (\omega - 0.011)(0.39490(V_r - 1.02355) - 0.93281(V_r -
        0.75464)\ln T_r]\frac{0.288}{Z_c}

        \lambda_1 = 1944 \eta_0

        \lambda_2 = \left\{b_1 + b_2\left[b_3 - \ln \left(\frac{T_0}{b_4}
        \right)\right]^2\right\}\rho_0

        \lambda_3 = \exp\left(a_1 + \frac{a_2}{T_0}\right)\left\{\exp[(a_3 +
        \frac{a_4}{T_0^{1.5}})\rho_0^{0.1} + (\frac{\rho_0}{0.1617} - 1)
        \rho_0^{0.5}(a_5 + \frac{a_6}{T_0} + \frac{a_7}{T_0^2})] - 1\right\}

        \lambda^{**} = [\lambda_1 + \lambda_2 + \lambda_3]H

        H = \left(\frac{16.04}{MW}\right)^{0.5}f^{0.5}/h^{2/3}

        X = \left\{\left[1 - \frac{T}{f}\left(\frac{df}{dT}\right)_v \right]
        \frac{0.288}{Z_c}\right\}^{1.5}

        \left(\frac{df}{dT}\right)_v = \frac{T_c}{190.4}\left(\frac{d\theta}
        {d T}\right)_v

        \left(\frac{d\theta}{d T}\right)_v = (\omega-0.011)\left[
        \frac{-0.86276}{T} + (V_r-0.5)\frac{0.46568T_c}{T^2}\right]

    Parameters
    ----------
    T : float
        Temperature of the gas [K]
    MW : float
        Molecular weight of the gas [g/mol]
    Tc : float
        Critical temperature of the gas [K]
    Vc : float
        Critical volume of the gas [m^3/mol]
    Zc : float
        Critical compressibility of the gas []
    omega : float
        Acentric factor of the gas [-]
    Cvm : float
        Molar contant volume heat capacity of the gas [J/mol/K]
    Vm : float
        Volume of the gas at T and P [m^3/mol]

    Returns
    -------
    kg : float
        Estimated dense gas thermal conductivity [W/m/k]

    Notes
    -----
    Reference fluid is Methane.
    MW internally converted to kg/g-mol.

    Examples
    --------
    >>> eli_hanley_dense(T=473., MW=42.081, Tc=364.9, Vc=1.81E-4, Zc=0.274,
    ... omega=0.144, Cvm=82.70, Vm=1.721E-4)
    0.06038475936515042

    References
    ----------
    .. [1] Ely, James F., and H. J. M. Hanley. "Prediction of Transport
       Properties. 2. Thermal Conductivity of Pure Fluids and Mixtures."
       Industrial & Engineering Chemistry Fundamentals 22, no. 1 (February 1,
       1983): 90-97. doi:10.1021/i100009a016.
    .. [2] Reid, Robert C.; Prausnitz, John M.; Poling, Bruce E.
       Properties of Gases and Liquids. McGraw-Hill Companies, 1987.