From 8969fea7e2fc2b1f87f18b7c7705a7bf6d59badb Mon Sep 17 00:00:00 2001 From: Max Luebke Date: Tue, 19 Mar 2024 22:03:22 +0000 Subject: [PATCH] Update grid and diffusion setup in test.R --- bench/barite/het/test.R | 22 +++++++++++++++++++--- 1 file changed, 19 insertions(+), 3 deletions(-) diff --git a/bench/barite/het/test.R b/bench/barite/het/test.R index 97c642b51..f7269e3f7 100644 --- a/bench/barite/het/test.R +++ b/bench/barite/het/test.R @@ -2,14 +2,30 @@ grid_def <- matrix(c(2, 3), nrow = 2, ncol = 5) # Define grid configuration for POET model grid_setup <- list( - pqc_in = "./barite_het.pqi", - pqc_db = "./db_barite.dat", # Path to the database file for Phreeqc + pqc_in_file = "./barite_het.pqi", + pqc_db_file = "./db_barite.dat", # Path to the database file for Phreeqc grid_def = grid_def, # Definition of the grid, containing IDs according to the Phreeqc input script grid_size = c(ncol(grid_def), nrow(grid_def)), # Size of the grid in meters constant_cells = c() # IDs of cells with constant concentration ) +diffusion_setup <- list( + boundaries = list( + "E" = list( + "type" = rep("constant", nrow(grid_def)), + "sol_id" = rep(4, nrow(grid_def)), + "cell" = seq_len(nrow(grid_def)) + ) + ), + alpha_x = 1e-6, + alpha_y = matrix(runif(10, 1e-8, 1e-7), + nrow = nrow(grid_def), + ncol = ncol(grid_def) + ) +) + # Define a setup list for simulation configuration setup <- list( - grid = grid_setup # Parameters related to the grid structure + Grid = grid_setup, # Parameters related to the grid structure + Diffusion = diffusion_setup # Parameters related to the diffusion process )