In the method, both electric and magnetic fields are expanded in terms of the vector spherical harmonic (VSH) functions in spherical geometry and the divergencelessness of magnetic field is analytically enforced by a projection method.