A data porting tool for coupling models with different discretization needs

The presented work is part of a larger research program dealing with developing tools for coupling biogeochemical models in contaminated landscapes. The specific objective of this article is to provide researchers with a data porting tool to build hexagonal raster using information from a rectangular raster data (e.g. GIS format). This tool involves a computational algorithm and an open source software (written in C). The method of extending the reticulated functions defined on 2D networks is an essential key of this algorithm and can also be used for other purposes than data porting. The algorithm allows one to build the hexagonal raster with a cell size independent from the geometry of the rectangular raster. The extended function is a bi-cubic spline which can exactly reconstruct polynomials up to degree three in each variable. We validate the method by analyzing errors in some theoretical case studies followed by other studies with real terrain elevation data. We also introduce and briefly present an iterative water routing method and use it for validation on a case with concrete terrain data. • We provide a method to resample piecewise constant functions. • The method recovers cubic polynomials and doesn't introduce spurious oscillations. • We construct a hexagonal cellular automaton model of flow routing. • The water flows from a cell to the neighbors pointed by the water velocity. • The velocity aligns to the projection of gravitational force on the water surface.