A Scalable and Adaptable Solution Framework within Components of the Community Climate System Model

  • Authors:

  • Katherine J. Evans;Damian W. Rouson;Andrew G. Salinger;Mark A. Taylor;Wilbert Weijer;James B. White, Iii

  • Affiliations:

  • Oak Ridge National Laboratory, Oak Ridge, TN 37831;Sandia National Laboratory, Livermore, CA 94551;Sandia National Laboratory, Albuquerque, NM 87185;Sandia National Laboratory, Albuquerque, NM 87185;Los Alamos National Laboratory, Los Alamos, NM 87545;Oak Ridge National Laboratory, Oak Ridge, TN 37831

  • Venue:
  • ICCS 2009 Proceedings of the 9th International Conference on Computational Science
  • Year:
  • 2009

Quantified Score

Hi-index 0.01

Visualization

Abstract

A framework for a fully implicit solution method is implemented into (1) the High Order Methods Modeling Environment (HOMME), which is a spectral element dynamical core option in the Community Atmosphere Model (CAM), and (2) the Parallel Ocean Program (POP) model of the global ocean. Both of these models are components of the Community Climate System Model (CCSM). HOMME is a development version of CAM and provides a scalable alternative when run with an explicit time integrator. However, it suffers the typical time step size limit to maintain stability. POP uses a time-split semi-implicit time integrator that allows larger time steps but less accuracy when used with scale interacting physics. A fully implicit solution framework allows larger time step sizes and additional climate analysis capability such as model steady state and spin-up efficiency gains without a loss in scalability. This framework is implemented into HOMME and POP using a new Fortran interface to the Trilinos solver library, ForTrilinos, which leverages several new capabilities in the current Fortran standard to maximize robustness and speed. The ForTrilinos solution template was also designed for interchangeability; other solution methods and capability improvements can be more easily implemented into the models as they are developed without severely interacting with the code structure. The utility of this approach is illustrated with a test case for each of the climate component models.