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

Erik Bernsen*, Henk A. Dijkstra, Jonas Thies, Fred W. Wubs

*Corresponding author for this work

    Research output: Contribution to journalArticleAcademicpeer-review

    15 Citations (Scopus)
    445 Downloads (Pure)


    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. (C) 2010 Elsevier Inc. All rights reserved.

    Original languageEnglish
    Pages (from-to)8167-8179
    Number of pages13
    JournalJournal of computational physics
    Issue number21
    Publication statusPublished - 20-Oct-2010


    • Ocean modelling
    • Jacobian-free Newton-Krylov methods
    • Preconditioning
    • SOLVER

    Cite this