Additionally, [1]} prescribe stellar wind flow
assigned at the wind boundary \(r_{\rm in}\)  to account for the asymmetry
function as,
\(v(r_{\rm in},\theta ) = \frac{ v_{\rm w} }{ \sqrt{ f(\theta ) } },\) 
