A Parallel Ocean Model for High Resolution Studies

Domain Decomposition Methods is used for the parallelization of the LODYC ocean general circulation model. The local dependencies problem is solved by using a pencil splitting and an overlapping strategy. Two different parallel solvers of the surface pressure gradient, a preconditioned conjugate gradient method and a Dual Schur Complement method, have been implemented. The code is now used for the high resolution study of the Atlantic circulation by the CLIPPER research project and is one of the components of the operational oceanography MERCATOR project.