Vertical interpolation (Vertical Velocity) — Vintp

Interpolates 3D vertical velocity from model levels to pressure levels.

Usage

Vintp_omega(omega_ml, ps, ap_full, b_full, ap_half, b_half, pl, nthreads = 1)

omega_ml: vertical velocity [lon, lat, lev, time] on model levels (Pa/s)
ps: surface pressure [lon, lat, time] (hPa).
ap_full: vertical coordinate formula term ap(k) or a(k) [lev] at layer midpoints (hPa).
b_full: vertical coordinate formula term b(k) [lev] at layer midpoints. b(k) is dimensionless fraction.
ap_half: vertical coordinate formula term ap(k) or a(k) [lev+1] at layer interfaces (hPa).
b_half: vertical coordinate formula term b(k) [lev+1] at layer interfaces. b(k) is dimensionless fraction.
pl: vector of the required pressure levels [plev] e.g., c(1000,950,850) in hPa.
nthreads: an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

vertical velocity [lon, lat, plev, time] on pressure levels (Pa/s).