def METHOD_NAME():
    mg = RasterModelGrid((10, 10))
    z = mg.add_zeros("topographic__elevation", at="node")
    z += mg.x_of_node + mg.y_of_node
    fa = FlowAccumulator(mg)
    fa.run_one_step()
    with pytest.raises(ValueError):
        ErosionDeposition(mg, phi=0)