The application of Jacobian-free Newton-Krylov methods to reduce the spin-up time of ocean general circulation models

  • Authors:

  • Erik Bernsen;Henk A. Dijkstra;Jonas Thies;Fred W. Wubs

  • Affiliations:

  • Institute for Marine and Atmospheric research Utrecht, Utrecht University, The Netherlands;Institute for Marine and Atmospheric research Utrecht, Utrecht University, The Netherlands;Johann Bernoulli Institute for Mathematics and Computer Science, Groningen University, The Netherlands;Johann Bernoulli Institute for Mathematics and Computer Science, Groningen University, The Netherlands

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2010

Quantified Score

Hi-index 31.45

Visualization

Abstract

In present-day forward time stepping ocean-climate models, capturing both the wind-driven and thermohaline components, a substantial amount of CPU time is needed in a so-called spin-up simulation to determine an equilibrium solution. In this paper, we present methodology based on Jacobian-Free Newton-Krylov methods to reduce the computational time for such a spin-up problem. We apply the method to an idealized configuration of a state-of-the-art ocean model, the Modular Ocean Model version 4 (MOM4). It is shown that a typical speed-up of a factor 10-25 with respect to the original MOM4 code can be achieved and that this speed-up increases with increasing horizontal resolution.