def METHOD_NAME():
    """
    Create a 6x5 test grid, but with two internal nodes closed.
    This is a sheet flow test.
    """
    mg = RasterModelGrid((6, 5), xy_spacing=(10.0, 10.0))
    mg.set_closed_boundaries_at_grid_edges(True, True, False, True)
    mg.status_at_node[7] = mg.BC_NODE_IS_CLOSED
    mg.status_at_node[16] = mg.BC_NODE_IS_CLOSED
    z = mg.node_x.copy()
    Q_in = np.full(25, 2.0)
    A_target = (
        np.array(
            [
                [0.0, 0.0, 0.0, 0.0, 0.0],
                [1.0, 1.0, 0.0, 1.0, 0.0],
                [6.0, 6.0, 3.0, 1.0, 0.0],
                [0.0, 0.0, 2.0, 1.0, 0.0],
                [3.0, 3.0, 2.0, 1.0, 0.0],
                [0.0, 0.0, 0.0, 0.0, 0.0],
            ]
        ).flatten()
        * 100.0
    )
    frcvr_target = np.array(
        [
            [0, 1, 2, 3, 4],
            [5, 5, 7, 12, 9],
            [10, 10, 11, 12, 14],
            [15, 16, 11, 17, 19],
            [20, 20, 21, 22, 24],
            [25, 26, 27, 28, 29],
        ]
    ).flatten()
    links2rcvr_target = np.full(mg.number_of_nodes, XX)
    links2rcvr_target[mg.core_nodes] = np.array([9, 62, 18, 19, 20, 67, 29, 36, 37, 38])
    steepest_target = np.array(
        [
            [0.0, 0.0, 0.0, 0.0, 0.0],
            [0.0, 1.0, 0.0, 0.0, 0.0],
            [0.0, 1.0, 1.0, 1.0, 0.0],
            [0.0, 0.0, 0.0, 1.0, 0.0],
            [0.0, 1.0, 1.0, 1.0, 0.0],
            [0.0, 0.0, 0.0, 0.0, 0.0],
        ]
    ).flatten()
    steepest_target[np.array([8, 17])] = 1.0 / np.sqrt(2.0)
    mg.add_field("topographic__elevation", z, at="node", units="-")
    class DansGrid:
        pass
    dans_grid = DansGrid()
    dans_grid.mg = mg
    dans_grid.z = z
    dans_grid.Q_in = Q_in
    dans_grid.A_target = A_target
    dans_grid.frcvr_target = frcvr_target
    dans_grid.steepest_target = steepest_target
    dans_grid.links2rcvr_target = links2rcvr_target
    return dans_grid