Merge pull request #260 from EcoExtreML/revise-IS

revise IS, relates to (#252)
EcoExtreML · Sep 25, 2024 · 56f2c81 · 56f2c81
2 parents 1c15e9a + 9db52b5
commit 56f2c81
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 2 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -4,6 +4,8 @@
 
 
 **Fixed:**
+- The bug in calculation of IS discussed in
+  [#252](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/252)
 
 
 <a name="1.6.0"></a>

diff --git a/run_model_on_snellius/exe/STEMMUS_SCOPE b/run_model_on_snellius/exe/STEMMUS_SCOPE
diff --git a/src/+init/soilHeteroSubroutine.m b/src/+init/soilHeteroSubroutine.m
@@ -47,12 +47,12 @@
             j = i - 1;
         end
 
-        SoilVariables.IS(i) = indexOfSoilType;
+        SoilVariables.IS(j) = indexOfSoilType;
         if subRoutine == 4
             SoilVariables.IS(5:8) = 5;
         end
 
-        J = SoilVariables.IS(i);
+        J = SoilVariables.IS(j);
         [SoilVariables, VanGenuchten] = init.updateSoilVariables(SoilVariables, VanGenuchten, SoilProperties, ModelSettings, j, J);
         SoilVariables.Imped(i) = ImpedF(J);