From f96ca1ca8176cab732b5ba440ec206a988aeca31 Mon Sep 17 00:00:00 2001
From: Mostafa Gomaa Daoud <54531356+MostafaGomaa93@users.noreply.github.com>
Date: Wed, 3 Jul 2024 16:48:54 +0200
Subject: [PATCH] adjust to let total soil thickness and max root depth not
 fixed in the code

---
 src/calc_rsoil.m | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/calc_rsoil.m b/src/calc_rsoil.m
index 446c0a65..a37cc034 100644
--- 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;