bugfix in merge on mesh

Deltares · Mar 1, 2024 · cf5ae11 · cf5ae11
1 parent 6379e20
commit cf5ae11
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/hydromt_sfincs/workflows/merge.py b/hydromt_sfincs/workflows/merge.py
@@ -340,9 +340,6 @@ def merge_multi_dataarrays_on_mesh(
     )
     uda_out = xu.UgridDataArray(da0, mesh2d)
 
-    # get extent of mesh
-    bbox = uda_out.ugrid.to_crs(4326).ugrid.total_bounds
-
     # combine with next dataset
     for i in range(0, len(da_list)):
         merge_method = da_list[i].get("merge_method", "first")
@@ -355,6 +352,9 @@ def merge_multi_dataarrays_on_mesh(
         reproj_method = da_list[i].get("reproj_method", resampling_method)
         da = da_list[i].get("da")
 
+        # get extent of mesh
+        bbox = uda_out.ugrid.to_crs(da.raster.crs).ugrid.total_bounds
+
         # clip before reproject
         da = da.raster.clip_bbox(bbox, buffer=2)
         if np.any(np.array(da.shape) <= 2):