MIRA-Datasets: Datasets from Metrics for Intercomparison of Remapping Algorithms

The Metrics for Intercomparison of Remapping Algorithms (MIRA) project provides the Python drivers for the intercomparison study to enable the computation of metrics for different remapping algorithms of interest in ESM.

The dataset repository contains three groups of artifacts: the original test cases used in the study and the output metrics data from four different remapping algorithms, along with some helpful scripts to compare the metrics data. Details are provided below.

All of the input meshes, sampled reference data on the meshes for several uniformly refined resolutions, and regionally refined cases are contained within the Meshes directory.

The uniformly refined meshes for Cubed-Sphere (CS), polygonal quasi-uniform MPAS and Regular Latitude-Longitude (RLL) meshes along with sampled field data for five different fields are provided in Meshes/UniformlyRefined/ directory.

The regionally refined meshes for CS and MPAS meshes around continental-US (CONUS) region with the sampled reference field data is available under Meshes/RegionallyRefined directory.

The input meshes provided under Meshes directory were used to perform a remapping intercomparison study that analyzed the key numerical metrics to gain better insight into the behavior of remapping algorithms, and to compare several key properties under a unified framework. Four different remapping algorithms were considered in this study.

Earth System Modeling Framework (ESMF) Regrid

TempestRemap high-order conservative maps

Generalized Moving-Least-Squares (GMLS) algorithm

A variation with the Clip-And-Assured-Sum (CAAS) algorithm to enforce bounds preservation

Weighted-Least-Squares Essentially Non-oscillatory Remap (WLS-ENOR) scheme

The metrics data collected for each of the cases and remapping algorithms are stored under the MetricsData directory. The metrics CSV files include details about:

Error convergence data in global norms $L_1, L_2, L_{\inf}, H_1$ and $\left|H_1\right|$

Global bounds preservation for determining monotonicity

Local feature preservation through repeated remapping cycles

Grid independence by using test cases with different mesh types and (uniformly refined/regionally refined) resolutions

A set of helpful Python scripts have also been provided to easily compare different aspects of the metrics data to gain more insight into the behavior of the remapping algorithms. These are under Scripts directory.