Skip to contents

Yonetani Index

Yonetani.Rd

Yonetani Index

Usage

Yonetani(ta_pl, hus_pl, zg_pl, pa_pl, modified = FALSE, nthreads = 1)

Arguments

ta_pl

air temperature [lon, lat, lev, time] on pressure level (K).

hus_pl

specific humidity [lon, lat, lev, time]on pressure levels (kg/kg).

zg_pl

geopotential [lon, lat, lev, time] on pressure levels (m^2/s^2)

pa_pl

pressure levels [lev] (hPa).

modified

logical value inidcating if to calculate the modified Yonetani Index (TRUE) or the typical Yonetani Index (FALSE). Default is FALSE.

nthreads

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

Value

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

References

  • Yonetani, T., 1979. Instability index which takes account of the stabilities of a lower atmospheric layer for the purpose of thunderstorm prediction in the northern Kanto Plain. Rep. Natl. Res. Center Disaster Prev., 35–44 (in Japanese).

  • Jacovides, C. P., & Yonetani, T. (1990). An evaluation of stability indices for thunderstorm prediction in Greater Cyprus. Weather and forecasting, 5(4), 559-569.

Examples

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

res <- Yonetani(ta_pl, hus_pl,
  zg_pl, pa_pl,
  modified = FALSE
)