PELMO_runs: Set up, execute and evaluate FOCUS PELMO runs

In jranke/rPELMO: Utilities for working with FOCUS PELMO

Description

Per default, the runs are not only set up but also executed with FOCUS PELMO, the results are processed and returned. Currently, only FOCUS PELMO as installed on Linux (or other Unix systems) using the install_PELMO from the PELMO.installeR package maintained on github is supported. In such installations, FOCUS PELMO is installed into the package installation directory of PELMO.installeR and run using wine.

Usage

PELMO_runs(runs, psm_dir = ".", version = "5.5.3",
  PELMO_base = "auto", execute = TRUE, cores = getOption("mc.cores",
  2L), evaluate = TRUE, overwrite = FALSE)

run_PELMO(runs, version = "5.5.3", PELMO_base = "auto",
  cores = getOption("mc.cores", 2L))

evaluate_PELMO(runs, version = "5.5.3", PELMO_base = "auto")

Arguments

runs	A list of lists. Each inner lists has an element named 'psm' that holds the psm string, and elements named using three letter crop acronyms, as used in FOCUS_PELMO_crops, that hold character vectors of three letter scenario acronyms as used in FOCUS_GW_scenarios_2012.
psm_dir	The directory where the psm files are located
version	The FOCUS PELMO version
PELMO_base	Where the FOCUS PELMO installation is located
execute	Should PELMO be executed directly?
cores	The number of cores to execute PELMO runs in parallel
evaluate	Should the results be returned?
overwrite	Should existing run directories be overwritten?

Details

As a side effect, an R data file (period_rPELMO.rda) is generated in each run directory, holding the results for all FOCUS periods, equivalent to the period.plm file generated by the FOCUS PELMO GUI.

Value

If evaluate is TRUE, a list of lists of matrices holding the PEC data.

References

PELMO.installeR https://pkgdown.jrwb.de/PELMO.installeR

Wine https://winehq.org

PELMO test results http://esdac.jrc.ec.europa.eu/public_path/projects_data/focus/gw/models/pelmo/test-results/test_results_FOCUS_PELMO_5_5_3.doc

Examples

# Reproduce the official test results for annual application of Pesticide D
# to winter cereals at the day before emergence
runs_1 <- list(
  list(psm = 'Pesticide_D',
       win = c("Cha", "Ham", "Jok", "Kre", "Oke", "Pia", "Por", "Sev", "Thi")),
  list(psm = 'Pesticide_D_1_day_pre_em_every_third_year',
       pot = c("Cha", "Ham", "Jok", "Kre", "Oke", "Pia", "Por", "Sev", "Thi")))
time_1 <- system.time(
  PECgw_1 <- PELMO_runs(runs_1, psm_dir = system.file("testdata", package = "rPELMO"),
    cores = 15, overwrite = TRUE)
)
print(PECgw_1)
# We get exactly the same PECgw values (on Linux, calling PELMO using Wine).
print(time_1)
if(!inherits(try(cpuinfo <- readLines("/proc/cpuinfo")), "try-error")) {
  cat(gsub("model name\t: ", "CPU model: ", cpuinfo[grep("model name", cpuinfo)[1]]))
}

# Demonstrate some results with metabolites.
runs_2 <- list(list(psm = 'Pesticide_D_1_May_every_other_year_mets',
                    win = c("Cha", "Ham", "Kre")))
PECgw_2 <- PELMO_runs(runs_2, psm_dir = system.file("testdata", package = "rPELMO"),
  cores = 3, overwrite = TRUE)
print(PECgw_2)

jranke/rPELMO documentation built on May 26, 2019, 3:33 p.m.