Adedokun Index I & II

Usage

Adedokun_Index(hus, pa, ta, ta500, nthreads = 1)

Arguments

hus: specific humidity [lon, lat, time], either at 850 hPa or at the surface for Adedokun_Index I, II, respectively (kg/kg).
pa: air pressure [lon, lat, time] (hPa).
ta: air temperature [lon, lat, time], either at 850 hPa or at the surface for Adedokun_Index I, II, respectively (K).
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 Adedokun Index I or II [lon, lat, time] (K).

References

Adedokun, J.A., 1981. Potential instability and precipitation occurrence within an inter-tropical discontinuity environment. Arch. Meteor. Geophys. Bioklimatol., A 30, 69–86.

Adedokun, J.A., 1982. On an instability index relevant to precipitation forecasting in West Africa. Arch. Meteor. Geophys. Bioklimatol., A 31, 221–230.

Examples

data("ERA5_pl")
plev <- ERA5_pl$plev / 100 # convert to hPa
# Level 850 hpa
l850 <- which(plev == 850)
hus <- ERA5_pl$hus[, , l850, ]
ta <- ERA5_pl$ta[, , l850, ]
# create 3D array
pa <- array(rep(850, prod(dim(hus))),
  dim = dim(hus)
)

l500 <- which(plev == 500)
ta500 <- ERA5_pl$ta[, , l500, ]
# Adedokun_Index I
res <- Adedokun_Index(hus, pa, ta, ta500, nthreads = 1)