We take the outer boundary \(R_{\rm out}\)  = 100.
Referring to the hydrodynamical results in our previous paper [1]},
we adopt \(\lambda \) = 1.5, \(v_{\rm out}\) = -0.072241 and \(a_{\rm out}\) = 0.037566
at the outer radial boundary.
Then, giving the input density \(\rho _{\rm out}\) , we get other variables at the radial outer boundary.
We examine six models of adiabatic case with \(\rho _{\rm out} =2 \)  (AD), optically thin case with \(\rho _{\rm out}\)  =2 (Thin1) and four optically thick cases with \(\rho _{\rm out}\)  = 20,
200, 2\(\times 10^3\)  and 2\(\times 10^4\)  (Thick2, Thick3, Thick4, and Thick5), respectively.
Hereafter the density is shown in unit of \({\rho }_0 =10^{-8}\)  g cm\(^{-3}\) .
In model Thin1, the input gas is given to be optically thin but in models
Thick2 to Thick5 the input accreting gas
is optically thick. We should notice that model Thin1 has the same
parameters as the adiabatic case AD but in the latter case, the radiation transport is not treated.
The flow parameters and other input variables at the outer radial boundary are listed in Table 1,
where \(\mathcal {E}\) = 2.328 \(\times 10^{-5}\)  and \(\lambda \) = 1.5 are used.
<TABLE>