A high-degree (degree l = 6 and order m = 0, 1, 2, ..., l. High-order model for short) and steady thermal free convective motion of an infinite Prandtl number and Boussinesq fluid in a spherical shell is calculated by a Galerkin method. Convection is driven by an imposed temperature drop across top rigid and bottom stress-free isothermal boundaries only heated from below of the shell. In this paper, the scalar poloidal and fluctuating temperature fields are expanded into associated Legendre polynomials with degree l= 6 and order m = 0, 1, 2, ..., l. Compared with zero-order model (degree 1= 6 and order m= 0), from which 2-D longitudinal (r-O) profiles can be obtained, high-order model can provide a series of southerly (r-0), easterly (r-cp) and radial (0-f) velocity profiles, which probably reveal more detail features of mass motion in the mantle. It is found that Rayleigh number has great effects on the patterns and velocities of thermal free convection and controls the relative ratio of hot and cold plume in the shell. Probably, the present results mainly reveal the mass motion in the lower mantle, while the striking differences of convection patterns from velocities at different positions have important geodynamical significances.