aiRthermo (version 1.1)

adiabatic_ascent: Properties of an air parcel after adiabatic ascent

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

# NOT RUN {
P0<-101325
T0<-273.15
w0<-0.0025
adiabEvov<-adiabatic_ascent(P0,T0,w0,50000,5)
# }

Run the code above in your browser using DataCamp Workspace