Merge pull request #200 from EcoExtreML/fix_issue_101

remove ObservationPoints script
EcoExtreML · Oct 17, 2023 · 8e2642f · 8e2642f
2 parents 9dc322f + 02556cf
commit 8e2642f
Show file tree

Hide file tree

Showing 4 changed files with 12 additions and 23 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -10,13 +10,16 @@
   [194](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/194)
 - All `Engry_*` functions are refcatored and moved to `+energy` folder in
 [197](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/197)
+- Remove `ObservationPoints` script
+[200](https://github.com/EcoExtreML/STEMMUS_SCOPE/pull/200)
 
 **Fixed:**
 
 - issue [#181](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/181)
 - issue [#98](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/98)
 - issue [#99](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/99)
 - issue [#100](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/100)
+- issue [#101](https://github.com/EcoExtreML/STEMMUS_SCOPE/issues/101)
 
 <a name="1.3.0"></a>
 # [1.3.0](https://github.com/EcoExtreML/STEMMUS_SCOPE/releases/tag/1.3.0) - 22 Jun 2023

diff --git a/run_model_on_snellius/exe/STEMMUS_SCOPE b/run_model_on_snellius/exe/STEMMUS_SCOPE
diff --git a/src/ObservationPoints.m b/src/ObservationPoints.m
diff --git a/src/STEMMUS_SCOPE.m b/src/STEMMUS_SCOPE.m
@@ -352,6 +352,12 @@
 Delt_t0 = Delt_t; % Duration of last time step
 TOLD_CRIT = [];
 
+% for soil moisture and temperature outputs
+monitorDepthTemperature = ModelSettings.NL:-1:1;
+monitorDepthSoilMoisture = ModelSettings.NL:-1:1;
+Sim_Theta_U = [];
+Sim_Temp = [];
+
 % Srt, root water uptake;
 Srt = InitialValues.Srt;  % will be updated!
 P_gg = InitialValues.P_gg;  % will be updated!
@@ -736,7 +742,9 @@
             SoilVariables.Theta_U = Theta_U;
             SoilVariables.Theta_I = Theta_I;
 
-            run ObservationPoints;
+            % replace run ObservationPoints, see issue 101
+            Sim_Theta_U(KT, 1:length(monitorDepthSoilMoisture)) = Theta_UUU(monitorDepthSoilMoisture, 1, KT);
+            Sim_Temp(KT, 1:length(monitorDepthTemperature)) = TTT(monitorDepthTemperature, KT);
         end
         if (TEND - TIME) < 1E-3
             for MN = 1:NN