MultiRun.B90: Make a multirun simulation using a set of variable input...

In pschmidtwalter/brook90r: Run the LWF-Brook90 hydrological model from within R.

Description

Repeatedly calls Run.B90 with varying input parameters in parallel, stores the single run results in subdirectories and returns the results as a named list.

Usage

MultiRun.B90(nRuns, param_var, param_const, options.b90, soil, climate,
  multirun.dir = "MultiRuns", keep.subdirs = FALSE,
  singlerun_names = paste0("RunNo", 1:nRuns), cores = 3,
  showProgress = TRUE, ...)

Arguments

nRuns	Number of single runs
param_var	data.frame of variable input parameters with realisations for each single run.
param_const	Named list of constant model input parameters to use in all single runs.
options.b90	Named list of constant model control options passed to Run.B90
soil	data.frame with soil properties passed to Run.B90
climate	data.frame with climate data passed to Run.B90
multirun.dir	The directory where to create the subdirectories for the the single runs.
singlerun_names	character vector with the names of the single runs. Used for creating the subdirectories and for naming the returned list of singlerun results.
cores	Number of cores to use for parallel processing
showProgress	Show progressbar?
...	Further arguments passed to Run.B90.

Value

A named list with the results for the single runs. Simulation errors are passed on.

Examples

#Set up lists containing model control options and model parameters:
param.b90 <- MakeParam.B90()
options.b90 <- MakeOptions.B90()

# Derive soil hydraulic properties from soil physical properties using pedotransfer functions
soil <- cbind(soil_slb1, hydpar_wessolek_mvg(soil_slb1$texture))

# Set start and end Dates for the simulation
options.b90$startdate <- as.Date("2000-01-01")
options.b90$enddate <- as.Date("2004-12-31")

# choose the 'Coupmodel' shape option for the annual lai dynamic, with fixed budburst and leaf fall dates:
options.b90$annuallaidyn <- 'Coupmodel'
options.b90$budburst <- 'fixed'
options.b90$leaffall <- 'fixed'

#set up data.frame with variable leaf area index parameters
n <- 10
vary_parms <- data.frame(maxlai = runif(n,4,8),
                         shapestart = runif(n, 0.1,1),
                         winlaifrac = runif(n, 0,0.5),
                         budburstdoy = runif(n,100,150),
                         optdoy =runif(n,180,240))

# Make a Multirun-Simulation
b90.multi <- MultiRun.B90(nRuns = n,
                     param_var = vary_parms,
                     param_const = param.b90,
                     options.b90 = options.b90,
                     soil = soil,
                     climate = meteo_slb1,
                     multirun.dir = "MultiRuns",
                     keep.subdirs = FALSE,
                     singlerun_names = paste0("result.",1:n),
                     cores = 2,
                     output.log = F,
                     output.param.options = F,
                     path_b90.exe = "b90.exe")

pschmidtwalter/brook90r documentation built on April 6, 2020, 6:35 p.m.