Convective Instability (Barber)

Usage

CI_Barber(ta_s, hus_s, pa_s, ta_pl, hus_pl, pa_pl, nthreads = 1)

Arguments

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

Value

the computed convective stability index [lon. lat, time] in K defined by Barber due to Peppler (1988).

References

Peppler, R. A., & Lamb, P. J. (1989). Tropospheric Static Stability and Central North American Growing Season Rainfall. Monthly Weather Review, 117(6), 1156–1180. doi:10.1175/1520-0493(1989)117<1156:TSSACN>2.0.CO;2
Barber, D. A., 1975: A contribution to the climatology of static stability and vertical wind shear. Preprints. Ninth Conf. Severe Local Storms, Norman, Amer. Meteor. S o c , 13-17.

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(td_s)

data("ERA5_pl")
ta_pl <- ERA5_pl$ta
hus_pl <- ERA5_pl$hus
pa_pl <- ERA5_pl$plev / 100

res <- CI_Barber(
  ta_s, hus_s, pa_s,
  ta_pl, hus_pl, pa_pl
)