the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Centroids in second-order conservative remapping schemes on spherical coordinates
Abstract. The transformation of data from one grid system to another is common in climate studies. Among the many schemes used for such transformations is second-order conservative remapping. In particular, a second-order conservative remapping scheme first introduced in 1987 and extended in 1999 to work on the general grids of a sphere has, either directly or indirectly, has served as an important base in a variety of studies.
In this study, the author describes a fundamental problem in the derivation of the method proposed by a pioneer study relating to the treatment of the centroid used as a reference point for the second-order terms in the longitudinal direction. In principle, use of the original formulation may cause damage to the entire remapping result. However, a method's native implementation software includes a preprocessing procedure that tends to minimize or even erase the error as a side effect in many, if not most, typical applications. In this study, three alternative formulations are proposed and tested and are shown to work in a simple application.
- Preprint
(1207 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 02 Jul 2024)
-
RC1: 'Comment on egusphere-2024-1101', Moritz Hanke, 24 Apr 2024
reply
Fuyuki Saito found an error in the formulation of the weights for second order conservative remapping paper by P. W. Jones from 1999 (referred to as J99). This error also made it into the SCRIP library, which is based on J99. This software library has been the basis for conservative interpolation in many climate models for many years, which makes the finding presented here a substantive contribution to the community. Unfortunately my limited understanding of the underlying math prevents me from commenting on the correctness of the presented formulas.
The author reproduces the results from J99 and describes in detail the error made therein. This is followed by an analysis of the SCRIP source code and the description of preprocessing step in SCRIP. This step avoids the error having an actual impact on the interpolation results for many common use-cases. Multiple solutions to problem are derived and explained in-depth. Afterwards the impact of error is analyzed in great detail and compared to presented solutions.
Overall the paper is very well written and after a minor revisit, I would recommend it for publishing.
General remarks:
Even through J99 and the SCRIP library have been widely used in the past, the author seems to wrongfully assume that this is still the case for current software (e.g. ESMF). This is apparent in the additional remarks and the summary of the paper. However, ESMF and other software (e.g. YAC or XIOS) use implementations for first and second order conservative remapping, which significantly differ from the SCRIP library and are therefore not prone to the presented error. Second order conservative remapping in these implementations is based on Kritsikis et al., 2017 and not J99.
In the summary of the paper the author encourages the further use of the SCRIP library. However, he fails to mention the various other drawbacks of it, which would lead me to a different conclusion. These drawbacks include inaccuracy for cells close to the poles (also mention by the author) due to how trigonometric functions are used for intersection computation and the misrepresentation of the true cell shapes for everything but RLL rectangular grids (Taylor, 2024). This limits the use of the SCRIP library in my opinion to remapping between two RLL rectangular grids, which could be implemented much simpler and more accurate.
Specific comments:
Line 3: missing reference for "1987" and "1999"
Line 3-4: duplicated "has"
Line 5: missing reference for "pioneer study"
Line 18: Conservative interpolation is just one of multiple methods being used in ESM’s. However, this sentence implies that it is the most commonly used method.
Line 34-37: The CDO’s use code extracted from YAC for the first order conservative interpolation. (see Taylor, 2024)
Line 42: Give Equation number from J99 ("(10)"?).
Section "Introduction": This sections could mention other issues of SCRIP and solutions to this implemented in other software (see General remarks).
Line 73 Eq. (2): "r" is not described
Figure 2: Instead of showing the actual source code, a mathematical description of what the code does might be easier to understand.
Line 245: This may not only happen at the poles. A cell with the longitude bounds of [179;181] may be represented by [179;-179] independent of that latitude.
Line 205: The title could be more concise. In general this paragraph contains in my view too much speculations and opinions of the author. It could probably be shortened without loosing significant information.
Line 256: The implementations of trigonometric function can be more accurate for small absolute values. This may be another reason for this code in SCRIP.
Line 366, 375, 381: repeated use of "naturally...within the cell"
Section "Additional remarks": In this section the discussion of the results for ESMF and other software tested by Valcke, 2022 assumes that it is based on J99, which is not the case.
Line 380: "less impact than a change in magnitude" this is not further explained and quantified. Is this due to numerical inaccuracy of SCRIP for cells close to the poles?
Line 399: "which has the maximum deviation from the pivot longitude within a source cell" this has already been explained above
Line 403: Maybe you could explicitly mention that Figure 4(c) shows deviations from the exact fields, which are independent from the issue discussed in this paper. Additionally, you could state that in order to be able to visualize this issue, you have to compute the difference (d-b/f-b) instead of (d-a/f-a).
Line 522: duplicated "this"
Line 578,585,594: duplication of "https://doi.org"
References:
Kritsikis, E., Aechtner, M., Meurdesoif, Y., and Dubos, T.: Conservative interpolation between general spherical meshes, Geosci. Model Dev., 10, 425–431, https://doi.org/10.5194/gmd-10-425-2017, 2017.
Taylor, K. E.: Truly conserving with conservative remapping methods, Geosci. Model Dev., 17, 415–430, https://doi.org/10.5194/gmd-17-415-2024, 2024.
Valcke S, Piacentini A, Jonville G. Benchmarking Regridding Libraries Used in Earth System Modelling. Mathematical and Computational Applications. 2022; 27(2):31. https://doi.org/10.3390/mca27020031
Citation: https://doi.org/10.5194/egusphere-2024-1101-RC1
Data sets
Resources of Saito (submitted to GMD) — software and experiment data archives Fuyuki Saito https://doi.org/10.5281/zenodo.10892796
SCRIP-p (p is for pivot) F. Saito https://github.com/saitofuyuki/scrip-p
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
151 | 41 | 9 | 201 | 4 | 3 |
- HTML: 151
- PDF: 41
- XML: 9
- Total: 201
- BibTeX: 4
- EndNote: 3
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1