adiabatic_ascent: Properties of an air parcel after adiabatic ascent

In aiRthermo: Atmospheric Thermodynamics and Visualization

Description

A particle located at Pstart pressure (Pa), Tstart temperature (K) and wstart mixing ratio (kg/kg) ascends (pseudo)adiabatically to Pend (Pa). The evolution is computed by numerically integrating the dT/dP ordinary differential equation (ODE) using a 4th order Runge-Kutta scheme, assuming hydrostatic equilibrium and that the particle is saturated after the Lifted Condensation Level (LCL).

Usage

adiabatic_ascent(Pstart, Tstart, wstart, Pend, deltaP = 1)

Arguments

Pstart	Initial value for pressure (Pa).
Tstart	Initial value for temperature (K).
wstart	Initial value for mixing ratio (kg/kg).
Pend	End value for pressure (Pa).
deltaP	deltaP (Pa) represents the numerical increment used for integrating the Ordinary Differential Equation (ODE) representing the vertical evolution.

Value

The function returns a list that includes Tend (final value of temperature) and mixRatioEnd (mixing ratio of the air parcel at the end of the evolution).

Tend	Temperature at the end (K).
mixRatioEnd	Mixing ratio at the end (kg/kg).

Examples

P0<-101325
T0<-273.15
w0<-0.0025
adiabEvov<-adiabatic_ascent(P0,T0,w0,50000,5)

Example output

aiRthermo documentation built on May 1, 2019, 9:24 p.m.