def METHOD_NAME(hydds, flwdir):
    flwda1, _ = flw.upscale_flwdir(hydds, flwdir, 3, uparea_name="uparea", method="dmm")
    assert hydds["flwdir"].shape[0] / 3 == flwda1.shape[0]
    assert np.all([v in flwda1.coords for v in ["x_out", "y_out", "idx_out"]])
    # errors
    with pytest.raises(ValueError, match="Flwdir and ds dimensions do not match"):
        flw.upscale_flwdir(hydds.isel(x=slice(1, -1)), flwdir, 3, method="dmm")