find_lcl {aiRthermo}

Calculation of the Lifted Condensation Level (LCL)

Description

For a particle with initial conditions P_0, T_0 and w_0, this function performs an adiabatic vertical evolution until it gets saturated at most when Ptop is reached.

Usage

find_lcl(Ptop, P0, T0, w0, deltaP)

Arguments

Ptop

Maximun level pressure selected (Pa).

P0

Initial value of pressure (Pa).

T0

Initial value of temperature (K).

w0

Initial value of mixing ratio (kg/kg).

deltaP

The width (Pa) of the layers used in the calculation of the numerical solution for the vertical evolution. A default value of 5 Pa is used.

Value

Returns a list including the following values:

Plcl

The pressure at LCL (Pa).

Tlcl

The temperature at LCL (K).

wlcl

The mixing ratio at LCL (kg/kg).

thetalcl

The potential temperature at LCL (K).

gotit

0 or 1 whether the particle arrived or not to saturation (LCL) before arriving to Ptop.

Examples

Ptop=50000
P0=101325
T0=273.15
w0=0.0025
deltaP=5
rh=100*w0/saturation_mixing_ratio(P0,T0,export_constants())
fndlcl=find_lcl(Ptop,P0,T0,w0,deltaP)