Skip to contents

Vertically Integrated Moisture Flux Divergence (VIMFD)

Usage

VIMFD(hus, ua, va, lon, lat, pa_pl, nthreads = 1)

Arguments

hus

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

ua

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

va

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

lon

longitude coordinates [lon] (degrees east).

lat

latitude coordinates [lat] (degrees north).

pa_pl

pressure levels [lev] (hPa).

nthreads

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

Value

array [lon, lat, time] the computed VIMFC (kg m^-2 s^-1).

Examples

data("ERA5_pl")
hus <- ERA5_pl$hus
ua <- ERA5_pl$ua
va <- ERA5_pl$va
lon <- ERA5_pl$lon
lat <- ERA5_pl$lat
# to hPa
pa_pl <- ERA5_pl$plev / 100
# VIMFD
res <- VIMFD(hus, ua, va, lon, lat, pa_pl, nthreads = 1)