code: def box_model_domain(num_points=2, **kwargs):
      
      ax = Axis(axis_type='abstract', num_points=num_points)
      boxes = _Domain(axes=ax, **kwargs)
      boxes.domain_type = 'box'
      return boxes
docs: Creates a box model domain (a single abstract axis).
  
      :param int num_points:      number of boxes [default: 2]
      :returns:                   Domain with single axis of type ``'abstract'``
                                  and ``self.domain_type = 'box'``
      :rtype:                     :class:`_Domain`
  
      :Example:
  
          ::
  
              >>> from climlab import domain
              >>> box = domain.box_model_domain(num_points=2)
  
              >>> print box
              climlab Domain object with domain_type=box and shape=(2,)