TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids

<p>We present a new method to identify connected components on triangular grids used in atmosphere and climate models to discretize the horizontal dimension. In contrast to structured latitude–longitude grids, triangular grids are unstructured and the neighbors of a grid cell do not simply fol...

Full description

Bibliographic Details
Main Authors: A. Voigt, P. Schwer, N. von Rotberg, N. Knopf
Format: Article
Language:English
Published: Copernicus Publications 2022-10-01
Series:Geoscientific Model Development
Online Access:https://gmd.copernicus.org/articles/15/7489/2022/gmd-15-7489-2022.pdf
_version_ 1811241925286559744
author A. Voigt
P. Schwer
N. von Rotberg
N. Knopf
author_facet A. Voigt
P. Schwer
N. von Rotberg
N. Knopf
author_sort A. Voigt
collection DOAJ
description <p>We present a new method to identify connected components on triangular grids used in atmosphere and climate models to discretize the horizontal dimension. In contrast to structured latitude–longitude grids, triangular grids are unstructured and the neighbors of a grid cell do not simply follow from the grid cell index. This complicates the identification of connected components compared to structured grids. Here, we show that this complication can be addressed by involving the mathematical tool of cubulation, which allows one to map the 2-D cells of the triangular grid onto the vertices of the 3-D cells of a cubical grid. Because the latter is structured, connected components can be readily identified by previously developed software packages for cubical grids. Computing the cubulation can be expensive, but, importantly, needs to be done only once for a given grid. We implement our method in a Python package that we name TriCCo and make available via pypi, gitlab, and zenodo. We document the package and demonstrate its application using simulation output from the ICON atmosphere model. Finally, we characterize its computational performance and compare it to graph-based identifications of connected components using breadth-first search. The latter shows that TriCCo is ready for triangular grids with up to 500 000 cells, but that its speed and memory requirement should be improved for its application to larger grids.</p>
first_indexed 2024-04-12T13:43:28Z
format Article
id doaj.art-f7edd7d7098942cfbb50e890ed2a2f33
institution Directory Open Access Journal
issn 1991-959X
1991-9603
language English
last_indexed 2024-04-12T13:43:28Z
publishDate 2022-10-01
publisher Copernicus Publications
record_format Article
series Geoscientific Model Development
spelling doaj.art-f7edd7d7098942cfbb50e890ed2a2f332022-12-22T03:30:46ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032022-10-01157489750410.5194/gmd-15-7489-2022TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular gridsA. Voigt0P. Schwer1N. von Rotberg2N. Knopf3Department of Meteorology and Geophysics, University of Vienna, Vienna, AustriaInstitute for Algebra and Geometry, Department of Mathematics, Otto-von-Guericke University, Magdeburg, GermanyInstitute for Algebra and Geometry, Department of Mathematics, Otto-von-Guericke University, Magdeburg, GermanyInstitute of Meteorology and Climate Research – Department Troposphere Research, Karlsruhe Institute of Technology, Karlsruhe, Germany<p>We present a new method to identify connected components on triangular grids used in atmosphere and climate models to discretize the horizontal dimension. In contrast to structured latitude–longitude grids, triangular grids are unstructured and the neighbors of a grid cell do not simply follow from the grid cell index. This complicates the identification of connected components compared to structured grids. Here, we show that this complication can be addressed by involving the mathematical tool of cubulation, which allows one to map the 2-D cells of the triangular grid onto the vertices of the 3-D cells of a cubical grid. Because the latter is structured, connected components can be readily identified by previously developed software packages for cubical grids. Computing the cubulation can be expensive, but, importantly, needs to be done only once for a given grid. We implement our method in a Python package that we name TriCCo and make available via pypi, gitlab, and zenodo. We document the package and demonstrate its application using simulation output from the ICON atmosphere model. Finally, we characterize its computational performance and compare it to graph-based identifications of connected components using breadth-first search. The latter shows that TriCCo is ready for triangular grids with up to 500 000 cells, but that its speed and memory requirement should be improved for its application to larger grids.</p>https://gmd.copernicus.org/articles/15/7489/2022/gmd-15-7489-2022.pdf
spellingShingle A. Voigt
P. Schwer
N. von Rotberg
N. Knopf
TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
Geoscientific Model Development
title TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
title_full TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
title_fullStr TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
title_full_unstemmed TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
title_short TriCCo v1.1.0 – a cubulation-based method for computing connected components on triangular grids
title_sort tricco v1 1 0 a cubulation based method for computing connected components on triangular grids
url https://gmd.copernicus.org/articles/15/7489/2022/gmd-15-7489-2022.pdf
work_keys_str_mv AT avoigt triccov110acubulationbasedmethodforcomputingconnectedcomponentsontriangulargrids
AT pschwer triccov110acubulationbasedmethodforcomputingconnectedcomponentsontriangulargrids
AT nvonrotberg triccov110acubulationbasedmethodforcomputingconnectedcomponentsontriangulargrids
AT nknopf triccov110acubulationbasedmethodforcomputingconnectedcomponentsontriangulargrids