Calculate the flux of the vector field $\mathbf{F} = \langle x, y, z \rangle$ through the boundary of the cylinder defined by $x^2 + y^2 \le 1$ and $0 \le z \le 1$. Use the Divergence Theorem.