try to fix reduced_ll grid for res 3601

ecmwf · Jan 28, 2025 · e4075b1 · e4075b1
1 parent 0860f11
commit e4075b1
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/polytope_feature/datacube/transformations/datacube_mappers/mapper_types/reduced_ll.py b/polytope_feature/datacube/transformations/datacube_mappers/mapper_types/reduced_ll.py
@@ -26,9 +26,11 @@ def first_axis_vals(self):
         end_lat = -90
         if self._resolution == 3601:
             start_lat = 89.9731
+            end_lat = -89.9731
         # resolution = 180 / (self._resolution - 1)
         # vals = [round(90 - i * resolution, 12) for i in range(self._resolution)]
-        resolution = 180 / (self._resolution - 1)
+        # resolution = 180 / (self._resolution - 1)
+        resolution = (start_lat - end_lat) / (self._resolution - 1)
         vals = [round(start_lat - i * resolution, 12) for i in range(self._resolution)]
         return vals