Brunt-Vaisalla (angular) frequency (\(s^{-2}\)) considering hydrostatic equilibrium. P is used as a vertical level.
bruntVaisallaOmegaSquared(Ps, Ts, ws, consts = export_constants())
A vector with pressure values (Pa)
A vector with temperature values (K)
A vector with mixing ratio values (kg/kg)
The constants defined in aiRthermoConstants data are necessary. The constants g and Rd are used.
omega_sq
The angular frequency (squared, (\(s^{-2}\))) is returned in order to avoid complex numbers.
PT2Theta and densityMoistAir are used inside bruntVaisallaOmegaSquared.
PT2Theta
densityMoistAir
bruntVaisallaOmegaSquared
# NOT RUN { data(RadiosondeD) dPs<-RadiosondeD[,1]*100 dTs<-C2K(RadiosondeD[,3]) dws<-RadiosondeD[,6]/1000 bruntVaisallaOmegaSquared(dPs,dTs,dws) # }
Run the code above in your browser using DataCamp Workspace