From 1e44e7eaa3822859fcad745f4ca18162e3e0c50a Mon Sep 17 00:00:00 2001 From: asjeb <151923475+asjeb@users.noreply.github.com> Date: Fri, 24 Jan 2025 17:47:47 +0100 Subject: [PATCH] DOC: imperviousness (#368) * math-num imperviousness section * remark of reviewers --------- Co-authored-by: Apolline Elbaz --- .../_static/gr4_structure_imperviousness.svg | 2154 +++++++++++++++++ .../forward_structure.rst | 47 +- smash/fcore/derived_type/mwd_physio_data.f90 | 1 + 3 files changed, 2201 insertions(+), 1 deletion(-) create mode 100644 doc/source/_static/gr4_structure_imperviousness.svg diff --git a/doc/source/_static/gr4_structure_imperviousness.svg b/doc/source/_static/gr4_structure_imperviousness.svg new file mode 100644 index 00000000..4c989223 --- /dev/null +++ b/doc/source/_static/gr4_structure_imperviousness.svg @@ -0,0 +1,2154 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/doc/source/math_num_documentation/forward_structure.rst b/doc/source/math_num_documentation/forward_structure.rst index 12015ed0..9df9dc4e 100644 --- a/doc/source/math_num_documentation/forward_structure.rst +++ b/doc/source/math_num_documentation/forward_structure.rst @@ -689,6 +689,8 @@ Hydrological processes can be described at pixel scale in `smash` with one of th q_t(x, t) = q_r(x, t) + + .. _math_num_documentation.forward_structure.hydrological_module.loieau: .. dropdown:: Génie Rural LoiEau (loieau) @@ -791,7 +793,6 @@ Hydrological processes can be described at pixel scale in `smash` with one of th q_t(x, t) = k_b(x)\left(q_r(x, t) + q_d(x, t)\right) - .. _math_num_documentation.forward_structure.hydrological_module.gr_rainfall_intensity: .. dropdown:: Génie Rural with rainfall intensity terms (gr4_ri, gr5_ri) @@ -1020,6 +1021,50 @@ Hydrological processes can be described at pixel scale in `smash` with one of th Same as ``gr4_ri`` transfer, see :ref:`GR4 Transfer ` +.. _math_num_documentation.forward_structure.hydrological_module.imperviousness: + +.. dropdown:: Génie Rural with imperviousness + :animate: fade-in-slide-down + + This imperviousness feature allows for the calculation of the impervious proportion of a pixel's surface and takes this into account when computing infiltration and evaporation fluxes applied to the GR type production reservoir. + The imperviousness coefficients :math:`imperv(x)` influence the fluxes of the production reservoir of each cell by being applied to the neutralized rainfall :math:`p_n(x,t)` and the evaporation :math:`e_s(x,t)`. + The imperviousness coefficients must range between 0 and 1 and be specified through an input map that is consistent with the model grid. This map can be obtained, for example, from soil occupation processing. + For instance, if the imperviousness coefficient is close to 1, the production part receives less neutralized rainfall :math:`p_n` and there is less evaporation :math:`e_s` from the impermeable soil. + This imperviousness accounting for the GR reservoir is applicable to GR model structures in SMASH. This is illustrated here on the GR4 structure. + + .. figure:: ../_static/gr4_structure_imperviousness.svg + :align: center + :width: 300 + + Diagram of the ``gr4`` hydrological operator with imperviousness, a simplified ``GR`` like model for spatialized modeling. + + + **Production** + + - Compute the neutralized precipitation :math:`p_n` on impermeable soil + + .. math:: + :nowrap: + + \begin{eqnarray} + + &p_n(x, t)& &=& & \left(1 - imperv(x)\right)\ p_n(x, t) + + \end{eqnarray} + + - Compute the production infiltrating precipitation :math:`p_s` and evaporation :math:`e_s` + + .. math:: + :nowrap: + + \begin{eqnarray} + + &p_s(x, t)& &=& &c_p(x) (1 - h_p(x, t - 1)^2) \frac{\tanh\left(\frac{p_n(x, t)}{c_p(x)}\right)}{1 + h_p(x, t - 1) \tanh\left(\frac{p_n(x, t)}{c_p(x)}\right)}\\ + + &e_s(x, t)& &=& &(1 - imperv(x)) \left(h_p(x, t - 1) c_p(x) (2 - h_p(x, t - 1)) \frac{\tanh\left(\frac{e_n(x, t)}{c_p(x)}\right)}{1 + (1 - h_p(x, t - 1)) \tanh\left(\frac{e_n(x, t)}{c_p(x)}\right)} \right) + \end{eqnarray} + + .. _math_num_documentation.forward_structure.hydrological_module.vic3l: diff --git a/smash/fcore/derived_type/mwd_physio_data.f90 b/smash/fcore/derived_type/mwd_physio_data.f90 index 8f24613e..f92f0618 100644 --- a/smash/fcore/derived_type/mwd_physio_data.f90 +++ b/smash/fcore/derived_type/mwd_physio_data.f90 @@ -10,6 +10,7 @@ !% `Variables` Description !% ======================== ======================================= !% ``descriptor`` Descriptor maps field [(descriptor dependent)] +!% ``imperviousness`` Imperviousness map !% ``l_descriptor`` Descriptor maps field min value [(descriptor dependent)] !% ``u_descriptor`` Descriptor maps field max value [(descriptor dependent)] !% ======================== =======================================