\bigg( 1 -  \frac{c(p) \, \rho_0^2 \, q_0^2}{A^2}   \bigg )  \frac{dp}{dl} = \rho(p) \, g \, \frac{dz}{dl}  - \frac{\rho_0^2 \, q_0^2 }{2 A^2 d} \frac{f(p)}{\rho(p)}

u(l) = \frac{\rho_0 \cdot q_0}{\rho(p) \cdot A(l)}

q(l) = \frac{\rho_0 \cdot q_0}{\rho(p)}

p(l) = p_0 + \rho \, g \, z(l) - \frac{\rho_0 \, q_0^2 }{2 A^2 d} \, f_0 \, l

\frac{dp}{dl} = \rho \, g \cos \theta(l) - \frac{\rho_0 \, q_0^2 }{2 A^2 d} \, f_0 

u(l) = \frac{q_0}{A(l)}

q(l) =q_0 = \rm const