calc_fdiag {meteoEVT}

F diagnostic

Description

Calculates the F diagnostic

Usage

calc_fdiag(
  t_fld,
  u_fld,
  v_fld,
  w_fld,
  lev_p,
  lat = NULL,
  dx = 0.25,
  dy = 0.25,
  mode = "lonlat"
)

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]

lev_p

vector containing pressure levels [Pa]

lat

only for lonlat mode: 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')

mode

the horizontal coordinate system, options are lonlat for a longitude-latitude-grid (default), or cartesian for an equidistant cartesian grid

Value

F diagnostic (dimensionless)

Examples

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