Humidity Index

Humidity Index

Usage

Humidity_index(ta850, hus850, ta700, hus700, ta500, hus500, nthreads = 1)

Arguments

ta850

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

hus850

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

ta700

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

hus700

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

ta500

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

hus500

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

nthreads

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

Value

the computed Humidity Index (K).

References

Litynska, Z., J. Parfiniewicz, and H. Pinkowski, 1976: The prediction of airmass thunderstorms and hails. W.M.O. Bull.,450, 128–130.

Examples

data("ERA5_pl")
pa <- ERA5_pl$plev / 100 # convert to hPa
hus <- ERA5_pl$hus
ta <- ERA5_pl$ta
z <- ERA5_pl$z / 9.80665 # convert to m
l500 <- which(pa == 500)
l700 <- which(pa == 700)
l850 <- which(pa == 850)
res <- Humidity_index(
  ta[, , l850, ],
  hus[, , l850, ],
  ta[, , l700, ],
  hus[, , l700, ],
  ta[, , l500, ],
  hus[, , l500, ]
)