From 2c1ecdbdb7441ebc56c66341fcba2c035663fbca Mon Sep 17 00:00:00 2001 From: Martijn Visser Date: Tue, 6 Feb 2024 22:28:55 +0100 Subject: [PATCH 1/2] Correct precipitation equation --- docs/core/equations.qmd | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/docs/core/equations.qmd b/docs/core/equations.qmd index 92d5d258f..163717eec 100644 --- a/docs/core/equations.qmd +++ b/docs/core/equations.qmd @@ -149,19 +149,21 @@ plt.show() The precipitation term is given by $$ - Q_P = P \cdot A(u). + Q_P = P \cdot A $$ {#eq-precip} -Here $P = P(t)$ is the precipitation rate and $A$ is the wetted area. $A$ is a -function of the storage $u = u(t)$: as the volume of water changes, the area of the free water -surface may change as well, depending on the slopes of the surface waters. +Here $P = P(t)$ is the precipitation rate and $A$ is the maximum area given in the `Basin / profile` table. +Precipitation in the Basin area is assumed to be directly added to the Basin storage. +The modeler needs to ensure all precipitation enters the model, and there is no overlap in the maximum profile areas, else extra water is created. +If a part of the catchment is not in any Basin profile, the modeler has to verify that water source is not forgotten. +It can for instance be converted to a flow rate and added to a Basin as a FlowBoundary. ## Evaporation The evaporation term is given by $$ - Q_E = E_\text{pot} \cdot A(u) \cdot \phi(d;0.1). + Q_E = E_\text{pot} \cdot A(u) \cdot \phi(d;0.1) $$ {#eq-evap} Here $E_\text{pot} = E_\text{pot}(t)$ is the potential evaporation rate and $A$ is the wetted area. $\phi$ is the [reduction factor](equations.qmd#sec-reduction_factor) which depends on the depth $d$. It provides a smooth gradient as $u \rightarrow 0$. From c6c002d46265bd69fd74ed39fc7ee64c6f893aef Mon Sep 17 00:00:00 2001 From: Martijn Visser Date: Fri, 9 Feb 2024 15:13:52 +0100 Subject: [PATCH 2/2] add back dots Except a misplaced one that breaks the quarto formatting. --- docs/core/equations.qmd | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/docs/core/equations.qmd b/docs/core/equations.qmd index 163717eec..2546c454e 100644 --- a/docs/core/equations.qmd +++ b/docs/core/equations.qmd @@ -149,7 +149,7 @@ plt.show() The precipitation term is given by $$ - Q_P = P \cdot A + Q_P = P \cdot A. $$ {#eq-precip} Here $P = P(t)$ is the precipitation rate and $A$ is the maximum area given in the `Basin / profile` table. @@ -163,7 +163,7 @@ It can for instance be converted to a flow rate and added to a Basin as a FlowBo The evaporation term is given by $$ - Q_E = E_\text{pot} \cdot A(u) \cdot \phi(d;0.1) + Q_E = E_\text{pot} \cdot A(u) \cdot \phi(d;0.1). $$ {#eq-evap} Here $E_\text{pot} = E_\text{pot}(t)$ is the potential evaporation rate and $A$ is the wetted area. $\phi$ is the [reduction factor](equations.qmd#sec-reduction_factor) which depends on the depth $d$. It provides a smooth gradient as $u \rightarrow 0$. @@ -194,7 +194,7 @@ MODFLOW 6 boundary conditions in the basin: $$ Q_\text{inf} = \sum_{i=1}^{n} \sum_{j=1}^{m} \max(Q_{\mathrm{mf6}_{i,j}}, 0.0) -$$ {#eq-inf}. +$$ {#eq-inf} Where $i$ is the index of the boundary condition, $j$ the MODFLOW 6 cell index, $n$ the number of boundary conditions, and $m$ the number of MODFLOW 6 cells in