I need to compute the integral of Integ(J^2 rdrdz) over an rz block (in axisymmetric geometry). Is there a convenient way to do this?