function [cosu,sinu]=sincos(umag,udir)
%function [cosu,sinu]=sincos(umag,udir)
dummy=exp(log(umag)+sqrt(-1)*udir*pi/180);
cosu=real(dummy);
sinu=imag(dummy);
