The coupling of highly turbulent convection with rotation within a full spherical shell geometry, such as in the solar convection zone, can be studied with the new anelastic spherical harmonic (ASH) code developed to exploit massively parallel architectures. Inter-processor transposes are used to ensure data locality in spectral transforms, a sophisticated load balancing algorithm is implemented and the Legendre transforms, which dominate the workload for large problems, are highly optimized by exploiting the features of cache memory and instruction pipelines. As a result, the ASH code achieves around 120 Mflop/s per node on the Cray T3E and scales nearly linearly for adequately large problem sizes.
['Gary A. Glatzmaier', 'Thomas Clune', 'Jeff Elliott', 'Mark S. Miesch', 'Juri Toomre']
Computational aspects of a code to study rotating turbulent convection in spherical shells