 Cold radial magnetohydrodynamic wind integral

For a cold equatorial axisymmetric magnetohydrodynamic wind with $B_p=C{R}^{-2}$ in the gravity of a point mass $M$, smooth crossing of the Alfvén surface gives ${\rho }_a={\mu }_0{k}^2$ and $\ell =\omega R_a^2$. Put $x=R/R_a$, $y=\rho /{\rho }_a$, $\alpha ={\omega }^2{R}_a^3/(GM)$ and $\beta =C^2/({\mu }_0{\rho }_{a}GM{R}_a^3)$. The corotating energy, scaled by $GM/R_a$, isThe azimuthal velocity relative to field-line rotation is $\omega R_a(x-x^{-1})/(y-1)$ and the poloidal speed is $C/(\sqrt{{\mu }_0{\rho }_a}{R}_a^2{x}^{2}y)$. The apparent singularity at $x=y=1$ is interpreted through a smooth limiting slope. A nondegenerate magnetosonic crossing on a constant-$f$ curve requires $f_x=f_y=0$ and a real admissible tangent slope.