code: def _insert_lcl_level(pressure, temperature, lcl_pressure):
      
      interp_temp = interpolate_1d(lcl_pressure, pressure, temperature)
  
      # Pressure needs to be increasing for searchsorted, so flip it and then convert
      # the index back to the original array
      loc = pressure.size - pressure[::-1].searchsorted(lcl_pressure)
      return np.insert(temperature.m, loc, interp_temp.m) * temperature.units
docs: Insert the LCL pressure into the profile.