Boyden Index

Boyden Index

Usage

Boyden_index(z700, z1000, ta700, nthreads = 1)

Arguments

z700

geopotential height [lon, lat, time] at 700 hPa (dam).

z1000

geopotential height [lon, lat, time] at 1000 hPa (dam).

ta700

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

nthreads

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

Value

the computed Boyden Index [lon, lat, time].

References

Boyden, C.J., 1963. A simple instability index for use as a synoptic parameter. Meteorol. Mag. 92, 198–210.

Examples

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

l700 <- which(plev == 700)
z700 <- ERA5_pl$z[, , l700, ]
z700 <- z700 / (9.80665 * 100) # convert to decametres

l1000 <- which(plev == 1000)
z1000 <- ERA5_pl$z[, , l1000, ]
z1000 <- z1000 / (9.80665 * 100) # convert to decametres

ta700 <- ERA5_pl$ta[, , l700, ]
res <- Boyden_index(z700, z1000, ta700, nthreads = 1)