Arrays passed to Fortran are no longer transposed

nens · Nov 19, 2024 · bab23cd · bab23cd
1 parent eff05d8
commit bab23cd
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/threedigrid_builder/grid/grid.py b/threedigrid_builder/grid/grid.py
@@ -992,11 +992,11 @@ def apply_cutlines(self, cutlines, dem_path):
         with Dataset(array, **dataset_kwargs) as dataset:
             gdal.RasterizeLayer(dataset, (1,), layer, options=["ATTRIBUTE=id"])
 
-        fortran_fragment_mask = np.asfortranarray(array[0].T)
+        fortran_fragment_mask = np.asfortranarray(array[0])
         assert fortran_fragment_mask.min() == no_data_value  # temp assert for mypy
 
         print(node_fragment_array)
-        fortran_node_fragment_array = np.asfortranarray(node_fragment_array.T)
+        fortran_node_fragment_array = np.asfortranarray(node_fragment_array)
         assert (
             fortran_node_fragment_array.min() == no_data_value
         )  # temp assert for mypy