correction fSL

emf-creaf · Jun 12, 2024 · d53cbbc · d53cbbc
1 parent 3453a72
commit d53cbbc
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 0 deletions.
diff --git a/src/lightextinction_advanced.cpp b/src/lightextinction_advanced.cpp
@@ -363,6 +363,7 @@ NumericVector layerSunlitFraction(NumericMatrix LAIme, NumericMatrix LAImd,
       s2 += kb[j]*ClumpingIndex[j]*0.5*(LAIme(i,j)+LAImd(i,j));
     }
     fSL[i] = exp(-1.0*s1)*exp(-1.0*s2);
+    if(fSL[i]<0.00001) fSL[i] = 0.0; //Avoids overestimation of absorbed radiation per leaf area when sunlit fraction is close to zero
   }
   return(fSL);
 }

diff --git a/src/spwb.cpp b/src/spwb.cpp
@@ -706,6 +706,7 @@ List spwbDay_advanced(List x, NumericVector meteovec,
                         _["ShadeLeaves"] = transp["ShadeLeaves"],
                         _["ExtractionInst"] = soilLayerExtractInst,
                         _["PlantsInst"] = PlantsInst,
+                        _["RadiationInputInst"] = transp["RadiationInputInst"],
                         _["SunlitLeavesInst"] = transp["SunlitLeavesInst"],
                         _["ShadeLeavesInst"] = transp["ShadeLeavesInst"],
                         _["LightExtinction"] = transp["LightExtinction"],

diff --git a/src/transpiration_advanced.cpp b/src/transpiration_advanced.cpp
@@ -1256,6 +1256,7 @@ List transpirationAdvanced(List x, NumericVector meteovec,
                         _["SunlitLeaves"] = Sunlit,
                         _["ShadeLeaves"] = Shade,
                         _["ExtractionInst"] = soilLayerExtractInst,
+                        _["RadiationInputInst"] = ddd,
                         _["PlantsInst"] = PlantsInst,
                         _["SunlitLeavesInst"] = SunlitInst,
                         _["ShadeLeavesInst"] = ShadeInst,