Conservative interpolation between general spherical meshes
An efficient, local, explicit, second-order, conservative interpolation algorithm between spherical meshes is presented. The cells composing the source and target meshes may be either spherical polygons or latitude–longitude quadrilaterals. Second-order accuracy is obtained by piece-wise linear fini...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2017-01-01
|
| Series: | Geoscientific Model Development |
| Online Access: | http://www.geosci-model-dev.net/10/425/2017/gmd-10-425-2017.pdf |
| _version_ | 1818601138887327744 |
|---|---|
| author | E. Kritsikis M. Aechtner Y. Meurdesoif T. Dubos |
| author_facet | E. Kritsikis M. Aechtner Y. Meurdesoif T. Dubos |
| author_sort | E. Kritsikis |
| collection | DOAJ |
| description | An efficient, local, explicit, second-order, conservative interpolation
algorithm between spherical meshes is presented. The cells composing the
source and target meshes may be either spherical polygons or
latitude–longitude quadrilaterals. Second-order accuracy is obtained by
piece-wise linear finite-volume reconstruction over the source mesh. Global
conservation is achieved through the introduction of a <q>supermesh</q>, whose cells
are all possible intersections of source and target cells. Areas and
intersections are computed exactly to yield a geometrically exact method. The
main efficiency bottleneck caused by the construction of the supermesh is
overcome by adopting tree-based data structures and algorithms, from which
the mesh connectivity can also be deduced efficiently.<br><br>The theoretical second-order accuracy is verified using a smooth test
function and pairs of meshes commonly used for atmospheric modelling.
Experiments confirm that the most expensive operations, especially the
supermesh construction, have <i>O</i>(<i>N</i><i>log</i><i>N</i>) computational cost. The method
presented is meant to be incorporated in pre- or post-processing atmospheric
modelling pipelines, or directly into models for flexible input/output. It
could also serve as a basis for conservative coupling between model
components, e.g., atmosphere and ocean. |
| first_indexed | 2024-12-16T12:46:37Z |
| format | Article |
| id | doaj.art-7cd486f7b4f347228d2883bee6dd931e |
| institution | Directory Open Access Journal |
| issn | 1991-959X 1991-9603 |
| language | English |
| last_indexed | 2024-12-16T12:46:37Z |
| publishDate | 2017-01-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj.art-7cd486f7b4f347228d2883bee6dd931e2022-12-21T22:31:17ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032017-01-0110142543110.5194/gmd-10-425-2017Conservative interpolation between general spherical meshesE. Kritsikis0M. Aechtner1Y. Meurdesoif2T. Dubos3Laboratoire d'analyse, géométrie et applications, université Paris 13, 93430 Villetaneuse, FranceLaboratoire de météorologie dynamique, École polytechnique – IPSL, 91128 Palaiseau, FranceLaboratoire des sciences du climat et de l'environnement, CEA – IPSL, 91191 Gif-sur-Yvette, FranceLaboratoire de météorologie dynamique, École polytechnique – IPSL, 91128 Palaiseau, FranceAn efficient, local, explicit, second-order, conservative interpolation algorithm between spherical meshes is presented. The cells composing the source and target meshes may be either spherical polygons or latitude–longitude quadrilaterals. Second-order accuracy is obtained by piece-wise linear finite-volume reconstruction over the source mesh. Global conservation is achieved through the introduction of a <q>supermesh</q>, whose cells are all possible intersections of source and target cells. Areas and intersections are computed exactly to yield a geometrically exact method. The main efficiency bottleneck caused by the construction of the supermesh is overcome by adopting tree-based data structures and algorithms, from which the mesh connectivity can also be deduced efficiently.<br><br>The theoretical second-order accuracy is verified using a smooth test function and pairs of meshes commonly used for atmospheric modelling. Experiments confirm that the most expensive operations, especially the supermesh construction, have <i>O</i>(<i>N</i><i>log</i><i>N</i>) computational cost. The method presented is meant to be incorporated in pre- or post-processing atmospheric modelling pipelines, or directly into models for flexible input/output. It could also serve as a basis for conservative coupling between model components, e.g., atmosphere and ocean.http://www.geosci-model-dev.net/10/425/2017/gmd-10-425-2017.pdf |
| spellingShingle | E. Kritsikis M. Aechtner Y. Meurdesoif T. Dubos Conservative interpolation between general spherical meshes Geoscientific Model Development |
| title | Conservative interpolation between general spherical meshes |
| title_full | Conservative interpolation between general spherical meshes |
| title_fullStr | Conservative interpolation between general spherical meshes |
| title_full_unstemmed | Conservative interpolation between general spherical meshes |
| title_short | Conservative interpolation between general spherical meshes |
| title_sort | conservative interpolation between general spherical meshes |
| url | http://www.geosci-model-dev.net/10/425/2017/gmd-10-425-2017.pdf |
| work_keys_str_mv | AT ekritsikis conservativeinterpolationbetweengeneralsphericalmeshes AT maechtner conservativeinterpolationbetweengeneralsphericalmeshes AT ymeurdesoif conservativeinterpolationbetweengeneralsphericalmeshes AT tdubos conservativeinterpolationbetweengeneralsphericalmeshes |