function sd = stdVel(S,a,b)
%STD velocity
sd=a+b*abs(S);
end

