def METHOD_NAME(
    phase,
    T="pore.temperature",
    salinity="pore.salinity",
):
    r"""
    Calculates thermal conductivity of pure water or seawater at atmospheric
    pressure using the correlation given [1].
    Values at temperature higher the normal boiling temperature are calculated
    at the saturation pressure.
    Parameters
    ----------
    %(phase)s
    %(T)s
    %(salinity)s
    Returns
    -------
    value : ndarray
        A numpy ndarray containing thermal conductivity of water/seawater in [W/m.K]
    Notes
    -----
    T must be in K, and S in g of salt per kg of phase, or ppt (parts per
    thousand). The correlation is valid for 273 < T < 453 K and
    0 < S < 160 g/kg within 3 percent accuracy.
    References
    ----------
    [1] D. T. Jamieson, and J. S. Tudhope, Desalination, 8, 393-401, 1970.
    """
    T = phase[T]
    if salinity in phase.keys():
        S = phase[salinity]
    else:
        S = 0
    T68 = 1.00024 * T  # convert from T_90 to T_68
    SP = S / 1.00472  # convert from S to S_P
    k_sw = 0.001 * (
        10
        ** (
            np.log10(240 + 0.0002 * SP)
            + 0.434
            * (2.3 - (343.5 + 0.037 * SP) / T68)
            * ((1 - T68 / (647.3 + 0.03 * SP))) ** (1 / 3)
        )
    )
    value = k_sw
    return value