 Solution

For a nondegenerate surface take $R>0$, and abbreviate $C=\cosh R$, $H=\sinh R$. Close the upper oblate spheroid by its equatorial disk $D$, obtaining a half-ellipsoid $V$. The divergence isThe integral of $2y$ vanishes by reflection symmetry. At height $z$, the cross-section is a disk of radius $q=C\sqrt{1-z^2/H^2}$. Using planar polar coordinates,ConsequentlyThe outward unit normal vector on $D$ is $-{\mathbf{e}}_z$, so its flux is $-2\pi C^2$. Applying the divergence theorem to the closed surface and subtracting this disk contribution givesThis is a special case of polynomial flux through a half-ellipsoid. Including the closing disk is essential: the constant vertical component has zero divergence but nonzero flux through the open curved surface.