def METHOD_NAME() -> CubeList:
    """Set up a single realization cube in parameter space"""
    realizations = [0, 1, 2, 3]
    data = np.ones((len(realizations), 2, 2), dtype=np.float32)
    times = [datetime(2017, 11, 10, hour) for hour in [4, 5, 6]]
    cubes = CubeList()
    for time in times:
        cubes.append(
            set_up_variable_cube(
                data,
                realizations=realizations,
                spatial_grid="equalarea",
                time=time,
                frt=datetime(2017, 11, 10, 1),
            )
        )
    cube = cubes.merge_cube()
    sliced_cubes = CubeList(cube.slices_over("realization"))
    [
        s.attributes.update({"history": f"20171110T{i:02d}00Z"})
        for i, s in enumerate(sliced_cubes)
    ]
    return sliced_cubes