def compute_Rn_c(albedo_c, T_air, T_c, T_s, e_atm, Rs_c, F):
    """Compute Canopy Net Radiation

    Parameters
    ----------
    albedo_c : ee.Image
    T_air : ee.Image
        Air temperature (Kelvin).
    T_c : ee.Image
        Canopy temperature (Kelvin).
    T_s : ee.Image
        Soil temperature (Kelvin).
    e_atm : ee.Image
    Rs_c : ee.Image
    F : ee.Image

    Returns
    -------
    Rn_c : ee.Image

    """
    kL = 0.95
    eps_s = 0.94
    eps_c = 0.99
    Lc = T_c.expression('eps_c * 5.67E-8 * (T_c ** 4)', {'eps_c': eps_c,
        'T_c': T_c})
    Ls = T_s.expression('eps_s * 5.67E-8 * (T_s ** 4)', {'eps_s': eps_s,
        'T_s': T_s})
    Rle = T_air.expression('e_atm * 5.67E-8 * (T_air ** 4)', {'e_atm':
        e_atm, 'T_air': T_air})
    Rn_c = albedo_c.expression(
        '((1 - albedo_c) * Rs_c) + ((1 - exp(-kL * F)) * (Rle + Ls - 2 * Lc))',
        {'albedo_c': albedo_c, 'F': F, 'kL': kL, 'Lc': Lc, 'Ls': Ls, 'Rle':
        Rle, 'Rs_c': Rs_c})
    return Rn_c
