since and .
After making the required substitutions your flux integral becomes , where is the region of the xy-plane defined by . I'd suggest switching to polar coordinates to get this double integral.
Note that the region is found by substituting z = 1 into .
Of course, you could always switch to cylindrical coordinates at the very start of the problem .......