docs: def get_pore_surface_parameters(surface_area):
      """ Get input parameters for pore surface binary.
  
      Get input parameters for pore_surface binary from zeo++ output,
      while keeping data provenance.
      """
      PoreSurfaceParameters = DataFactory('phtools.surface')
      d = {
          'accessible_surface_area': surface_area.get_dict()['ASA_A^2'],
          'target_volume': 40e3,
          'sampling_method': 'random',
      }
      return PoreSurfaceParameters(dict=d)