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

feat: update input scripts to be used with DHT

Browse Source
This commit is contained in:
Max Lübke 2022-11-03 12:26:19 +01:00
parent d70073724e
commit aa05eaad76
2 changed files with 10 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
data/SimDol1D_diffu.R
View File

@ -146,7 +146,10 @@ selout <- c(
"-kinetic_reactants Calcite Dolomite", "-equilibrium O2g"
)
# TODO: dt and iterations
# Needed when using DHT
signif_vector <- c(7, 7, 7, 7, 7, 7, 7, 5, 5)
prop_type <- c("act", "act", "act", "act", "logact", "logact", "ignore", "act", "act")
prop <- names(init_cell)
iterations <- 500
@ -158,7 +161,7 @@ setup <- list(
# Cf = 1,
grid = grid,
diffusion = diffusion,
prop = names(init_cell),
prop = prop,
immobile = c(7, 8, 9),
kin = c(8, 9),
ann = list(O2g = -0.1675),
@ -175,7 +178,3 @@ setup <- list(
iterations = iterations,
timesteps = rep(10, iterations)
)
# not needed yet
# signif_vector <- c(7, 7, 7, 7, 7, 7, 7, 5, 5)
# prop_type <- c("act", "act", "act", "act", "logact", "logact", "ignore", "act", "act")

11
data/SimDol2D_diffu.R
View File

@ -154,7 +154,10 @@ selout <- c(
"-kinetic_reactants Calcite Dolomite", "-equilibrium O2g"
)
# TODO: dt and iterations
# Needed when using DHT
signif_vector <- c(7, 7, 7, 7, 7, 7, 7, 5, 5)
prop_type <- c("act", "act", "act", "act", "logact", "logact", "ignore", "act", "act")
prop <- names(init_cell)
iterations <- 10
@ -166,7 +169,7 @@ setup <- list(
# Cf = 1,
grid = grid,
diffusion = diffusion,
prop = names(init_cell),
prop = prop,
immobile = c(7, 8, 9),
kin = c(8, 9),
ann = list(O2g = -0.1675),
@ -183,7 +186,3 @@ setup <- list(
iterations = iterations,
timesteps = rep(10, iterations)
)
# not needed yet
# signif_vector <- c(7, 7, 7, 7, 7, 7, 7, 5, 5)
# prop_type <- c("act", "act", "act", "act", "logact", "logact", "ignore", "act", "act")