R/beta.R

In aj-sykes92/soilc.ipcc: The IPCC (2019) Tier 2 Steady-State Method for Soil Carbon

#' @title Calculate \code{beta} parameter
#' @description Calculates the \code{beta} model parameter, which is derived from carbon inputs,
#' nitrogen and lignin fractions.
#' @param c_input Numeric vector indicating combined carbon inputs for the model, in tonnes C ha^-1.
#' @param lignin_frac Numeric vector indicating lignin fraction of the carbon inputs, fractional.
#' @param n_frac Numeric vector indicating nitrogen fraction of the carbon inputs, fractional.
#' @param params A named list containing model parameters. Defaults to the package default parameter
#' set \code{soilc_params}.
beta <- function(c_input, lignin_frac, n_frac, params = soilc_params) {
  c_input * (params$sp1$be - params$sp2$be * lignin_frac / n_frac)
}