Abstract. I present a new family of analytical flow solutions to the incompressible Stokes equation in a spherical shell. The velocity is tangential to both inner and outer boundaries, the viscosity is radial and of the power-law type, and the solution has been designed so that the expressions for velocity, pressure, and body force are simple polynomials and therefore simple to implement in (geodynamics) codes. Various flow average values, e.g., the root mean square velocity, are analytically computed. This forms the basis of a numerical benchmark for convection codes and I have implemented it in two finite-element codes: ASPECT and ELEFANT. I report error convergence rates for velocity and pressure.