code: def w(self, units=None):
          
  
          if units is None:
              if self.hamiltonian is None:
                  units = dimensionless
              else:
                  units = self.hamiltonian.units
  
          return super(Orbit, self).w(units=units)
docs: This returns a single array containing the phase-space positions.
  
          Parameters
          ----------
          units : `~gala.units.UnitSystem` (optional)
              The unit system to represent the position and velocity in
              before combining into the full array.
  
          Returns
          -------
          w : `~numpy.ndarray`
              A numpy array of all positions and velocities, without units.
              Will have shape ``(2*ndim,...)``.