Correct precipitation equation #1088
Conversation
docs/core/equations.qmd
Outdated
|
|
||
| ## 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) |
There was a problem hiding this comment.
In a scientific writing course I was thought that equations in this context must be treated as part of the sentence, and so there should be a period at the end of the equation.
There was a problem hiding this comment.
Interesting, I didn't know. I added them back. There do seems to be quite some equations without dots though.
| 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. |
There was a problem hiding this comment.
Is it really the maximum area that is used? Considering
fixed_area = basin.area[i][end]it depends on the ordering of these vectors.
There was a problem hiding this comment.
Indeed it depends on the ordering, but that is enforced to be increasing. We rely on these orderings in many places in the code, so I don't think an explicit max is needed.
https://deltares.github.io/Ribasim/core/usage.html#basin-profile
Except a misplaced one that breaks the quarto formatting.
Before we used the current wetted area to convert the precipitation$m/s$ input to a flow rate $m^3/s$ . It was deemed simpler and more stable to use the maximum profile area instead, otherwise empty Basins with 0 wetted area could never fill back up due to rain. This updates the equations to reflect that.