import math


def volatilization_rate(k, mass, dz, dt):
    """
    Compute the transfer of the volatilization (kg/s)
    source : (“CHEMMAP technical User’s manual 6.10,” 2014)

    params
    ------
    k : volatilization coefficient [m/s]
    dz : vertical diffusivity [m²/s]
    dt : length of a timestep [s]
    mass: mass of the pollutant [kg]
    """

    return k * mass / math.sqrt(2 * dz * dt)