It is extremely uncommon to be able to predict the velocity profile of a turbulent flow. In twodimensional flows, atmosphere dynamics, and plasma physics, large scale coherent jets are created through inverse energy transfers from small scales to the largest scales of the flow. We prove that in the limits of vanishing energy injection, vanishing friction, and small scale forcing, the velocity profile of a jet obeys an equation independent of the details of the forcing. We find another general relation for the maximal curvature of a jet and we give strong arguments to support the existence of an hydrodynamic instability at the point with minimal jet velocity. Those results are the first computations of Reynolds stresses and self consistent velocity profiles from the turbulent dynamics, and the first consistent analytic theory of zonal jets in barotropic turbulence.Theoretical prediction of velocity profiles of inhomogeneous turbulent flows is a long standing challenge, since the nineteenth century. It involves closing hierarchy for the velocity moments, and for instance obtaining a relation between the Reynolds stress and the velocity profile. Since Boussinesq in the nineteenth century, most of the approaches so far have been either empirical or phenomenological. Even for the simple case of a three dimensional turbulent boundary layer, plausible but so far unjustified similarity arguments may be used to derive von Kármán logarithmic law for the turbulent boundary layer (see for instance [1]), but the related von Kármán constant [2] has never been computed theoretically. Still this problem is a crucial one and has some implications in most of scientific fields, in physics, astrophysics, climate dynamics, and engineering. Equations (6-7), (9), and (11) are probably the first prediction of the velocity profile for turbulent flows, and relevant for barotropic flows.In this paper we find a way to close the hierarchy of the velocity moments, for the equation of barotropic flows with or without effect of the Coriolis force. This two dimensional model is relevant for laboratory experiments of fluid turbulence [3], liquid metals [4], plasma [5], and is a key toy model for understanding planetary jet formation [6] and basics aspects of plasma dynamics on Tokamaks in relation with drift waves and zonal flow formation [7]. It is also a relevant model for Jupiter troposphere organization [8]. Moreover, our approach should have future implications for more complex turbulent boundary layers, which are crucial in climate dynamics in order to quantify momentum and energy transfers between the atmosphere and the ocean.It has been realized since the sixties and seventies in the atmosphere dynamics and plasma communities that in some regimes two dimensional turbulent flows are * Eric.Woillez@ens-lyon.fr † Freddy.Bouchet@ens-lyon.fr strongly dominated by large scale coherent structures. Jets and large vortices are often observed in numerical simulations or in experiments, but the general mechanism leading to such an organization of t...