diff --git a/src/lightextinction_advanced.cpp b/src/lightextinction_advanced.cpp
index 8a7123e50..c3ca59e4b 100644
--- a/src/lightextinction_advanced.cpp
+++ b/src/lightextinction_advanced.cpp
@@ -353,7 +353,7 @@ NumericVector layerSunlitFraction(NumericMatrix LAIme, NumericMatrix LAImd,
     double s2=0.0;
     for(int j = 0;j<ncoh;j++) {
       s1 += kb[j]*ClumpingIndex[j]*(LAIme(i,j)+LAImd(i,j));
-      s2 += kb[j]*ClumpingIndex[j]*(LAIme(i,j)+LAImd(i,j));
+      s2 += kb[j]*ClumpingIndex[j]*0.5*(LAIme(i,j)+LAImd(i,j));
     }
     fSL[i] = exp(-1.0*s1)*exp(-1.0*s2);
   }
diff --git a/vignettes/evaluation/StandLevelEvaluation.Rmd b/vignettes/evaluation/StandLevelEvaluation.Rmd
index f0ff1f6e7..fe3907d0e 100644
--- a/vignettes/evaluation/StandLevelEvaluation.Rmd
+++ b/vignettes/evaluation/StandLevelEvaluation.Rmd
@@ -1,6 +1,6 @@
 ---
 title: "Model evaluation in experimental plots"
-author: "Miquel De Cáceres (CREAF), Víctor Granda (CREAF),  Rafael Poyatos (CREAF), Teresa Gimeno (CREAF), Nicolas Martin-StPaul (INRAE), Antoine Cabon (WSL)"
+author: "Miquel De Cáceres (CREAF), Víctor Granda (CREAF),  Rafael Poyatos (CREAF), Teresa Gimeno (CREAF), Nicolas Martin-StPaul (INRAE), Antoine Cabon (WSL), Vincenzo Saponaro (U. Tuscia)"
 date: "`r Sys.Date()`"
 output: 
   rmarkdown::html_document:
@@ -13,10 +13,11 @@ output:
 ```{r setup, include=FALSE}
 params = list(wd='/home/miquel/OneDrive/mcaceres_work/model_development/medfate_development/StandLevelEvaluation/',
               model='spwb',
-              # names = c("Font-Blanche"),
-              # codes = c("FONBLA"),
+              # names = c("Collelongo"),
+              # codes = c("COLLELONGO"),
               names=c("Wombat","Euc-FACE",
                       "Puéchabon", "Hesse","Fontainebleau-Barbeau",  "Font-Blanche",
+                      "Collelongo",
                       "Yatir",
                       "Mitra",
                       "Rinconada", "Vallcebre (Barrol)", "Vallcebre (Sort)", "Prades", "Can Balasc", "Alto-Tajo Armallones", "Ronda",
@@ -24,6 +25,7 @@ params = list(wd='/home/miquel/OneDrive/mcaceres_work/model_development/medfate_
                       "Morgan-Monroe", "Sevilleta"),
               codes=c("WOMBAT", "EUCFACE",
                       "FRAPUE", "FRAHES", "FONTAINEBLEAU","FONBLA",
+                      "COLLELONGO",
                       "ISRYAT",
                       "MITRA",
                       "ESPRIN", "QVALLCEBRE", "PVALLCEBRE", "PRADES", "CANBALASC", "ESPALTARM", "RONDA",
diff --git a/vignettes/evaluation/_Stand_Evaluation_child.Rmd b/vignettes/evaluation/_Stand_Evaluation_child.Rmd
index 09f192d03..f41891d92 100644
--- a/vignettes/evaluation/_Stand_Evaluation_child.Rmd
+++ b/vignettes/evaluation/_Stand_Evaluation_child.Rmd
@@ -117,7 +117,7 @@ eval_SWC <- "SWC" %in% names(df_site$SiteData[[1]]$measuredData)
 
 <!-- Plant transpiration -->
 ```{r, echo=FALSE, results='asis', warning = FALSE, message=FALSE}
-eval_Eplant <- TRUE
+eval_Eplant <- any(startsWith(names(df_site$SiteData[[1]]$measuredData),"E_"))
  if(eval_Eplant) {
    res <- knitr::knit_child(file.path(this_path,'_Eplant_Evaluation_child.Rmd'), quiet = TRUE)
    cat(res, sep = '\n')
diff --git a/vignettes/runmodels/AdvancedWaterEnergyBalance.Rmd b/vignettes/runmodels/AdvancedWaterEnergyBalance.Rmd
index 76c0a4252..d077c56c8 100644
--- a/vignettes/runmodels/AdvancedWaterEnergyBalance.Rmd
+++ b/vignettes/runmodels/AdvancedWaterEnergyBalance.Rmd
@@ -61,7 +61,7 @@ Apart from data inputs, the behaviour of simulation models is controlled using a
 control <- defaultControl("Sperry")
 ```
 
-To use the advanced water balance model we must change the values of `transpirationMode` to switch from `"Granier"` to either `"Sperry"` or `"Cochard"`. 
+To use the advanced water balance model we must change the values of `transpirationMode` to switch from `"Granier"` to either `"Sperry"` or `"Sureau"`. 
 
 ### Water balance input object
 A last object is needed before calling simulation functions, called `spwbInput`. It consists in the compilation of aboveground, belowground parameters and the specification of additional parameter values for each plant cohort. This is done by calling function `forest2spwbInput()`: