Add flushing (doorspoeling) based on 'old' LHM data

pixi.toml

@@ -66,7 +66,6 @@ quarto = "*"
 quartodoc = "*"
 rasterstats = "*"
 requests = "*"
-ribasim = "==2025.4.0"
 ruff = "*"
 seaborn = "*"
 shapely = ">=2"
@@ -82,3 +81,4 @@ bokeh_helpers = { path = "src/bokeh_helpers", editable = true }
 hydamo = { path = "src/hydamo", editable = true }
 peilbeheerst_model = { path = "src/peilbeheerst_model", editable = true }
 ribasim_nl = { path = "src/ribasim_nl", editable = true }
+ribasim = { git = "https://github.com/Deltares/Ribasim.git", branch = "main", subdirectory = "python/ribasim" }

src/peilbeheerst_model/Parametrize/AmstelGooienVecht_parametrize.py

@@ -113,6 +113,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # check basin area
 ribasim_param.validate_basin_area(ribasim_model)

src/peilbeheerst_model/Parametrize/Delfland_parametrize.py

@@ -11,6 +11,7 @@
 import peilbeheerst_model.ribasim_parametrization as ribasim_param
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -114,6 +115,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # check basin area
 ribasim_param.validate_basin_area(ribasim_model)
@@ -338,6 +340,10 @@
 )
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/HollandsNoorderkwartier_parametrize.py

@@ -11,6 +11,7 @@
 import peilbeheerst_model.ribasim_parametrization as ribasim_param
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -111,6 +112,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # check basin area
 ribasim_param.validate_basin_area(ribasim_model)
@@ -369,6 +371,10 @@
 )
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/HollandseDelta_parametrize.py

@@ -11,6 +11,7 @@
 import peilbeheerst_model.ribasim_parametrization as ribasim_param
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -112,6 +113,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # model specific tweaks
 # merge small basins into larger basins for numerical stability
@@ -738,6 +740,10 @@
     ribasim_model.basin.time.df.precipitation /= meteo_factor  # decrease meteo
     ribasim_model.basin.time.df.potential_evaporation /= meteo_factor  # decrease meteo
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 ribasim_model.use_validation = True
 ribasim_model.starttime = starttime

src/peilbeheerst_model/Parametrize/Rijnland_parametrize.py

@@ -12,6 +12,7 @@
 from peilbeheerst_model import supply
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -110,6 +111,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 inlaat_pump = []
 
@@ -302,6 +304,10 @@
 
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/Rivierenland_parametrize.py

@@ -12,6 +12,7 @@
 import peilbeheerst_model.ribasim_parametrization as ribasim_param
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -110,6 +111,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # check basin area
 ribasim_param.validate_basin_area(ribasim_model)
@@ -467,6 +469,10 @@
 )
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/Scheldestromen_parametrize.py

@@ -11,6 +11,7 @@
 import peilbeheerst_model.ribasim_parametrization as ribasim_param
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -109,6 +110,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # check basin area
 ribasim_param.validate_basin_area(ribasim_model)
@@ -325,6 +327,9 @@
 )
 ribasim_model = assign.assign_authorities()
 
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/SchielandendeKrimpenerwaard.py

@@ -12,6 +12,7 @@
 from peilbeheerst_model import supply
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -113,6 +114,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # model specific tweaks
 # change unknown streefpeilen to a default streefpeil
@@ -450,6 +452,10 @@
 )
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # set numerical settings
 # write model output
 ribasim_model.use_validation = True

src/peilbeheerst_model/Parametrize/WetterskipFryslan_parametrize.py

@@ -12,6 +12,7 @@
 from peilbeheerst_model import supply
 from peilbeheerst_model.add_storage_basins import AddStorageBasins
 from peilbeheerst_model.assign_authorities import AssignAuthorities
+from peilbeheerst_model.assign_flushing import Flushing
 from peilbeheerst_model.assign_parametrization import AssignMetaData
 from peilbeheerst_model.controle_output import Control
 from peilbeheerst_model.ribasim_feedback_processor import RibasimFeedbackProcessor
@@ -111,6 +112,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # Uitlaat toevoegen at Ter Schelling
 level_boundary_node = ribasim_model.level_boundary.add(
@@ -580,6 +582,10 @@
 )
 ribasim_model = assign.assign_authorities()
 
+# Add flushing data
+flush = Flushing(ribasim_model)
+flush.add_flushing()
+
 # TEMP CHANGES! VERY IMPORTANT TO REMOVE THIS AFTERWARDS! #@TODO ##################################
 reduce_computation_time = False
 if reduce_computation_time:

src/peilbeheerst_model/Parametrize/Zuiderzeeland_parametrize.py

@@ -111,6 +111,7 @@
 with warnings.catch_warnings():
     warnings.simplefilter(action="ignore", category=FutureWarning)
     ribasim_model = Model(filepath=ribasim_work_dir_model_toml)
+    ribasim_model.set_crs("EPSG:28992")
 
 # merge the smallest basins together
 ribasim_model.merge_basins(node_id=30, to_node_id=29)  # 4363 m2