code: def rpole_to_pot_aligned(rpole, sma, q, F, d, component=1):
      
      q = q_for_component(q, component=component)
      rpole_ = np.array([0, 0, rpole/sma])
      logger.debug("libphoebe.roche_Omega(q={}, F={}, d={}, rpole={})".format(q, F, d, rpole_))
      pot = libphoebe.roche_Omega(q, F, d, rpole_)
      return pot_for_component(pot, component, reverse=True)
docs: Transforms polar radius to surface potential