def METHOD_NAME():
    mg = RasterModelGrid((10, 10))
    mg.add_field("topographic__elevation", mg.node_y + mg.node_x, at="node")
    fa = FlowAccumulator(mg, flow_director="FlowDirectorMFD")
    fa.run_one_step()
    # this one should flow half to the west and half to the south
    w_links = mg.adjacent_nodes_at_node[:, 2]
    s_links = mg.adjacent_nodes_at_node[:, 3]
    node_ids = np.arange(mg.number_of_nodes)
    true_receivers = -1 * np.ones(fa.flow_director._receivers.shape)
    true_receivers[mg.core_nodes, 2] = w_links[mg.core_nodes]
    true_receivers[mg.core_nodes, 3] = s_links[mg.core_nodes]
    true_receivers[mg.boundary_nodes, 0] = node_ids[mg.boundary_nodes]
    true_proportions = np.zeros(fa.flow_director._proportions.shape)
    true_proportions[mg.boundary_nodes, 0] = 1
    true_proportions[mg.core_nodes, 2:] = 0.5
    assert_array_equal(true_receivers, fa.flow_director._receivers)
    assert_array_equal(true_proportions, fa.flow_director._proportions)