Vertical Velocity (hl)

`r lifecycle::badge('experimental')` This function estimate the vertical velocity from model output. It works only with models on hybrid height coordinate.

Usage

omega_hl(ua, va, ps, ap_full, b_full, ap_half, b_half, lon, lat, nthreads = 1)

Arguments

ua

zonal wind component [lon, lat, lev, time] on model levels (m/s).

va

meridional wind component [lon, lat, lev, time] on model levels (m/s).

ps

surface pressure [lon, lat, time] (Pa).

ap_full

vertical coordinate formula term a(k) [lev] at layer midpoints (Pa).

b_full

vertical coordinate formula term b(k) [lev] at layer midpoints. b(k) is dimensionless fraction.

ap_half

vertical coordinate formula term a(k) [lev+1] at layer interfaces (Pa).

b_half

vertical coordinate formula term b(k) [lev+1] at layer interfaces. b(k) is dimensionless fraction.

lon

longitude coordinates [lon] (degrees east).

lat

latitude coordinates [lat] (degrees north).

nthreads

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

Value

the vertical velocity [lon, lat, lev, time] in Pa/s. Lev is from the top of the atmosphere to the surface.

Note

The function handles model output organised top of the atmosphere to the surface and vice versa.