From b7ffc33ee86b1b47978885112b5416cd79811d2b Mon Sep 17 00:00:00 2001
From: Uwe Schulzweida <uwe.schulzweida@mpimet.mpg.de>
Date: Thu, 20 Mar 2025 20:29:27 +0100
Subject: [PATCH] gme_grid: check if calculation of coordinates failed
---
src/grid_gme.cc | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/grid_gme.cc b/src/grid_gme.cc
index 864dc87d6..4696e13b9 100644
--- a/src/grid_gme.cc
+++ b/src/grid_gme.cc
@@ -1198,8 +1198,9 @@ gme_grid(int withBounds, size_t gridsize, double *rlon, double *rlat, double *bl
size_t inull = 0;
for (size_t i = 0; i < gridsize; ++i)
{
- if ((xn[i * 3 + 0] == 0.0 || std::isnan(xn[i * 3 + 0])) && (xn[i * 3 + 1] == 0.0 || std::isnan(xn[i * 3 + 1]))
- && (xn[i * 3 + 2] == 0.0 || std::isnan(xn[i * 3 + 2])))
+ if ((std::fabs(xn[i * 3 + 0]) <= 0 || std::isnan(xn[i * 3 + 0]))
+ && (std::fabs(xn[i * 3 + 1]) <= 0 || std::isnan(xn[i * 3 + 1]))
+ && (std::fabs(xn[i * 3 + 2]) <= 0 || std::isnan(xn[i * 3 + 2])))
inull++;
}
if (inull > gridsize / 4)
--
GitLab