Use the Divergence Theorem to evaluate $\iint_{s} \mathbf{F} \cdot d \mathbf{S}$, where
$$
\mathbf{F}(x, y, z)=z^{2} x \mathbf{i}+\left(\frac{1}{3} y^{3}+\tan ^{-1} z\right) \mathbf{j}+\left(x^{2} z+y^{2}\right) \mathbf{k}
$$
and $S$ is the top half of the sphere $x^{2}+y^{2}+z^{2}=1$. [Hint: Note that $S$ is not a closed surface. First compute integrals over $S_{1}$ and $S_{2}$, where $S_{1}$ is the disk $x^{2}+y^{2} \leq 1$, oriented downward, and $S_{2}=S \cup S_{1 .}$ ]