The flow equilibrium problem with prescribed adiabatic constraints may be solved by simultaneous calculations of flux surface geometry and original profile functions. A numerical technique, alternative to Grad's well-known ADM method has been proposed to deal with the slow adiabatic evolution of a toroidal plasma with flows. In the case of field-aligned sub-Alfvénic flow the system is in the second elliptic regime if β < A2/(1 – A2) at the magnetic axis, where A is the Mach Alfvén number of the flow.