def massic_rate_sutton(k, conc_evaporating_fluid,
                       wind_speed, pool_area, schmdt_nmbr, r):
    """
    Return the evaporation rate [kg/s m²] of sutton's model
    source : (Fingas, 2015)

    Parameters
    ----------
    k : mass transsfert coefficent [1/m²]
    conc_evaporating_fluid : Concentration of the evaporating fluid [kg/m³]
    wind_speed : Wind speed 10 meters above the surface [m/s]
    pool_area : Area of the pool []
    schmdt_nmbr :  Schmidt number []
    r : empirical (0->2/3) []

    """
    return (k * conc_evaporating_fluid
         * (wind_speed**(7/9)) * (pool_area**(1/9)) * (schmdt_nmbr**r))