Added lake ids in lakes warning msg

Author: doc78

src/lisflood/hydrological_modules/lakes.py

@@ -242,7 +242,9 @@ def dynamic_inloop(self, NoRoutingExecuted):
             # Check LakeOutflowCC for negative values
             if any(LakeStorageIndicator < 0):
                 # this can happen with large oscillations in the lake inflow
-                msg = f"Negative outflow from lakes. (Total negative values: {sum(LakeStorageIndicator < 0)}) \n" \
+                msg = f"Negative outflow from lakes. \n" \
+                      f"(Total negative values: {sum(LakeStorageIndicator < 0)}, " \
+                      f"Lakes IDs: {self.var.LakeSitesCC[LakeStorageIndicator < 0]}) \n" \
                       "Consider increasing computation time step for routing (DtSecChannel) \n"
                 warnings.warn(LisfloodWarning(msg))
                 LakeStorageIndicator[LakeStorageIndicator < 0] = 0
@@ -274,8 +276,8 @@ def dynamic_inloop(self, NoRoutingExecuted):
             # Check LakeStorageM3CC for negative values and set them to zero
             if any(np.isnan(self.var.LakeStorageM3CC)) or any(self.var.LakeStorageM3CC < 0):
                 msg = "Negative or NaN volume for lake storage set to 0. \n" \
-                      f"(Total negative values: {sum(self.var.LakeStorageM3CC < 0)}, " \
-                      f"Total NaN values: {sum(np.isnan(self.var.LakeStorageM3CC))})\n" \
+                      f"Total negative values: {sum(self.var.LakeStorageM3CC < 0)}, Lakes IDs: {self.var.LakeSitesCC[self.var.LakeStorageM3CC < 0]} \n" \
+                      f"Total NaN values: {sum(np.isnan(self.var.LakeStorageM3CC))}, Lakes IDs: {self.var.LakeSitesCC[np.isnan(self.var.LakeStorageM3CC)]} \n" \
                       "Increase computation time step for routing (DtSecChannel) \n"
                 warnings.warn(LisfloodWarning(msg))
                 self.var.LakeStorageM3CC[self.var.LakeStorageM3CC < 0] = 0