powered by
Calculates the potential temperature
calc_theta(t_fld, lev_p)
density [kg/m^3]
temperature field [K]
vector containing pressure levels [Pa]
myfile=system.file("extdata", "era5_storm-zeynep.nc", package = "meteoEVT") data = readin_era5(myfile) theta=calc_theta(data$temp,data$lev)
Run the code above in your browser using DataLab