Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 10, issue 11 | Copyright
Geosci. Model Dev., 10, 4175-4186, 2017
© Author(s) 2017. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model description paper 17 Nov 2017

Model description paper | 17 Nov 2017

Parcels v0.9: prototyping a Lagrangian ocean analysis framework for the petascale age

Michael Lange1 and Erik van Sebille2,3 Michael Lange and Erik van Sebille
  • 1Grantham Institute & Department of Earth Science and Engineering, Imperial College London, London, UK
  • 2Institute for Marine and Atmospheric research Utrecht, Utrecht University, Utrecht, the Netherlands
  • 3Grantham Institute & Department of Physics, Imperial College London, London, UK

Abstract. As ocean general circulation models (OGCMs) move into the petascale age, where the output of single simulations exceeds petabytes of storage space, tools to analyse the output of these models will need to scale up too. Lagrangian ocean analysis, where virtual particles are tracked through hydrodynamic fields, is an increasingly popular way to analyse OGCM output, by mapping pathways and connectivity of biotic and abiotic particulates. However, the current software stack of Lagrangian ocean analysis codes is not dynamic enough to cope with the increasing complexity, scale and need for customization of use-cases. Furthermore, most community codes are developed for stand-alone use, making it a nontrivial task to integrate virtual particles at runtime of the OGCM. Here, we introduce the new Parcels code, which was designed from the ground up to be sufficiently scalable to cope with petascale computing. We highlight its API design that combines flexibility and customization with the ability to optimize for HPC workflows, following the paradigm of domain-specific languages. Parcels is primarily written in Python, utilizing the wide range of tools available in the scientific Python ecosystem, while generating low-level C code and using just-in-time compilation for performance-critical computation. We show a worked-out example of its API, and validate the accuracy of the code against seven idealized test cases. This version 0.9 of Parcels is focused on laying out the API, with future work concentrating on support for curvilinear grids, optimization, efficiency and at-runtime coupling with OGCMs.

Publications Copernicus
Short summary
Here, we present version 0.9 of Parcels (Probably A Really Computationally Efficient Lagrangian Simulator). Parcels is an experimental prototype code aimed at exploring novel approaches for Lagrangian tracking of virtual ocean particles in the petascale age. The modularity, flexibility and scalability will allow the code to be used to track water, nutrients, microbes, plankton, plastic and even fish.
Here, we present version 0.9 of Parcels (Probably A Really Computationally Efficient Lagrangian...