calc_lamb {meteoEVT} | R Documentation |
Lamb vector (sometimes called vortex energy)
Description
Calculates the Lamb vector (cross product of wind vector and vorticity vector)
Usage
calc_lamb(
u_fld,
v_fld,
w_fld,
lev_p,
lat = NULL,
dx = 0.25,
dy = 0.25,
relative = TRUE,
zvort_fld = NULL,
mode = "lonlat"
)
Arguments
u_fld |
zonal velocity field [m/s] |
v_fld |
meridional velocity field [m/s] |
w_fld |
vertical velocity field [m/s] |
lev_p |
vector containing pressure levels [Pa] |
lat |
vector containing latitude |
dx |
x resolution in the corresponding unit (e.g. 0.25 degree for ERA5 with |
dy |
y resolution in the corresponding unit (e.g. 0.25 degree for ERA5 with |
relative |
logical, TRUE: only relative vorticity, FALSE: whole (absolute) vorticity should be used for calculation of enstrophy, default: TRUE |
zvort_fld |
optional zvort field (if e.g., zvort is directly taken from ERA5 and not calculated separately) |
mode |
use 'lonlat' if the data is given on a lon-lat-grid or 'cartesian' if the data is given on an equidistant cartesian grid |
Value
lamb vector [m/s^2]
Examples
myfile=system.file("extdata", "era5_storm-zeynep.nc", package = "meteoEVT")
data = readin_era5(myfile)
lamb=calc_lamb(data$u,data$v,data$w,data$lev,lat=data$lat)