Merge pull request #238 from EcoExtreML/SSM_v.0.3.3

Adjust to let total soil thickness and max root depth be varied and not fixed

CHANGELOG.md

@@ -16,9 +16,12 @@
 - The bug in the QVT calculations discussed in
   [#230](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/230) and fixed in
   [#234](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/234)
-- The bug in  activating the dry air calculations discussed in
+- The bug in activating the dry air calculations discussed in
   [#227](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/230) and fixed in
   [#234](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/227)
+- Defining the indices of the first four layers discussed in
+  [#237](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/237) and fixed in
+  [#238](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/238)
 
 
 <a name="1.5.0"></a>

src/+soilmoisture/calculateBoundaryConditions.m

@@ -72,7 +72,7 @@
         % Check applied infiltration doesn't exceed infiltration capacity
         topThick = 5; % first 5 cm of the soil
         satCap = SoilProperties.theta_s0 * topThick; % saturation capacity represented by saturated water content of the top 5 cm of the soil
-        actTheta = ModelSettings.DeltZ(51:54) * SoilVariables.Theta_UU(51:54, 1); % actual moisture of the top 5 cm of the soil
+        actTheta = ModelSettings.DeltZ(end - 3:end) * SoilVariables.Theta_UU(end - 4:end - 1, 1); % actual moisture of the top 5 cm of the soil
         infCap = (satCap - actTheta) / TimeProperties.DELT; % (cm/sec)
         infCap_min = min(Ks0, infCap);
 

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;
 

src/calc_sfactor.m

@@ -1,5 +1,5 @@
 function [sfactor] = calc_sfactor(Rl, Theta_s, Theta_r, Theta_LL, bbx, Ta, Theta_f)
-    SMC = Theta_LL(1:54, 1); % soil surface moisture
+    SMC = Theta_LL(1:end - 1, 1); % soil surface moisture
     wfrac = 1 ./ (1 + exp((-100 .* Theta_s') .* (SMC - (Theta_f' + Theta_r') / 2))) .* bbx; % The soil water stress factor in each layer
     RL = Rl .* bbx;
     RLfrac = RL ./ (sum(sum(RL))); % root fraction in each layer