def METHOD_NAME(
    phase,
    rhos='pore.density.*',
    ks='pore.thermal_conductivity.*',
    MWs='param.molecular_weight.*',
):
    r"""
    Calculates the thermal conductivity of a liquid mixture using the
    correlation in [1]
    Parameters
    ----------
    %(phase)s
    %(rhos)s
    %(ks)s
    %(MWs)s
    Returns
    -------
    References
    ----------
    [1] DIPPR9I
    """
    raise NotImplementedError("This function is not ready yet")
    from chemicals import rho_to_Vm
    xs = phase['pore.mole_fraction']
    kLs = phase.get_comp_vals(ks)
    # kL = numba_vectorized.DIPPR9I(xs, kLs)  # Another one that doesn't work
    Vm = [rho_to_Vm(c.get_comp_vals(rhos), c.get_comp_vals(MWs))
          for c in phase.components.keys()]
    denom = np.sum([xs[i]*Vm[i] for i in range(len(xs))], axis=0)
    phis = np.array([xs[i]*Vm[i] for i in range(len(xs))])/denom
    kij = 2/np.sum([1/kLs[i] for i in range(len(xs))], axis=0)
    kmix = np.zeros_like(xs[0])
    N = len(xs)
    for i in range(N):
        for j in range(N):
            kmix += phis[i]*phis[j]*kij
    return kmix