code: def pot_to_rpole_aligned(pot, sma, q, F, d, component=1):
      
      q = q_for_component(q, component=component)
      Phi = pot_for_component(pot, q, component=component)
      logger.debug("libphobe.roche_pole(q={}, F={}, d={}, Omega={})".format(q, F, d, pot))
      return libphoebe.roche_pole(q, F, d, pot) * sma
docs: Transforms surface potential to polar radius