The aim is to improve exascale readiness for the Ocean model FESOM2. FESOM2 implements unstructured meshes with variable resolution, employing a finite volume method. The mesh flexibility allows to increase resolution in dynamically active regions, while keep a relatively coarse-resolution setup elsewhere.
While performing a simulation, FESOM2 has to take into account vast datasets of atmospheric forcings, stored as files in netCDF format. By developing an independent, unit tested module to distribute the load on the HPC system when reading and synchronizing the forcing data, the impact of huge forcing data sets on the simulation is minimized and workload is more equally spread on the file system and the HPC network. In addition, actual file reading could be offloaded and done asynchronously by using a dedicated C++ library. This results in a speed improvement of up to factor 11 when reading big ERA5 reanalysis forcing data sets.
Starting from a stand-alone version of FESIM, the Finite-Element Sea Ice Model used in FESOM2, work has started to compare the performance of the currently used explicit pseudo-time stepping (mEVP model, Link: ) with the VP model. The latter, combined with provenly scalable Schwarz-type preconditioning is both significantly faster and reduces the number of communication steps during the simulations.