ETCCDI canonical index R99p. Annual total precipitation on days where
daily precipitation exceeds the 99th percentile of wet-day precipitation
in a reference period (default 1961 to 1990). Same convention as
ck_r95p().
Arguments
- precip
Numeric vector of daily precipitation (mm).
- dates
Date vector of the same length as
precip. Must contain data covering the reference period.- ref_start, ref_end
Integer. Reference period boundary years (inclusive). Defaults to 1961 and 1990.
- period
Character. Aggregation period:
"annual"(default) or"monthly".
Examples
set.seed(1)
dates <- seq(as.Date("1961-01-01"), as.Date("1991-12-31"), by = "day")
precip <- pmax(rgamma(length(dates), shape = 0.4, scale = 8) - 1, 0)
tail(ck_r99p(precip, dates))
#> period value index unit
#> 26 1986-01-01 0.00000 r99p mm
#> 27 1987-01-01 79.45814 r99p mm
#> 28 1988-01-01 62.00729 r99p mm
#> 29 1989-01-01 68.95095 r99p mm
#> 30 1990-01-01 37.26143 r99p mm
#> 31 1991-01-01 0.00000 r99p mm