code: def north_pole_uvw(self):
          
          # TODO: is this rpole scaling true for all distortion_methods??
          rpole = self.instantaneous_rpole*self.sma
          return self.polar_direction_uvw*rpole+self.mesh._pos
docs: location of the north pole in the global/system frame