function [alpha] = WingAngle(Vw)
%Returns angle of wing relative to spine (alpha) given the velocity of the water
%relative to the spine (Vw)
%
%Relies on many assumptions...

theta = atan2(-Vw(2,:),-Vw(1,:));
Mag = (Vw(1,:).^2+Vw(1,:).^2).^.5;

alpha = theta;
idx = find(alpha < -30*pi/180)
    alpha(idx) = -30*pi/180;
idx = find(alpha > 10*pi/180)
    alpha(idx) = 10*pi/180;


end

