points and the Rankine-Hugoniot conditions at the shock.

Results
=======

![(a) Mach number $M=v_r/a$ on the equatorial plane with $r$ at $t=317.39~t_g$, (b) $M$ on the equatorial plane with $r$, after the solution has reached steady state. (c) Specific angular momentum $\lambda$ on the equatorial plane with $r$, and (d) the distribution of SPH particles in steady state in $r-z$ plane. The injection parameters are injection velocity $v_{inj} = -0.06436$, sound speed $a_{\rm inj} = 0.06328$ and angular momentum $\lambda_{\rm inj} = 1.63$ at $r_{inj}=50.4$. See text for details.](fig1.eps){width="50.00000%"}

![ Time evolution of shock location. Viscosity parameter chosen are (a) $\alpha = 0.0$ (b) $0.005$ and (c) $0.011$, for injection parameters $r_{\rm inj} = 50.4~r_g$ with injection velocity $v_{\rm inj} = -0.06436$, sound speed $a_{\rm inj} = 0.06328$ and angular momentum $\lambda_{\rm inj} = 1.63$. And the panels on the right are plotted for viscosity (d) $\alpha=0$, (e) $0.003$ and (f) $0.007$, for injection parameters $r_{\rm inj} = 50.8~r_g$ with injection velocity $v_{\rm inj} = -0.06532$, sound speed $a_{\rm inj} = 0.06221$ and angular momentum $\lambda_{\rm inj} = 1.7$. Here, the adiabatic index $\gamma = 4/3$.](fig2.eps){width="50.00000%"}

In order to obtain the time dependent axisymmetric, viscous accretion solution, we adopt Smooth 