MIR-650 Exception terminates mir-get-data: HEALPix nested ordering wi…

…th N=1
ecmwf · Mar 11, 2024 · 04178bd · 04178bd
1 parent c7f5c24
commit 04178bd
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/mir/reorder/HEALPix.cc b/src/mir/reorder/HEALPix.cc
@@ -185,8 +185,8 @@ int HEALPixReorder::ring_to_nest(int r) const {
             p -= 8 * Nside_;
         }
 
-        int i = (r + p) >> 1;
-        int j = (r - p) >> 1;
+        int i = std::max(0, (r + p)) >> 1;
+        int j = std::max(0, (r - p)) >> 1;
 
         ASSERT(f < 12 && i < Nside_ && j < Nside_);
         return CodecFijNest::fij_to_nest(f, i, j, k_);