def METHOD_NAME(hydds, flwdir):
    # simple test because deprecated
    with pytest.warns(DeprecationWarning, match="basin_shape is deprecated"):
        gdf = flw.basin_shape(hydds, flwdir)
        assert isinstance(gdf, gpd.GeoDataFrame)
        with pytest.raises(ValueError, match="flwdir and ds dimensions do not match"):
            flw.basin_shape(hydds.isel(x=slice(1, -1)), flwdir)