maximum temperature change

emf-creaf · May 31, 2024 · b51ce9b · b51ce9b
1 parent 6af83a5
commit b51ce9b
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/transpiration_advanced.cpp b/src/transpiration_advanced.cpp
@@ -1040,7 +1040,7 @@ List transpirationAdvanced(List x, NumericVector meteovec,
 
  //Soil temperature changes
  NumericVector soilTchange = temperatureChange(widths, Tsoil, sand, clay, Ws, Theta_SAT, Theta_FC, Ebalsoil[n], tstep);
- for(int l=0;l<nlayers;l++) Tsoil[l] = Tsoil[l] + soilTchange[l];
+ for(int l=0;l<nlayers;l++) Tsoil[l] = Tsoil[l] + std::max(-3.0, std::min(3.0, soilTchange[l]));
  if(n<(ntimesteps-1)) Tsoil_mat(n+1,_)= Tsoil;
 
  } else { //Multilayer canopy balance