ECC-1708: grib_get_data returns offset longitudes in some cases

ecmwf · Nov 23, 2023 · 1229481 · 1229481
1 parent a3ec82a
commit 1229481
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/src/grib_iterator_class_regular.cc b/src/grib_iterator_class_regular.cc
@@ -233,10 +233,14 @@ static int init(grib_iterator* iter, grib_handle* h, grib_arguments* args)
         self->los[loi] = lon1;
         lon1 += idir;
     }
-    /* ECC-1406: Due to rounding, errors can accumulate.
-     * So we ensure the last longitude is longitudeOfLastGridPointInDegrees
-    */
-    self->los[Ni-1] = normalise_longitude_in_degrees(lon2); // Also see ECC-1671
+
+    // ECC-1406: Due to rounding, errors can accumulate.
+    // So we ensure the last longitude is longitudeOfLastGridPointInDegrees
+    // Also see ECC-1671, ECC-1708
+    if (lon2 > 0) {
+        lon2 = normalise_longitude_in_degrees(lon2);
+    }
+    self->los[Ni-1] = lon2;
 
     return ret;
 }