fnc_read_climdb: Reading from Climate Database

In rhabel/modLWFB90: R-package for applying LWF-Brook90 directly from FVA-BW servers

Reading from Climate Database

Description

This function reads in the Project Climate Data files in "R:/klima/whh/brook90_input/rds/" plus the respective tranche, and processes the data correctly, so that it can be used via the function parameters climate and climate_args in run_multisite_LWFB90. Function arguments must be passed to run_multisite_LWFB90 as a list of lists, which is created here. See example on how to do this.

Usage

fnc_read_climdb(IDs, clim_dir, id_standard, mindate = NA, maxdate = NA)

Arguments

clim_dir	path to project climate data file directory
id_standard	climate-cell where point is located, use fnc_relateCoords to create this information, see example
mindate	R-date-object, optional. Set if only a selection of the data (1961-2018) in the climate-database shall be retrieved.
maxdate	R-date-object, optional, same as mindate. Must be provided with mindate if sub-selection on date is desired.
ID_custom	ID-name of point in project climate data base

Value

returns climate data in correct format

Examples

# set options
options.b90 <- set_optionsLWFB90(startdate = as.Date("2009-01-01"),
                                 enddate = as.Date("2010-12-31"),
                                 root_method = "soilvar")

# set clim_args
df.ids <- left_join(test_ids_bze, fnc_relateCoords(test_ids_bze)[c("ID_custom", "id_standard", "tranche")])

# climate arguments for on-the-fly-processing
climlist <- split(df.ids, seq(nrow(df.ids)))
names(climlist) <- as.character(df.ids$ID_custom)

clim_args <-
  lapply(climlist,
         function(x) list(IDs = as.character(x$ID_custom),
                          id_standard = x$id_standard,
                          clim_dir = paste0(input_clim, "tr", x$tranche, "/"),
                          mindate = options.b90$startdate,
                          maxdate = options.b90$enddate))

# create soil
soil.test <- fnc_get_soil(df.ids = test_ids_bze,
                          soil_option = "BZE",
                          PTF_to_use = "WESSOLEK")
# create parameters
parms.test <- fnc_get_params(df.ids = test_ids_bze,
                             tree_species = "spruce")

# remove poins without soil information
with_soil_data <- which(df.ids$ID_custom %in% names(soil.test))
parms.test <- parms.test[with_soil_data]
clim_args <- clim_args[with_soil_data]

res <- run_multisite_LWFB90(options_b90 = options.b90,
                            param_b90 = parms.test,
                            soil = soil.test,

                            climate = fnc_read_climdb,
                            climate_args = clim_args,

                            rtrn_input = F,

                            cores = parallel::detectCores())

# single df
clim_args[[1]]
test <- fnc_read_climdb(IDs = clim_args[[1]]$IDs,
                        id_standard = clim_args[[1]]$id_standard,
                        clim_dir = clim_args[[1]]$clim_dir,
                        mindate = clim_args[[1]]$mindate,
                        maxdate = clim_args[[1]]$maxdate)

rhabel/modLWFB90 documentation built on Nov. 21, 2024, 3:28 a.m.