diff --git a/docs/couple/modflow-demo.qmd b/docs/couple/modflow-demo.qmd
index 628532baf..b407a0450 100644
--- a/docs/couple/modflow-demo.qmd
+++ b/docs/couple/modflow-demo.qmd
@@ -60,11 +60,11 @@ and volume-discharge relationships.
:::
A visual representation of this simplified conceptual schematization is given
-in @fig-volume-depth and @fig-volume-discharge.
+in Figure 1 and Figure 2.
-{#fig-volume-depth}
+alt="Figure 1: Distribution of water depths over the primary, secondary, and tertiary system."
-{#fig-volume-discharge}
+alt="Figure 2: Discharge as a function of basin storage volume."
An example of the resulting parameters for a single cell is shown in
@tbl-hupsel-v-h. The first row shows the water levels when the basin is empty.
@@ -85,21 +85,21 @@ in the test cases for the sake of simplicity.
: Volume-level table for a single cell in the Hupsel basin. {#tbl-hupsel-v-h}
-@fig-grid-volume shows the volume of the first row of the cell based
+Figure 3 shows the volume of the first row of the cell based
input for the primary system. Symbology is set to unique values. While water
levels differ per cell in this parametrization, the "normative volume" defined
above is shared by all cells in a basin.
-{#fig-grid-volume}
+alt="Figure 3: Basin normative volume of the primary system."
-@fig-grid-volume shows the water level corresponding to the normative storage
+Figure 3 shows the water level corresponding to the normative storage
volume based input for the primary system (it corresponds to the value shown in
the first row of the primary column in @tbl-hupsel-v-h). We see a clear gradient
from west to east: as our simplified parametrization assumes a constant water
depth for all cells in a single system, water levels spatially fall and rise
with the bottom elevation.
-{#fig-grid-level}
+alt="Figure 4: Water level corresponding to the normative basin volume of the primary system."
## Example: Configuration
@@ -169,9 +169,9 @@ From these tests, we expect the following behavior:
2. In case of negative recharge (evapotranspiration), infiltration occurs in the surface waters.
Infiltration should be zero when the basin volume is 0.
-{#fig-hupsel-gwb-steady}
+alt="Figure 5: Water balance of the MODFLOW 6 boundary conditions for the Hupsel basin for a standalone MODFLOW 6 run. The four sequential steady states (01, 02, 03, 04) use net groundwater recharge values of 1.0, 0.5, -0.05 and -0.1 mm/d."
-@fig-hupsel-gwb-steady shows the water balance of steady-state for submodel of the
+Figure 5 shows the water balance of steady-state for submodel of the
LHM that has been by selecting the cells belonging to the district containing
the Hupsel catch, the Berkel.
@@ -194,9 +194,9 @@ steady-state model), the surface waters provide mostly inflow, and recharge is
a negative term. In this case, the secondary system provides a small amount of
infiltration; most of the water is drawn from the surroundings instead.
-{#fig-hupsel-gwb-volume0}
+alt="Figure 6: Water balance of the MODFLOW 6 boundary conditions for the Hupsel basin for a zero volume run. The four sequential steady states (01, 02, 03, 04) use net groundwater recharge values of 1.0, 0.5, -0.05 and -0.1 mm/d."
-@fig-hupsel-gwb-volume0 shows the same model, with 0-basin volume which causes
+Figure 6 shows the same model, with 0-basin volume which causes
water levels to be set equal to bed elevation. Consequently, primary and
secondary outflow terms are larger for positive groundwater recharge as they
drain at a lower level and intercept the water before the tertiary system does.
@@ -215,9 +215,9 @@ process, but without Ribasim. This results in volumes of 0.0, so all MODFLOW 6
water levels are set equal to bed elevation.
3. A coupled run where the water levels are updated by Ribasim.
-{#fig-tol-gwb-steady}
+alt="Figure 7: Water balance of the MODFLOW 6 boundary conditions for De Tol basin for a standalone MODFLOW 6 run. The four sequential steady states (01, 02, 03, 04) use net groundwater recharge values of 1.0, 0.5, -0.05 and -0.1 mm/d."
-@fig-tol-gwb-steady shows the water balance of steady-state for a submodel of
+Figure 7 shows the water balance of steady-state for a submodel of
the LHM for the Polder de Tol and its surroundings. While groundwater recharge
is the dominant ingoing flow, lateral groundwater flow (over the entire depth
of the groundwater model) is a sizable inflow for the area; the larger lateral
@@ -228,9 +228,9 @@ ephemeral tertiary ditches, but by the permanently water-bearing ditches of the
primary and secondary system. Unlike the Hupsel, the water balance does not
shrink to very small discharges, as there is sizable regional groundwater flow.
-{#fig-tol-gwb-volume0}
+alt="Figure 8: Water balance of the MODFLOW 6 boundary conditions for De Tol basin for a zero volume run. The four sequential steady states (01, 02, 03, 04) use net groundwater recharge values of 1.0, 0.5, -0.05 and -0.1 mm/d."
-@fig-tol-gwb-volume0 shows the same model, with 0-basin volume which causes
+Figure 8 shows the same model, with 0-basin volume which causes
water levels to be set equal to bed elevation. The total discharge is larger:
the primary and secondary systems are set to lower levels, and so the head
difference is larger. While De Tol's evapotranspiration excess can be fed by
diff --git a/docs/index.qmd b/docs/index.qmd
index 8597df9a7..fb3b7434e 100644
--- a/docs/index.qmd
+++ b/docs/index.qmd
@@ -52,9 +52,9 @@ a set of symbolic equations, and can be connected to each other. From this a sim
system of equations is generated automatically. We use solvers with adaptive time stepping
from [DifferentialEquations.jl](https://diffeq.sciml.ai/stable/) to get results.
-
+alt="Figure 1: Example timeseries of a single basin, the Hupselse Beek, with the input and output fluxes on the top, and the storage volume (the state) below."
-
+alt="Figure 2: Example bar plot of the daily waterbalance for the Hupselse Beek, comparing results of Mozart (left) and Ribasim (right)."
# Introduction
@@ -108,7 +108,7 @@ Local Surface Water (LSW)). Each basin has an associated polygon, and the set of
connected to each other as described by a graph, which we call the network. Below is a
representation of both on the map.
-
+alt="Figure 3: Mozart Local Surface Water polygons and their drainage."
The network is described as graph. Flow can be bi-directional, and the graph does not have
to be acyclic.
@@ -131,16 +131,16 @@ lines between A, B, C, and D) instead.
Multiple basins may exist within the same spatial polygon, representing
different aspects of the surface water system (perennial ditches, ephemeral
-ditches, or even surface ponding). @fig-p, @fig-s, @fig-t show the 25.0 m
+ditches, or even surface ponding). Figure 4, Figure 5, Figure 6 show the 25.0 m
rasterized primary, secondary, and tertiary surface waters as identified by BRT
TOP10NL [@pdoktopnl] in the Hupsel basin (as defined in the Mozart LSW's).
These systems may represented in multiple ways.
-{#fig-p}
+alt="Figure 4: Hupsel: primary surface water."
-{#fig-s}
+alt="Figure 5: Hupsel: secondary surface water."
-{#fig-t}
+alt="Figure 6: Hupsel: tertiary surface water."
As a single basin (A) containing all surface water, discharging to its
downstream basin to the west (B):
diff --git a/docs/qgis/index.qmd b/docs/qgis/index.qmd
index ee1d4a7f7..f4ebc1aad 100644
--- a/docs/qgis/index.qmd
+++ b/docs/qgis/index.qmd
@@ -14,18 +14,18 @@ Download `ribasim_qgis.zip`, see the [download section](../index.qmd#sec-downloa
Plugins menu > Manage and Install Plugins...
-{fig-align="left"}
+
Select "Install from ZIP":
- Browse to the `ribasim_qgis.zip` file containing the plugin that was downloaded earlier
- Click "Install Plugin"
-{fig-align="left"}
+
Start the Ribasim plugin.
-{fig-align="left"}
+
### Add iMOD plugin
@@ -41,20 +41,20 @@ Documentation on the Time Series widget can be found in the [iMOD documentation]
Open example model `database.gpkg` or create a new model.
-{fig-align="left"}
+
Check if your coordinate reference system (CRS) is set correctly.
-{fig-align="left"}
+
If you are working with an unknown CRS, right click the model database group in Layers,
and click "Set Group CRS...".
-{fig-align="left"}
+
If you are modeling the Netherlands, select "Amersfoort / RD New" (EPSG:28992).
-{fig-align="left"}
+
## Edit nodes
@@ -62,35 +62,35 @@ If you are modeling the Netherlands, select "Amersfoort / RD New" (EPSG:28992).
Select the Node layer.
-{fig-align="left"}
+
Turn on the edit mode to be able to add nodes on the map.
-{fig-align="left"}
+
Add nodes to the map with a left click and select the node type.
-{fig-align="left"}
+
Turn the edit mode off and save the edits to the Nodes layer.
-{fig-align="left"}
+
### Edit tables
Right click a layer and select "Open Attribute Table".
-{fig-align="left"}
+
Click the yellow pencil icon on the top left to enable editing, and copy and paste a record.
A record can be selected by clicking on the row number.
-{fig-align="left"}
+
Adjust the content. If you prefer, it also works to copy data with the same columns from
Excel. Turn off edit mode and save changes to the layer.
-{fig-align="left"}
+
## Connect nodes
@@ -99,21 +99,21 @@ Excel. Turn off edit mode and save changes to the layer.
Make sure the Snapping Toolbar is visible, by going to the View > Toolbars menu. Turn on
snapping mode by clicking the magnet and set the snapping distance to 25 pixels.
-{fig-align="left"}
+
### Create connecting edges
Select the Edge layer and turn on the edit mode.
-{fig-align="left"}
+
Select "Add line feature".
-{fig-align="left"}
+
Create a connection by left clicking a source node and right clicking the destination node.
-{fig-align="left"}
+
Now leave the edit mode and save the results to the layer.
@@ -140,15 +140,15 @@ database = "database.gpkg"
In QGIS select the model group.
-{fig-align="left"}
+
In the Ribasim plugin widget, select the Results tab and click "Associate Results".
-{fig-align="left"}
+
Select `results/basin.arrow`.
-{fig-align="left"}
+
This adds metadata to the model that the iMOD plugin can use to find the timeseries data
that is associated to the model nodes.
@@ -157,17 +157,17 @@ that is associated to the model nodes.
Click the "Time Series" button of the iMOD plugin.
-{fig-align="left"}
+
Select the variables that you want to plot.
-{fig-align="left"}
+
Click "Select points" and select a node by dragging a rectangle around it on the map.
Hold the Ctrl key to select multiple nodes. Currently only the Basin nodes can be plotted.
-{fig-align="left"}
+
The associated time series are shown the the graph.
-{fig-align="left"}
+
