I have to solve the the following integral $\int_{D}^{} A \cdot dS$, where $D$ is the surface of the hemiellipsoid $(\frac{x}{a})^2+ (\frac{y}{b})^2 + (\frac{z}{c})^2= 1$ a,b and c are known. It is also given from the task that $A$ can be written as $A=(B \cdot r)C$, where $B$ and $C$ are constant vectors, and $r$ is the position vector. From what I see, it seems easiest to use Gauss law for integration, where we instead of integrating over the surface of the hemiellipsoid, we integrate over the volume. However, I don't seem to know what to do with the expression $\nabla \cdot A = \nabla \cdot ((B \cdot r)C)$.
EDIT
I've tried to expand the expression above to $\nabla \cdot A = \nabla \cdot ((B \cdot r)C)= \nabla \cdot (BC \cdot rC)= (\nabla \cdot BC)\cdot rC = 0$, since $BC$ is a constant vector, which should mean that the divergence is 0. But that means that the entire integral will be $0$, which seems odd.
As Ted Shifrin noted, your expression for the divergence of A is incorrect. Resolving it into components, $$r=\begin{bmatrix}x\\y\\z\end{bmatrix}\\A=(B_xx+B_yy+B_zz)\begin{bmatrix}C_x\\C_y\\C_z\end{bmatrix}\\\nabla\cdot{A}=C_xB_x+C_yB_y+C_zB_z=C\cdot{B}$$ Thus, $$\int_DA\cdot{dS}=C\cdot{B}\iiint_VdV$$ Using the formula for the volume of an ellipse and dividing by two, you get: $$\int_DA\cdot{dS}=C\cdot{B}\frac4 6\pi{abc}$$