find_lcl: 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
# NOT RUN {Ptop=50000P0=101325T0=273.15w0=0.0025deltaP=5rh=100*w0/saturation_mixing_ratio(P0,T0,export_constants())
fndlcl=find_lcl(Ptop,P0,T0,w0,deltaP)
# }