A novel method for finding allowed regions in the space of CFT-data, coined navigator method, was recently proposed in [1]. Its efficacy was demonstrated in the simplest example possible, i.e. that of the mixed-correlator study of the 3D Ising Model. In this paper, we would like to show that the navigator method may also be applied to the study of the family of d-dimensional O(N ) models. We will aim to follow these models in the (d, N ) plane. We will see that the "sailing" from island to island can be understood in the context of the navigator as a parametric optimization problem, and we will exploit this fact to implement a simple and effective path-following algorithm. By sailing with the navigator through the (d, N ) plane, we will provide estimates of the scaling dimensions (∆ φ , ∆ s , ∆ t ) in the entire range (d, N ) ∈ [3, 4] × [1, 3]. We will show that to our level of precision, we cannot see the non-unitary nature of the O(N ) models due to the fractional values of d [2] or N [3] in this range. We will also study the limit N − → 1, and see that we cannot find any solution to the unitary mixed-correlator crossing equations below N = 1.