We solved analytically viscous two-dimensional (2D) fluid equations for accretion and outflows in spherical polar coordinates (r, θ, φ) and obtained explicitly flow variables in r− and θ−directions around black holes (BHs). We investigated global transonic advection-dominated accretion flow (ADAF) solutions in r−direction on an equatorial plane with using Paczyński-Wiita potential. We used radial flow variables of ADAFs with symmetric conditions on the equatorial plane, as initial values for integration in θ−direction. In the study of 2D disk structure, we used twoazimuthal components of viscous stress tensors namely, τ rφ and τ θφ . Interestingly, we found that the whole advective disk is not participating in outflow generation and the outflows form close to the BHs. Normally, outflow strength increased with increasing viscosity parameter (α 1 ), mass-loss parameter (s) and decreasing gas pressure ratio (β). Outflow region increased with increasing s, α 1 for τ rφ and decreasing α 2 for τ θφ . The τ θφ is effective in angular momentum transportation at high latitude and outflows collimation along an axis of symmetry since it changes polar velocity (v θ ) of the flow. The outflow emission is also affected by the ADAF size and decreased with decreasing it. Transonic surfaces formed for both inflows (v r < 0, very close to BH) and outflows (v r > 0). We also explored no outflows, outflows and failed outflows regions, which mainly depend on the viscosity parameters.