data-raw/make_old_new.R

# Set location where to write the comparison data out to.
BASE_DIR <- here::here()
TEST_DIR <- file.path(BASE_DIR, "tests", "testthat")

time <- floor(seq(0, 36500, length.out = 25))

# Run the different outputs.
mend_out <- memc_solve(mod = MEND_config, time)
names(mend_out) <-
  c("time", "variable", "old_value", "units", "name")

comisssion_out <- memc_solve(mod = COMISSION_config, time)
names(comisssion_out) <-
  c("time", "variable", "old_value", "units", "name")

corpse_out <- memc_solve(CORPSE_config, time)
names(corpse_out) <-
  c("time", "variable", "old_value", "units", "name")

corpse_out <- memc_solve(CORPSE_config, time)
names(corpse_out) <-
  c("time", "variable", "old_value", "units", "name")

mimcs_out <- memc_solve(MIMCS_config, time)
names(mimcs_out) <-
  c("time", "variable", "old_value", "units", "name")

mems_out <- memc_solve(MEMS_config, time)
names(mems_out) <-
  c("time", "variable", "old_value", "units", "name")

bams_out <- memc_solve(BAMS_config, time)
names(bams_out) <-
  c("time", "variable", "old_value", "units", "name")

out <- rbind(mend_out, comisssion_out, corpse_out, mimcs_out, mems_out, bams_out)

write.csv(out,
          file = file.path(TEST_DIR, "old-new.csv"),
          row.names = FALSE)
Microbial-Explicit-Model/MEMC documentation built on April 12, 2025, 12:50 p.m.