docs: def lifting_condensation_level(T, RH):
      '''Compute the Lifiting Condensation Level (LCL) for a given temperature and relative humidity
  
      Inputs:  T is temperature in Kelvin
              RH is relative humidity (dimensionless)
  
      Output: LCL in meters
  
      This is height (relative to parcel height) at which the parcel would become saturated during adiabatic ascent.
  
      Based on approximate formula from Bolton (1980 MWR) as given by Romps (2017 JAS)
  
      For an exact formula see Romps (2017 JAS), doi:10.1175/JAS-D-17-0102.1
      '''
      Tadj = T-55.  # in Kelvin
      return cp/g*(Tadj - (1/Tadj - log(RH)/2840.)**(-1))