The application code under investigation simulates dynamos in the Earth's core. It uses pseudo-spectral methods and a discretization of the physical properties of the core to solve the Boussinesq equations of momentum, heat diffusion, and magnetic induction for the fluid velocity, codensity, and magnetic field in a spherical shell.
The most costly computational component of the simulation involves the repeated calculation of discrete Legendre transforms.
The resolution of the discretization can be changed by modifying the number of grid lines in the radial dimension of the spherical shell, and/or number of divisions in the angular dimensions (theta and phi).
For more background on the details of the simulation code, the numerical methods involved, and its parallel implementation, see the supplementary material in the Appendix.
