def METHOD_NAME( F, x, y, z, t, amplitude, length_scale ):
    """
    Function to be called at each timestep on the particles
    """
    return( F + amplitude*math.cos( 2*np.pi*(z-c*t)/length_scale ) )