def METHOD_NAME(vslice_geo):
    da = vslice_geo.read().Temperature
    da2 = da.sel(x=10.8, y=55.6, z=-3)
    assert type(da2.geometry) == GeometryPoint3D
    assert da2.geometry.x == pytest.approx(10.802878)
    assert da2.geometry.y == pytest.approx(55.603096)
    assert da2.geometry.z == pytest.approx(-2.3)
    assert da2.n_timesteps == da.n_timesteps