function f = floor_force_param ( t, n, un, vn, k, A, Evfloor, cvfloor )
    f = zeros(n,1);
    if (un(k) < 0.00)
        f(k) = (-Evfloor*un(k) -cvfloor*vn(k))*A;
    end
end
