Towards petascale simulation of atmospheric circulations with soundproof equations

This paper highlights progress with the development of a petascale implementation of general-purpose high-resolution (nonoscillatory) hydrodynamical simulation code EULAG [Prusa et al. 2008, Comput. Fluids 37, 1193]. The applications addressed are anelastic atmospheric flows in the range of scales from micro to planetary. The new modeldomain decomposition into a three dimensional processor array has been implemented to increase model performance and scalability. The performance of the new code is demonstrated on the IBM BlueGene/L and Cray XT4/XT5 supercomputers. The results show significant improvement of the model efficacy compared to the original decomposition into a two-dimensional processor array in the horizontal - a standard in meteorological models.