docs: def distance(self,value):
          """New distance value must be a ``Quantity`` object
          """
          self.stars['distance'] = value.to('pc').value
  
          old_distmod = self.stars['distmod'].copy()
          new_distmod = distancemodulus(self.stars['distance'])
  
          for m in self.bands:
              self.stars[m] += new_distmod - old_distmod
  
          self.stars['distmod'] = new_distmod
  
          logging.warning('Setting the distance manually may have screwed up your constraints.  Re-apply constraints as necessary.')