Density at the maximum speed of an isothermal breeze ($r_M=GM/(2a^2),{\rho }_M=\dot{M}{a}^3/(\pi G^2{M}^2{\mathcal{M}}_M)$)

A steady spherical subsonic isothermal breeze with isothermal sound speed $a$ obeys $(u-a^2/u)u^{\prime }=2a^2/r-GM/r^2$. Since $0<u<a$, its finite-radius speed maximum lies at $r_M=GM/(2a^2)$. Spherical mass conservation then gives the displayed density, where ${\mathcal{M}}_M=u(r_M)/a\in (0,1)$ is the peak Mach number. This free parameter is selected by a boundary condition, not by $T$, $M$ and $\dot{M}$ alone. Therefore ${\rho }_M>\dot{M}{a}^3/(\pi G^2{M}^2)$, while a peak Mach number of order unity gives that expression as a characteristic scale.

A direct demonstration of the free parameter comes from the first integral $u^2/2-a^2\ln (u/a)=2a^2\ln (r/r_M)+GM/r+C$. Its additive constant remains free for a breeze; imposing the regular transonic branch would be a different condition.