Deep Convective Index

Usage

DCI(ta_s, hus_s, pa_s, ta500, ta850, hus850, nthreads = 1)

Arguments

ta_s: parcel temperature [lon, lat, time] at 2 meter above the surface (K).
hus_s: parcel specific humidity [lon, lat, time] at 2 meter above the surface (kg/kg).
pa_s: surface pressure [lon, lat, time] (hPa).
ta500: air temperature [lon, lat, time] at 500 hPa (K).
ta850: air temperature [lon, lat, time] at 850 hPa (K).
hus850: specific humidity [lon, lat, time] at 850 hPa (kg/kg).
nthreads: an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

Value

the computed deep convective index [lon, lat, time] (K).

References

Barlow, W.R., 1993. A new index for the prediction of deep convection. Preprints, 17th Conf. on Severe Local Storms. Amer. Meteor., St. Louis, MO, pp. 129–132.

Examples

data("ERA5_sfc")
ta_s <- ERA5_sfc$t2m
td_s <- ERA5_sfc$t2d
pa_s <- ERA5_sfc$sp / 100
hus_s <- hus_from_td(td_s, pa_s)
dim(hus_s) <- dim(ta_s)

data("ERA5_pl")
ta <- ERA5_pl$ta
hus <- ERA5_pl$hus
#' convert to hPa
pa <- ERA5_pl$plev / 100
l500 <- which(pa == 500)
l850 <- which(pa == 850)
res <- DCI(
  ta_s,
  hus_s,
  pa_s,
  ta[, , l500, ],
  ta[, , l850, ],
  hus[, , l850, ]
)