Fixed issue when adding multiple rivers to the same DEM

src/geofabrics/dem.py

@@ -1257,11 +1257,18 @@ def interpolate_rivers(
             region_to_rasterise.dissolve().buffer(self.catchment_geometry.resolution),
             drop=True,
         )
+        self._write_netcdf_conventions_in_place(edge_dem, self.catchment_geometry.crs)
+        edge_dem["z"] = edge_dem.z.rio.interpolate_na(method="nearest")
+        edge_dem = self._dem.rio.clip(
+            region_to_rasterise.dissolve().buffer(self.catchment_geometry.resolution),
+            drop=True,
+        )
         edge_dem = edge_dem.rio.clip(
             region_to_rasterise.dissolve().geometry,
             invert=True,
             drop=True,
         )
+
         # Define the river and mouth edge points
         grid_x, grid_y = numpy.meshgrid(edge_dem.x, edge_dem.y)
         flat_x = grid_x.flatten()

src/geofabrics/processor.py

@@ -1289,7 +1289,7 @@ def add_hydrological_features(
                     ),
                     cache_path=temp_folder,
                 )
-                temp_file = temp_folder / "dem_added_rivers.nc"
+                temp_file = temp_folder / f"dem_added_{index + 1}_rivers.nc"
                 self.logger.info(
                     f"Save temp DEM with rivers added to netCDF: {temp_file}"
                 )