code: def polar_direction_xyz(self):
          
          return mesh.spin_in_roche(self.polar_direction_uvw,
                                    self.true_anom, self.elongan, self.eincl).astype(float)
docs: get current polar direction in Roche (xyz) coordinates