A 15.2 TFlops Simulation of Geodynamo on the Earth Simulator

For realistic geodynamo simulations, one must solve the magnetohydrodynamic equations to follow time development of thermal convection motion of electrically conducting fluid in a rotating spherical shell. We have developed a new geodynamo simulation code by combining the finite difference method with the recently proposed spherical overset grid called Yin-Yang grid. We achieved performance of 15.2 Tflops (46% of theoretical peak performance) on 4096 processors of the Earth Simulator.