Usage
Showalter_Index(ta850, hus850, ta500, nthreads = 1)
Arguments
- ta850
air temperature [lon, lat, time] at 850 hPa (K).
- hus850
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 Showalter Index [lon, lat, time] (K).
References
Showalter, A.K., 1953. A stability index for thunderstorm forecasting.
Bull. Am. Meteorol. Soc. 34, 250–252.
Examples
data("ERA5_pl")
pa <- ERA5_pl$plev / 100 # convert to hPa
hus <- ERA5_pl$hus
ta <- ERA5_pl$ta
l500 <- which(pa == 500)
l850 <- which(pa == 850)
res <- Showalter_Index(
ta[, , l850, ],
hus[, , l850, ],
ta[, , l500, ]
)
