code: def _concatenate_coils(coils):
      """Concatenate MEG coil parameters."""
      rmags = np.concatenate([coil["rmag"] for coil in coils])
      cosmags = np.concatenate([coil["cosmag"] for coil in coils])
      ws = np.concatenate([coil["w"] for coil in coils])
      n_int = np.array([len(coil["rmag"]) for coil in coils])
      if n_int[-1] == 0:
          # We assume each sensor has at least one integration point,
          # which should be a safe assumption. But let's check it here, since
          # our code elsewhere relies on bins[-1] + 1 being the number of sensors
          raise RuntimeError("not supported")
      bins = np.repeat(np.arange(len(n_int)), n_int)
      return rmags, cosmags, ws, bins
  
imports: from copy import deepcopy
  import numpy
  from .._fiff.constants import FIFF
  from ..bem import _import_openmeeg, _make_openmeeg_geometry
  from ..fixes import bincount, jit
  from ..parallel import parallel_func
  from ..surface import _jit_cross, _project_onto_surface
  from ..transforms import apply_trans, invert_transform
  from ..utils import _check_option, _pl, fill_doc, logger, verbose, warn
label_desc: Use image node_id 5 for calling a global function with example usage: _concatenate_coils(coils) and returns: rmags,  cosmags,  ws,  bins