moist_lapse_rate

typhon.physics.moist_lapse_rate(p, T, e_eq=None)[source]

Calculate the moist-adiabatic temperature lapse rate.

Bohren and Albrecht (Equation 6.111, note the sign change):

\[\frac{dT}{dz} = \frac{g}{c_p} \frac{1 + l_v w_s / RT}{1 + l_v^2 w_s/c_p R_v T^2}\]

Parameters:

p (float or ndarray) – Pressure [Pa].
T (float or ndarray) – Temperature [K].
e_eq (callable) – Function to calculate the equilibrium vapor pressure of water in Pa. The function must implement the signature e_eq = f(T) where T is temperature in Kelvin. If None the function e_eq_water_mk() is used.

Returns:

Moist-adiabatic lapse rate [K/m].

Return type:

float or ndarray

Examples

>>> moist_lapse_rate(1013.25e2, 288.15)
0.004728194612232855

References

Bohren C. and Albrecht B., Atmospheric Thermodynamics, p. 287-92