Insertion of PETSc in the NEMO stack software driving NEMO towards exascale computing

This paper addresses the scientific challenges related to high level implementation strategies which steer the NEMO (Nucleus for European Modelling of the Ocean) code toward the effective exploitation of the opportunities offered by exascale systems. We consider, as case studies, two components of the NEMO ocean model (OPA-Ocean PArallelization): the Sea Surface Height equation solver and the Variational Data Assimilation module. The advantages rising from the insertion of consolidated scientific libraries in the NEMO code are highlighted: such advantages concern both the “software quality” improvement (see the software quality parameters like robustness, portability, resilience, etc.) and the reduction of time spent for software development and maintenance. Finally, we consider the Shallow Water equations as a toy model for NEMO ocean model to show how the use of PETSc objects predisposes the application to gain a good level of scalability and efficiency when the most suitable level of abstraction is used.

[1]  Almerico Murli,et al.  Monitoring and Migration of a PETSc-based Parallel Application for Medical Imaging in a Grid computing PSE , 2006, Grid-Based Problem Solving Environments.

[2]  Jack Dongarra,et al.  High-performance computing systems: Status and outlook* , 2012, Acta Numerica.

[3]  G. Madec.  NEMO ocean engine , 2008 .

[4]  Rossella Arcucci,et al.  Data Assimilation Achievements in HPC sysyems: experiments on OceanVar in the Mediterranean Sea , 2013 .

[5]  Almerico Murli,et al.  HPC computation issues of the incremental 3D variational data assimilation scheme in OceanVar software , 2012 .

[6]  Kristian Mogensen,et al.  The NEMOVAR ocean data assimilation system as implemented in the ECMWF ocean analysis for System 4 , 2012 .

[7]  Luisa Carracciuolo,et al.  HADAB: Enabling Fault Tolerance in Parallel Applications Running in Distributed Environments , 2011, PPAM.

[8]  Barry F. Smith,et al.  PETSc ’ s Software Strategy for the Design Space of Composable Extreme-Scale Solvers , 2012 .

[9]  Pier Luca Maffettone,et al.  Computational simulations of 3D large-scale time-dependent viscoelastic flows in high performance co , 2011 .

[10]  N. Pinardi,et al.  An oceanographic three-dimensional variational data assimilation scheme , 2008 .

[11]  Italo Epicoco,et al.  The performance model for a parallel SOR algorithm using the red-black scheme , 2012, Int. J. High Perform. Syst. Archit..

[12]  Almerico Murli,et al.  A Scalable Approach for Variational Data Assimilation , 2014, J. Sci. Comput..

[13]  Almerico Murli,et al.  OceanVAR Software for Use with NEMO: Documentation and Test Guide , 2012 .

[14]  Margaret H. Wright,et al.  The opportunities and challenges of exascale computing , 2010 .

[15]  Ecmwf Newsletter,et al.  EUROPEAN CENTRE FOR MEDIUM-RANGE WEATHER FORECASTS , 2004 .

[16]  John Shalf,et al.  The International Exascale Software Project roadmap , 2011, Int. J. High Perform. Comput. Appl..

[17]  Jed Brown,et al.  The Portable Extensible Toolkit for Scientific computing , 2012 .

[18]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .