calc_tfp {meteoEVT}

Thermic Front Parameter (TFP)

Description

Calculates the thermic front parameter based on the potential temperature

Usage

calc_tfp(t_fld, lev_p, lat = NULL, dx = 0.25, dy = 0.25, mode = "lonlat")

Arguments

t_fld

temperature field [K]

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

thermic front parameter [K/m^2]

Examples

myfile=system.file("extdata", "era5_storm-zeynep.nc", package = "meteoEVT")
data = readin_era5(myfile)
tfp=calc_tfp(data$temp,data$lev,data$lat)