A Solovev-like solution describing equilibria with field aligned incompressible flows [G. N. Throumoulopoulos and H. Tasso, Phys. Plasmas 19, 014504 (2012)] is extended to non parallel flows. The solution expressed as a superposition of Bessel functions contains an arbitrary number of free parameters which are exploited to construct a variety of configurations including ITER shaped ones. For parallel flows, application of a sufficient condition for linear stability shows that this condition is satisfied in an appreciable part of the plasma region on the high-field side mostly due to the variation of the magnetic field perpendicular to the magnetic surfaces. Also, the results indicate that depending on the shape of the Mach-function profile and the values of the free parameters the flow and flow shear may have either stabilizing or destabilizing effects.