def METHOD_NAME(fluxes, noise_level=0.1):
    '''
    Adds random noise to a light curve.
    Parameters
    ----------
    noise_level: float
        The standard deviation of the normal distribution used to sample the added noise.
    Returns
    -------
    fluxes+noise 
    '''
    noise = np.random.normal(loc=0, scale=noise_level, size = fluxes.shape[0])
    return fluxes + noise, noise_level*np.ones_like(fluxes)