Yearly maximum daily flow for a return period

The maximum daily flow is calculated from an adjustement on a Gumbel law.

Functions named calcQJXA_[return_period] are also defined for commun indicators (2, 10, 20 years).

Usage

calcQJXAn(x, return_period)

# S3 method for QJXA
calcQJXAn(x, return_period)

# S3 method for data.frame
calcQJXAn(x, return_period)

# S3 method for GumbelParams
calcQJXAn(x, return_period)

calcQJXA2(x)

calcQJXA10(x)

calcQJXA20(x)

Arguments

x: data.frame of flows (the first column should be a POSIXct), or a QJXA object (See calcQJXA)
return_period: numeric Return period. Its unit is given by the aggregation time step used for the Gumbel law.

Value

A numeric with the calculated flow

Examples

#! load data
data(L0123001, package = "airGR")

#! define return period
return_period <- 10

# Calculation in detail
QJXA <- calcQJXA(x = BasinObs$DatesR, flows = BasinObs$Qls / 1000)
gumbel <- getGumbelParams(QJXA)
calcQJXAn(gumbel, return_period)
#> [1] 75.19146

# Or directly
dfQ <- data.frame(DatesR = BasinObs$DatesR, Q = BasinObs$Qls / 1000)
calcQJXAn(dfQ, return_period)
#>        Q 
#> 73.43942