This function calculates (harmonic-) means of the soilphys dataframe per layer for the flux calculation.
Arguments
- soilphys
(dataframe) the soilphys dataframe
- layers_map
(dataframe) containing the following parameters:
- upper
upper limit of the layer in cm
- lower
lower limit of the layer in cm
- layer
the name of the layer
- ...
any
id_cols
necessary
- param
(character vector) A vector containing the names of the variables in soilphys to be carried over
- funs
(character vector) A vector defining the type of mean to be used. One of "arith" or "harm"
- id_cols
(character vector) A list of all columns, that, together, uniquely identify one profile.
Value
soilphys, reduced to any id_cols
and param
columns and
summarised to the layers provided in layers_map
#' @examples data(soildiff) layers_map <- data.frame(site = rep(c("site_a", "site_b"),each = 2), upper = c(5,0,7,0), lower = c(0,-100,0,-100), layer = rep(c("HU","M1"),times = 2) )soildiff$depth <- (soildiff$upper + soildiff$lower) / 2soilphys_layered(soilphys = soildiff, layers_map = layers_map, param = c("TPS","a"), funs = c("harm","arith"), id_cols = c("site"))
See also
Other soilphys:
check_soilphys()
,
complete_soilphys()
,
discretize_depth()