Pre-calculated turning points of Stability-curve for stationary Shakura-Sunyaev accretion disc

<p>The pre-calculated&nbsp;turning points for&nbsp;surface density <span class="math-tex">\(\Sigma\)</span>, effective temperature&nbsp;<span class="math-tex">\(T_{\rm eff}\)</span>&nbsp;and accretion rate&nbsp;<span class="math-tex">\(\dot{M}\)</span>&nbsp;of S-shaped stability curve for stationary Shakura-Sunyaev viscous accretion disc. Calculations performed&nbsp;for 20&nbsp;linearly scaled values of central source mass&nbsp;from&nbsp;<span class="math-tex">\(1\, M_{\odot}\)</span>&nbsp;to <span class="math-tex">\(20\, M_{\odot}\)</span>, 20&nbsp;logarithmically scaled values of viscous parameter&nbsp;<span class="math-tex">\(\alpha\)</span>&nbsp;from <span class="math-tex">\(3\cdot10^{-4}\)</span>&nbsp;to 0.7, 20&nbsp;logarithmically scaled values of radius&nbsp;from <span class="math-tex">\(7\cdot10^7\)</span>&nbsp;cm to <span class="math-tex">\(5\cdot10^{11}\)</span>&nbsp;cm.</p> <p><a href="https://github.com/AndreyTavleev/DiscVerSt">The new open-source Python3 code</a> is used to calculate the&nbsp;vertical structure of the disc with tabular opacity (solar chemical composition) and both convective and radiative energy transport taking into account. As a result, one can obtain&nbsp;S-shaped stability&nbsp;curve &mdash;&nbsp;a&nbsp;graphically depicted sequence of solutions of the vertical-structure equations, obtained at a single disc radius, in the coordinates&nbsp;of accretion rate or effective temperature versus the surface density. S-curve turn points&nbsp;allow&nbsp;to investigate the stability of disc with respect to temperature and surface density fluctuations.&nbsp;</p> <p>Reading of the data (Python3):</p> <pre><code class="language-python">import numpy as np data = np.genfromtxt('Turn_points_array.dat', names=True)</code></pre> <p>&nbsp;</p>