// Check for windmilling.
double radius = Diameter * 0.375; // 75% of radius
double windmill_cutoff = tan(Pitch * 1.745329E-2) * omega * radius;
- if (Vel > windmill_cutoff) {
- cout << "Windmilling: " << Vel << " > " << windmill_cutoff << endl;
+ if (Vel > windmill_cutoff)
Thrust = -Thrust;
- }
vFn(1) = Thrust;