Rackliff Index

Rackliff Index

Usage

Rackliff_Index(ta900, hus900, ta500, nthreads = 1)

Arguments

ta900

air temperature [lon, lat, time] at 900 hPa (K).

hus900

specific humidity [lon, lat, time] at 900 hPa (kg/kg).

ta500

air temperature [lon, lat, time] at 500 hPa (K).

nthreads

an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

Value

the computed Rackliff Index (K).

References

Haklander, A. J., & Van Delden, A. (2003). Thunderstorm predictors and their forecast skill for the Netherlands. Atmospheric Research, 67–68, 273–299. doi:10.1016/S0169-8095(03)00056-5 .

Examples

data("ERA5_pl")
pa <- ERA5_pl$plev / 100 # convert to hPa
hus <- ERA5_pl$hus
ta <- ERA5_pl$ta

l500 <- which(pa == 500)
l900 <- which(pa == 900)
res <- Rackliff_Index(
  ta[, , l900, ],
  hus[, , l900, ],
  ta[, , l500, ]
)