Rename ksatprofile

to `saturated_hydraulic_conductivity_profile`. Move this key from `config.vertical` to `config.model`.
Deltares · Dec 10, 2024 · ef64378 · ef64378
1 parent 26d80ab
commit ef64378
Show file tree

Hide file tree

Showing 4 changed files with 9 additions and 9 deletions.
diff --git a/src/routing/subsurface.jl b/src/routing/subsurface.jl
@@ -57,7 +57,7 @@ function LateralSsfParameters(
     (; theta_s, theta_r, soilthickness) = soil
     soilthickness = soilthickness .* 0.001
 
-    kh_profile_type = get(config.input, "ksat_profile", "exponential")::String
+    kh_profile_type = get(config.model, "saturated_hydraulic_conductivity_profile", "exponential")::String
     dt = Second(config.timestepsecs) / basetimestep
     if kh_profile_type == "exponential"
         (; kv_0, f) = soil.kv_profile

diff --git a/src/sbm_model.jl b/src/sbm_model.jl
@@ -147,7 +147,7 @@ function initialize_sbm_model(config::Config)
             y_length,
         )
         # update variables `ssf`, `ssfmax` and `kh` (layered profile) based on ksat_profile
-        kh_profile_type = get(config.input, "ksat_profile", "exponential")::String
+        kh_profile_type = get(config.model, "saturated_hydraulic_conductivity_profile", "exponential")::String
         if kh_profile_type == "exponential" || kh_profile_type == "exponential_constant"
             initialize_lateralssf!(subsurface_flow, subsurface_flow.parameters.kh_profile)
         elseif kh_profile_type == "layered" || kh_profile_type == "layered_exponential"

diff --git a/src/soil/soil.jl b/src/soil/soil.jl
@@ -211,7 +211,7 @@ function sbm_kv_profiles(
     sumlayers,
     dt,
 )
-    kv_profile_type = get(config.input, "ksat_profile", "exponential")::String
+    kv_profile_type = get(config.model, "saturated_hydraulic_conductivity_profile", "exponential")::String
     n = length(indices)
     if kv_profile_type == "exponential"
         kv_profile = KvExponential(kv_0, f)
@@ -268,9 +268,9 @@ function sbm_kv_profiles(
             )
         end
     else
-        error("""An unknown "ksat_profile" is specified in the TOML file ($ksat_profile).
-              This should be "exponential", "exponential_constant", "layered" or
-              "layered_exponential".
+        error("""An unknown "saturated_hydraulic_conductivity_profile" is specified in the
+              TOML file ($ksat_profile). This should be "exponential",
+              "exponential_constant", "layered" or "layered_exponential".
               """)
     end
     return kv_profile

diff --git a/test/run_sbm.jl b/test/run_sbm.jl
@@ -476,7 +476,7 @@ Wflow.close_files(model; delete_output = false)
     end
 
     @testset "exponential constant profile" begin
-        config.input.ksat_profile = "exponential_constant"
+        config.model.saturated_hydraulic_conductivity_profile = "exponential_constant"
         model = Wflow.initialize_sbm_model(config)
         (; soil) = model.vertical
         (; kv_profile) = soil.parameters
@@ -497,7 +497,7 @@ Wflow.close_files(model; delete_output = false)
     end
 
     @testset "layered profile" begin
-        config.input.ksat_profile = "layered"
+        config.model.saturated_hydraulic_conductivity_profile = "layered"
         model = Wflow.initialize_sbm_model(config)
         (; soil) = model.vertical
         (; kv_profile) = soil.parameters
@@ -513,7 +513,7 @@ Wflow.close_files(model; delete_output = false)
     end
 
     @testset "layered exponential profile" begin
-        config.input.ksat_profile = "layered_exponential"
+        config.model.saturated_hydraulic_conductivity_profile = "layered_exponential"
         model = Wflow.initialize_sbm_model(config)
         (; soil) = model.vertical
         (; kv_profile) = soil.parameters