adjust to let total soil thickness and max root depth not fixed in th…

…e code
EcoExtreML · Jul 3, 2024 · f96ca1c · f96ca1c
1 parent 3fa65aa
commit f96ca1c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/calc_rsoil.m b/src/calc_rsoil.m
@@ -2,7 +2,7 @@
     % load Constants
     Constants = io.define_constants();
 
-    SMC = SoilVariables.Theta_LL(1:54, 2);
+    SMC = SoilVariables.Theta_LL(1:end - 1, 2);
     Se  = (SMC - VanGenuchten.Theta_r') ./ (VanGenuchten.Theta_s' - VanGenuchten.Theta_r');
     Ksoil = SoilVariables.Ks' .* Se.^Constants.l .* (1 - (1 - Se.^(1 ./ VanGenuchten.m')).^(VanGenuchten.m')).^2;