Nothing
\donttest{
data("slb1_meteo")
data("slb1_soil")
opts <- set_optionsLWFB90(budburst_method = "Menzel", enddate = as.Date("2002-12-31"))
# define parameter sets
param_l <- list(spruce = set_paramLWFB90(maxlai = 5,
budburst_species = "Picea abies (frueh)",
winlaifrac = 0.8),
beech = set_paramLWFB90(maxlai = 6,
budburst_species = "Fagus sylvatica",
winlaifrac = 0))
soil <- cbind(slb1_soil, hydpar_wessolek_tab(slb1_soil$texture))
# define list of soil objects
soils <- list(soil1 = soil, soil2 = soil)
# define list of climate objects
climates <- list(clim1 = slb1_meteo, clim2 = slb1_meteo)
# run two parameter sets on a series of climate and soil-objects
res <- run_multisite_LWFB90(param_b90 = param_l,
options_b90 = opts,
soil = soils,
climate = climates)
names(res)
# set up and run individual parameter sets for individual locations
# set up location parameters
loc_parm <- data.frame(loc_id = names(climates),
coords_y = c(48.0, 54.0),
eslope = c(30,0),
aspect = c(180,0))
# create input list of multiple param_b90 list objects
param_l <- lapply(names(climates), function(x, loc_parms) {
parms <- set_paramLWFB90()
parms[match(names(loc_parm),names(parms), nomatch = 0)] <-
loc_parm[loc_parm$loc_id == x, which(names(loc_parm) %in% names(parms))]
parms
}, loc_parm = loc_parm)
names(param_l) <- c("locpar1", "locpar2")
res <- run_multisite_LWFB90(param_b90 = param_l,
options_b90 = opts,
soil = soils,
climate = climates)
names(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.