Calculate the curve(s) of the Bagrov relation

Calculate the curve(s) of the Bagrov relation

Usage

calculate_bagrov_curve(
  effectivity,
  P_over_Ep_max = 4,
  Ep = 650,
  delta_Ea = 0.1,
  dbg = FALSE
)

Arguments

effectivity

Bagrov effectivity value (= n-value)

P_over_Ep_max

maximum value of the ratio P/Ep with precipitation P and potential evaporation Ep.

Ep

potential evaporation in mm/a (default: 650)

delta_Ea

step of actual evaporation (default: 0.1)

dbg

logical indicating whether or not to show debug messages

Value

data frame with columns P_over_Ep (P/Ep), Ea_over_Ep

(Ea/Ep), effectivity (n-value)

Examples

if (FALSE) {
bagrov_data <- calculate_bagrov_curve(1:2)
plot(bagrov_data$P_over_Ep, bagrov_data$Ea_over_Ep)
}