Handling variable length input/output arrays, terminating computation when a condition is met (Boundary layer calculations)

Question

I\'m trying to write code to calculate boundary layer development over an aerofoil.

I start with a set of points, S, that describe the shape of the aerof