Surface Integral of vector field bounded by two spheres .

Question

enter image description here

A vector field : $$D\mapsto \mathbf{a}_r\frac{\cos^2\phi}{r^3}+\mathbf{a}_\theta\sin\theta$$ Exists in the region between two spherical shells defined by r=1 cm and r=2 cm . Evaluate : 1) Plot the field lines inside the mentioned region

2. the closed surface integral $\oint D\cdot ds$
3. the divergence theorem $\iiint \nabla\cdot DdV$

please help me to solve it , I can solve it using divergence theorem but I don't no how can I solve it by the closed surface integral $\oint D\cdot ds$