A　spectral　method　is　used　to　study　numerically　flow　between　two　concentric　rotating　spheres.　Analytical　expressions　of　eigenfunctions　of　a　Stokes　operator　in　the　spherical　gap　region　are　presented,　numerical　experiments　are　carried　out　by　applying　these　expressions　to　the　spectral　approximation　of　spherical　Couette　flow,　and　numerical　results　are　given.　(C)　2004　Elsevier　Ltd.　All　rights　reserved.