Update diffusion boundary definitions and variable names

Max Lübke 2024-03-19 15:09:04 +01:00
parent f63f7be7dc
commit 9c6e1bbc14

@ -50,23 +50,23 @@ Diffusion <- list(
"cell" = numeric(), # Define the cell ID; counting begins @ 1
"sol_id" = numeric(), # Define the Phreeqc solution ID for the boundary
),
# Only those boundaries needs to be defined, which are constant
# Only those boundaries needs to be defined, which are not "closed"
"E" = list(...),
# You can also leave sides out -> closed boundary on complete side
),
# alphas can be set by a list of matrices for ALL species to transport
alphax = list(
alpha_x = list(
"element1" = numeric(nrow * ncol), # column-major order
# OR
"element2" = numeric(1), # for homogenous diffusion coefficients
...
),
# or by defining a single matrix for all species, assuming homogenous diffusion coefficients among species at the same cell index
alphax = numeric(nrow * ncol), # column-major order
alpha_x = numeric(nrow * ncol), # column-major order
# or by defining a single value, which acts as a homogeneous diffusion coefficient for all species and cells
alphax = numeric(1),
alpha_x = numeric(1),
# ----
alphay = alphax # See documentation for alphax
alpha_y = alpha_x # See documentation for alphax
)
```