function [s,d]=uv2sd(u,v);

if u~=0;
  d=atan(v/u)*180/pi;

% four cases
  if u>0&v>0;
    d=90-d;
  elseif u>0&v<=0;
    d=90+abs(d);
  elseif u<=0&v>0;
    d=270+abs(d);
  elseif u<0&v<0;
    d=270-abs(d);
  end;
else;
  if v<0;
    d=180;
  else;
    d=0;
  end;
end;

s=sqrt(u^2+v^2);

