code: def dsksrf(dsk, bodyid):
      
      dsk    = stypes.stringToCharP(dsk)
      bodyid = ctypes.c_int(bodyid)
      srfids = stypes.SPICEINT_CELL(10000)
      libspice.dsksrf_c(dsk, bodyid, ctypes.byref(srfids))
      return srfids
docs: Find the set of surface ID codes for all surfaces associated with 
      a given body in a specified DSK file. 
      
      https://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/dsksrf_c.html
      
      :param dsk: Name of DSK file.
      :type dsk: str
      :param bodyid: Integer body ID code.
      :type bodyid: int
      :return: Set of ID codes of surfaces in DSK file.