Total unbound mass flux as a function of time, normalized to the total injected mass flux. The fluxes are calculated at the simulation domain up to r = 40.0 (i.e., integrated over all of the surface of the volume (27, 40) × (0.42π, 0.58π) × (0.81π, 0.96π)). We also calculate the mass flux at r = 45.0, 50.0 and do not find a significant difference in the quasi-steady state.