function [u,v]=rect(umag,theta);
%function [u,v]=rect(umag,theta);
%     theta should be compass degrees
u=umag.*cos( (90-theta)*pi/180 );
v=umag.*sin( (90-theta)*pi/180 );
