code: def configure_box(self, boxsize, root_nx=1, root_ny=1, root_nz=1):
          
          clibrebound.reb_configure_box(byref(self), c_double(boxsize), c_int(root_nx), c_int(root_ny), c_int(root_nz))
          return
docs: Initialize the simulation box.
  
          This function only needs to be called it boundary conditions other than "none"
          are used. In such a case the boxsize must be known and is set with this function.
  
          Parameters
          ----------
          boxsize : float, optional
              The size of one root box.
          root_nx, root_ny, root_nz : int, optional
              The number of root boxes in each direction. The total size of the simulation box
              will be ``root_nx * boxsize``, ``root_ny * boxsize`` and ``root_nz * boxsize``.
              By default there will be exactly one root box in each direction.