Vertical Wind Shear Magnitude and Direction

Usage

wind_shear(
  ua_s,
  va_s,
  zg_s,
  pa_s,
  ua_pl,
  va_pl,
  zg_pl,
  pa_pl,
  h1,
  h2,
  nthreads = 1
)

ua_s: 10 metre zonal wind component [lon, lat, time] (m/s).
va_s: 10 metre meridional wind component [lon, lat, time] (m/s).
zg_s: surface Geopotential [lon, lat] (m^2/s^2).
pa_s: surface pressure [lon, lat, time] (hPa).
ua_pl: zonal wind component on pressure level [lon, lat, lev, time] (m/s).
va_pl: meridional wind component on pressure level [lon, lat, lev, time] (m/s).
zg_pl: geopotential on pressure levels [lon, lat, lev, time] (m^2/s^2).
pa_pl: pressure profile [lev] in descending manner (hPa).
h1: lower height level for shear calculation (m).
h2: upper height level for shear calculation (m).
nthreads: an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

a list with the computed wind shear magnitude [lon, lat, time] (m/s), and wind shear direction [lon, lat, time] (degree).