max/poet
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/poet.git synced 2025-12-16 04:48:23 +01:00
Code Issues Packages Projects Releases Wiki Activity

Update grid and diffusion setup in test.R

Browse Source
This commit is contained in:
Max Luebke 2024-03-19 22:03:22 +00:00
parent 793c2bb4ff
commit 8cd69bcde4
1 changed files with 19 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
bench/barite/het/test.R
View File

@ -2,14 +2,30 @@ grid_def <- matrix(c(2, 3), nrow = 2, ncol = 5)
# Define grid configuration for POET model # Define grid configuration for POET model
grid_setup <- list( grid_setup <- list(
pqc_in = "./barite_het.pqi", pqc_in_file = "./barite_het.pqi",
pqc_db = "./db_barite.dat", # Path to the database file for Phreeqc 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_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 grid_size = c(ncol(grid_def), nrow(grid_def)), # Size of the grid in meters
constant_cells = c() # IDs of cells with constant concentration 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 # Define a setup list for simulation configuration
setup <- list( 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
) )