find_lcl {aiRthermo} | R Documentation |
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.
find_lcl(Ptop, P0, T0, w0, deltaP)
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. |
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. |
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)