def METHOD_NAME(loadsfile1):
    cube1 = loadsfile1
    surf1 = xtgeo.surface_from_cube(cube1, 1100.0)
    surf2 = xtgeo.surface_from_cube(cube1, 2900.0)
    # values taken from IB spreadsheet except variance
    attributes = {
        "min": -85.2640,
        "max": 250.0,
        "rms": 196.571,
        "mean": 157.9982,
        "var": 13676.7175,
        "sumneg": -2160.80,
        "sumpos": 73418,
        "meanabs": 167.5805,
        "meanpos": 194.7427,
        "meanneg": -29.2,
        "maxneg": -85.264,
        "maxpos": 250.0,
    }
    attrs = surf1.slice_cube_window(
        cube1,
        other=surf2,
        other_position="below",
        attribute=list(attributes.keys()),
        sampling="cube",
        snapxy=True,
        ndiv=None,
        algorithm=2,
        showprogress=False,
    )
    for name in attributes.keys():
        print(attrs[name].values.mean())
        assert attributes[name] == pytest.approx(attrs[name].values.mean(), abs=0.001)
    attrs = surf1.slice_cube_window(
        cube1,
        other=surf2,
        other_position="below",
        attribute=list(attributes.keys()),
        sampling="trilinear",
        snapxy=True,
        ndiv=None,
        algorithm=2,
    )
    # less strict abs using relative 1%:
    for name in attributes.keys():
        assert attributes[name] == pytest.approx(attrs[name].values.mean(), rel=0.01)