Show simple item record

dc.contributor.authorNordam, Tor
dc.contributor.authorDuran, Rodrigo
dc.identifier.citationGeoscientific Model Development. 2020, 13 (12), 5935-5957.en_US
dc.description.abstractA common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field pre-calculated with an ocean model. Mathematically, this is simply numerical integration of an Ordinary Differential Equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge-Kutta integrators of orders 1-4, as well as three variable-step Runge-Kutta methods (Bogacki-Shampine 3(2), Dormand-Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, that are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand-Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods, at a reduced computational cost. The best results are seen for coarser resolutions (4 km and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.en_US
dc.publisherCopernicus Publications on behalf of the European Geosciences Union.en_US
dc.rightsNavngivelse 4.0 Internasjonal*
dc.titleNumerical integrators for Lagrangian oceanographyen_US
dc.typeJournal articleen_US
dc.typePeer revieweden_US
dc.rights.holder© Author(s) 2020en_US
dc.subject.nsiVDP::Matematikk og naturvitenskap: 400en_US
dc.subject.nsiVDP::Mathematics and natural scienses: 400en_US
dc.source.journalGeoscientific Model Developmenten_US
dc.relation.projectNorges forskningsråd: 262741en_US
dc.relation.projectNorges forskningsråd: 267793en_US

Files in this item


This item appears in the following Collection(s)

Show simple item record

Navngivelse 4.0 Internasjonal
Except where otherwise noted, this item's license is described as Navngivelse 4.0 Internasjonal