function power=wind2power_stormcontrol(wind)
% Enercon 44
% https://pdf.directindustry.com/pdf/enercon/enercon-product-overview/20877-243513.html

% unit: kW
u = [0 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 21 25 25.6 26.2 26.8 27.4 28 28.6 29.2 29.8 30.4 31 31.6 32.2 32.8 33.4 34];
P = [0 0 4 20 50 96 156 238 340 466 600 710 790 850 880 905 910 910 910 905 880 850 790 710 600 466 340 238 156 96 50 20 4 0];
power = nan(size(wind));

for i = 1:length(u)-1
	index = find(wind>u(i) & wind<=u(i+1));
    power(index) = P(i) + (P(i+1)-P(i))*(wind(index)-u(i))/(u(i+1)-u(i));
end
index = find(wind==0|wind>u(end));
power(index)=0;