Added riverbed roughness parameter for MCT

ec-jrc · Sep 13, 2024 · c2e8e15 · c2e8e15
1 parent 1a46d6c
commit c2e8e15
Show file tree

Hide file tree

Showing 3 changed files with 16 additions and 3 deletions.
diff --git a/src/lisflood/hydrological_modules/routing.py b/src/lisflood/hydrological_modules/routing.py
@@ -49,7 +49,7 @@ class routing(HydroModule):
                                 'ChanSdXdY', 'TotalCrossSectionAreaInitValue', 'PrevDischarge'],
                         'SplitRouting': ['CrossSection2AreaInitValue', 'PrevSideflowInitValue', 'CalChanMan2'],
                         'dynamicWave': ['ChannelsDynamic'],
-                        'MCTRouting': ['ChannelsMCT', 'ChanGradMaxMCT', 'PrevCmMCTInitValue', 'PrevDmMCTInitValue']}
+                        'MCTRouting': ['ChannelsMCT', 'ChanGradMaxMCT', 'PrevCmMCTInitValue', 'PrevDmMCTInitValue','CalChanMan3']}
     module_name = 'Routing'
 
     def __init__(self, routing_variable):
@@ -568,7 +568,7 @@ def initialMCT(self):
 
             # cmcheck
             # This could become a calibration parameter if we want to use MCT+SplitRouting
-            self.var.ChanManMCT = (self.var.ChanMan / self.var.CalChanMan) * loadmap('CalChanMan2')
+            self.var.ChanManMCT = (self.var.ChanMan / self.var.CalChanMan) * loadmap('CalChanMan3')
             # Mannings coefficient for MCT pixels (same as second line of split routing)
 
             PrevCmMCT = loadmap('PrevCmMCTInitValue')

diff --git a/src/lisfloodSettings_reference.xml b/src/lisfloodSettings_reference.xml
@@ -439,6 +439,13 @@ Multiplier [-] applied to Channel Manning's n for second line of routing
 </comment>
 </textvar>
 
+<textvar name="CalChanMan3" value="$(PathParams)/params_CalChanMan3">
+<comment>
+default: 3.0 [-]
+Multiplier [-] applied to Channel Manning's n for MCT routing
+</comment>
+</textvar>
+
 <textvar name="LakeMultiplier" value="$(PathParams)/params_LakeMultiplier">
 <comment>
 default: 1.0 [-]
@@ -1446,7 +1453,7 @@ Minimum slope gradient [-] (for kin. wave: slope cannot be 0)
 </comment>
 </textvar>
 
-<textvar name="ChanGradMin" value="0.0001">
+<textvar name="ChanGradMin" value="0.00001">
 <comment>
 Minimum channel gradient [-] (for kin. wave: slope cannot be 0)
 </comment>
@@ -2130,6 +2137,12 @@ Multiplier applied to Channel Manning's n for second routing line
 </comment>
 </textvar>
 
+<textvar name="CalChanMan3" value="$(CalChanMan3)">
+<comment>
+Multiplier applied to Channel Manning's n for MCT routing
+</comment>
+</textvar>
+
 <textvar name="QSplitMult" value="$(QSplitMult)">
 <comment>
 Multiplier applied to average Q to split into a second line of routing

diff --git a/tests/data/MCT_UseCase/maps/parameters/params_CalChanMan3.nc b/tests/data/MCT_UseCase/maps/parameters/params_CalChanMan3.nc