meteoEVT (version 0.1.0)

calc_dsi: Dynamic State Index (DSI)

Description

Calculates the dynamic state index DSI

Usage

calc_dsi(
  t_fld,
  u_fld,
  v_fld,
  w_fld,
  phi_fld,
  lev_p,
  lat = NULL,
  dx = 0.25,
  dy = 0.25,
  zvort_only = FALSE,
  relative = FALSE,
  pv_fld = NULL,
  mode = "lonlat"
)

Value

dynamic state index [K^2*m^4/(kg^2*s^3)]

Arguments

t_fld

temperature field [K]

u_fld

zonal velocity field [m/s]

v_fld

meridional velocity field [m/s]

w_fld

vertical velocity field [m/s]

phi_fld

geopotential height [gpm]

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 mode='lonlat' or e.g. 1000 m in cartesian coordinates with mode='cartesian')

dy

y resolution in the corresponding unit (e.g. 0.25 degree for ERA5 with mode='lonlat' or e.g. 1000 m in cartesian coordinates with mode='cartesian')

zvort_only

logical, TRUE: if only the vertical vorticity (zvort) should be calculated, FALSE: for the whole vorticity vector, default: FALSE

relative

logical, TRUE: only relative vorticity, FALSE: whole (absolute) vorticity, default: FALSE

pv_fld

optional pv field (if e.g., PV 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

Examples

myfile=system.file("extdata", "era5_storm-zeynep.nc", package = "meteoEVT")
data = readin_era5(myfile)
dsi=calc_dsi(data$temp,data$u,data$v,data$w,data$z,lev_p=data$lev,lat=data$lat)