Fix grc structure (model for reunion). remove leaking from slow trans…

…fert reservoir cl, because the model calibration was done with this structure
DassHydro · Nov 13, 2024 · b0fff3f · b0fff3f
1 parent 9a5587a
commit b0fff3f
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/doc/source/math_num_documentation/forward_structure.rst b/doc/source/math_num_documentation/forward_structure.rst
@@ -666,7 +666,7 @@ Hydrological processes can be described at pixel scale in `smash` with one of th
 
         \begin{eqnarray}
             &p_{rr}(x, t)& &=& &0.6 \times 0.9(p_r(x, t) + p_{erc}(x, t)) + l_{exc}(x, t)\\
-            &p_{rl}(x, t)& &=& &0.4 \times 0.9(p_r(x, t) + p_{erc}(x, t)) + l_{exc}(x, t)\\
+            &p_{rl}(x, t)& &=& &0.4 \times 0.9(p_r(x, t) + p_{erc}(x, t)) \\
             &p_{rd}(x, t)& &=& &0.1(p_r(x, t) + p_{erc}(x, t))
         \end{eqnarray}
 

diff --git a/smash/fcore/operator/md_gr_operator.f90 b/smash/fcore/operator/md_gr_operator.f90
@@ -1577,7 +1577,7 @@ subroutine grc_time_step(setup, mesh, input_data, options, returns, time_step, a
                 end if
 
                 prr = 0.6_sp*0.9_sp*(pr + perc) + l
-                prl = 0.4_sp*0.9_sp*(pr + perc) + l
+                prl = 0.4_sp*0.9_sp*(pr + perc)
                 prd = 0.1_sp*(pr + perc)
 
                 call gr_transfer(5._sp, ac_prcp(k), prr, ac_ct(k), ac_ht(k), qr)