changing roots tensile strength with additional cohesion

ARPA-SIMC · Oct 26, 2023 · 97a54e0 · 97a54e0
1 parent 1352c3f
commit 97a54e0
Show file tree

Hide file tree

Showing 5 changed files with 4 additions and 4 deletions.
diff --git a/DATA/PROJECT/montue/data/crop.db b/DATA/PROJECT/montue/data/crop.db
diff --git a/agrolib/criteriaModel/criteria1DCase.cpp b/agrolib/criteriaModel/criteria1DCase.cpp
@@ -566,7 +566,7 @@ bool Crit1DCase::computeDailyModel(Crit3DDate &myDate, std::string &error)
             else
                 rootDensityNorm = crop.roots.rootDensity[l] / rootDensityMax;
 
-            double rootCohesion = crop.roots.rootTensileStrength * rootDensityNorm;
+            double rootCohesion = crop.roots.rootsAdditionalCohesion * rootDensityNorm;
 
 
             soilLayers[l].factorOfSafety = soilLayers[l].computeSlopeStability(unit.slope, rootCohesion);

diff --git a/agrolib/crop/cropDbTools.cpp b/agrolib/crop/cropDbTools.cpp
@@ -66,7 +66,7 @@ bool loadCropParameters(const QSqlDatabase &dbCrop, QString idCrop, Crit3DCrop &
     myCrop.roots.rootDepthMin = query.value("root_depth_zero").toDouble();
     myCrop.roots.rootDepthMax = query.value("root_depth_max").toDouble();
     myCrop.roots.actualRootDepthMax = myCrop.roots.rootDepthMax;
-    myCrop.roots.rootTensileStrength = query.value("root_tensile_strength").toDouble();
+    myCrop.roots.rootsAdditionalCohesion = query.value("roots_additional_cohesion").toDouble();
 
     getValue(query.value("degree_days_root_increase"), &(myCrop.roots.degreeDaysRootGrowth));
     if (myCrop.roots.degreeDaysRootGrowth == NODATA)

diff --git a/agrolib/crop/root.cpp b/agrolib/crop/root.cpp
@@ -60,7 +60,7 @@ void Crit3DRoot::clear()
     rootLength = NODATA;
     rootDepth = NODATA;
     rootDensity.clear();
-    rootTensileStrength = NODATA;
+    rootsAdditionalCohesion = NODATA;
 }
 
 

diff --git a/agrolib/crop/root.h b/agrolib/crop/root.h
@@ -33,7 +33,7 @@
         int firstRootLayer;                 /*!< [-]  */
         int lastRootLayer;                  /*!< [-]  */
         std::vector<double> rootDensity;    /*!< [-]  */
-        double rootTensileStrength;         /*!< [kPa] tensile strength of a single root */
+        double rootsAdditionalCohesion;      /*!< [kPa] Cr = roots reinforcement (RR) derived from a model */
 
         /*! state variables */
         double rootDepth;                   /*!<  [m]  current root depth */