t90p: Calculate T90p Index

In rIACI: Iberian Actuarial Climate Index Calculations

Calculate T90p Index

Description

Calculates the combined percentage of days when temperature is above the 90th percentile.

Usage

t90p(ci, freq = c("monthly", "annual"))

Arguments

ci	List. Climate input object.
freq	Character. Frequency of calculation, either "monthly" or "annual".

Value

Data frame with dates and calculated T90p values.

Examples

# 1. Generate a daily date sequence from 1960-01-01 to 2020-12-31
dates <- seq.Date(
  from = as.Date("1960-01-01"),
  to   = as.Date("2020-12-31"),
  by   = "day"
)

# 2. Create random weather data for each date
n <- length(dates)
tmax <- runif(n, min = 5, max = 40)
tmin <- runif(n, min = -10,  max = 5)
# Example: use a Poisson distribution to simulate precipitation
prec <- rpois(n, lambda = 2)
# Random wind speeds, e.g., 0 to 10 m/s
wind <- runif(n, min = 0, max = 10)

# 3. Construct the climate_input object
ci <- climate_input(
  tmax  = tmax,
  tmin  = tmin,
  prec  = prec,
  wind  = wind,
  dates = dates
)
# Then:
# 4. Calculate monthly T90p index
t90p_values <- t90p(ci, freq = "monthly")

rIACI documentation built on April 12, 2025, 9:16 a.m.