code: def _get_pores(self, sampling_points):
          
          pores = []
          for point in sampling_points:
              pores.append(
                  Pore(
                      self.system['elements'],
                      self.system['coordinates'],
                      com=point))
          return pores
docs: Under development.