Add inner boundaries for diffusion

Max Luebke 2024-04-04 14:49:40 +02:00
parent ade856e395
commit fdf6e6cb44

@ -54,6 +54,11 @@ Diffusion <- list(
"E" = list(...),
# You can also leave sides out -> closed boundary on complete side
),
inner_boundaries = list( # define inner constant boundaries for diffusion
"row" = numeric(), # Define the row ID; counting begins @ 1
"col" = numeric(), # Define the column ID; counting begins @ 1
"sol_id" = numeric() # Define the Phreeqc solution ID for the boundary
),
# alphas can be set by a list of matrices for ALL species to transport
alpha_x = list(
"element1" = numeric(nrow * ncol), # column-major order