Method for finding sonic points. Solutions are presented in terms of Mach number M (=v/a) vs log r plot. Θ_pin = 7.162 × 10⁻² for all iterations. Panel a: iterations to obtain inner sonic point r_ci (black circle) and panel b: iterations to obtain outer sonic point r_co (black star). Various branches plotted are multivalued branches (MB; green dashed-dot), transonic (TS; red dashed), and supersonic (SB; blue dashed-dot). Respective Θ_eins are mentioned inside the panels. Panel c: full set of transonic solutions. Global accretion solution (red solid) through r_co and accretion solution through r_ci (red, dashed) is plotted. Equatorial global wind (through r_co) and non-global wind (through r_ci) are represented using red dotted curve. The accretion disc flow parameters used are λ = 2.5, E = 1.000045, Ṁ = 0.001 Ṁ_Edd, and M_BH = 10 M_⊙.