Bradbury Index

Bradbury Index

Usage

Bradbury_Index(hus850, ta850, hus500, ta500, nthreads = 1)

Arguments

hus850

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

ta850

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

hus500

specific humidity [lon, lat, time] at 850 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 Bradbury Index [lon, lat, time] (K).

References

Bradbury, T.A.M., 1977. The use of wet-bulb potential temperature charts. Meteorol. Mag. 106, 233–251.

Examples

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

l850 <- which(plev == 850)
l500 <- which(plev == 500)

hus850 <- ERA5_pl$hus[, , l850, ]
ta850 <- ERA5_pl$ta[, , l850, ]

hus500 <- ERA5_pl$hus[, , l500, ]
ta500 <- ERA5_pl$ta[, , l500, ]

res <- Bradbury_Index(hus850, ta850, hus500, ta500, nthreads = 1)