## Time-stamp: "Last modified 2023-04-24 15:12:02 mluebke"
init_cell = as.data.frame(init_cell, check.names = FALSE),
##         Diffusion parameters and boundary conditions         ##
## list of boundary conditions/inner nodes
init = as.data.frame(init_diffu, check.names = FALSE),
##                  Chemistry module (Phreeqc)                 ##
dht_species <- c("H", "O", "Charge", "C(4)", "Ca", "Cl", "Mg", "Calcite",
##              Putting all those things together              ##