by Warnecke et al. (2016) suggest that radial pumping is downward throughout the convection zone below $45^{\circ}$ and upward above $45^{\circ}$, while latitudinal pumping is poleward at the surface and equatorward at the base of the convection zone. Our generated turbulent pumping profiles in the northern hemisphere (defined within $0 \leq \theta \leq \pi/2$) based on the suggestions of Warnecke *et al.* (2006) are: $$\begin{aligned}
\gamma_r = -  \gamma_{0r} \left[ 1 + \rm{erf}\left( \frac{r - 0.715R_\odot}{0.05R_\odot}\right) \right] \left[ 1 - \rm{erf} \left( \frac{r-0.98R_\odot}{0.08R_\odot}\right) \right] \nonumber \\
\times \sin(4 \theta), ~~~~\end{aligned}$$ $$\begin{aligned}
\gamma_\theta = \left\{\begin{array}{cc}
\gamma_{0\theta} \sin \left[\frac{2 \pi (r- R_p)}{R_0-R_p}\right]
\times \cos \theta \sin^4 \theta ~~~~  & r \geq R_p\\
  0  & r<R_p
\end{array}\right.,\end{aligned}$$

where $R_p= 0.76R_\odot$ i.e. the penetration depth of the latitudinal pumping. The amplitudes of $\gamma_r$ and $\gamma_\theta$ are determined by the value of $\gamma_{0r}$ and $\gamma_{0\theta}$ respectively. Turbulent pumping profiles in the southern hemisphere are generated by replacing colatitude $\theta$ by $(\pi- \theta)$. Fig. 8 (the top and middle panels) show that our generated turbulent pumping profiles capture the basic essence of the suggestions made by Warnecke *et al.* (2016). Our simulations (with shallow meridional flow) and the more complex turbulent pumping profile gleaned from Warnecke 