To use the divergence theorem you need to close the surface with the plane that is gives. The flux though the plane is easy to calculate(you don't need an integral for a constant vector field though a flat surface). So find the flux with the divergence theorem then subtract off the plane to get the flux of the parabolic part.

P.S you limits of integration are incorrect on the triple integral z goes from the plane to the parabolic cylinder.