data("slb1_soil")
data("slb1_meteo")
opts <- set_optionsLWFB90(startdate = as.Date("2002-06-01"), enddate = as.Date("2002-06-05"))
parms <- set_paramLWFB90()
soil <- cbind(slb1_soil, hydpar_wessolek_tab(texture = slb1_soil$texture))
outsel <- set_outputLWFB90()
outsel[,] <- 1L
res <- run_LWFB90(options_b90 = opts,
param_b90 = parms,
climate = slb1_meteo,
soil = soil)
# Calculate output-aggregations using the returned object
process_outputs_LWFB90(res, selection = outsel)
# or calculate aggregations at run time by passing the function via output_fun-arg
run_LWFB90(options_b90 = opts,
param_b90 = parms,
climate = slb1_meteo,
soil = soil,
rtrn_input = FALSE,
output_fun = process_outputs_LWFB90,
selection = outsel)$output_fun
Run the code above in your browser using DataLab