def METHOD_NAME():
    X = -1
    neighbors = sgrid.neighbor_node_array(
        (3, 3),
        closed_boundary_nodes=[3, 5, 6, 7, 8],
        open_boundary_nodes=[0, 1, 2],
        inactive=X,
    )
    expected = np.array(
        [
            [X, X, X, X],
            [X, 4, X, X],
            [X, X, X, X],
            [X, X, X, X],
            [X, X, X, 1],
            [X, X, X, X],
            [X, X, X, X],
            [X, X, X, X],
            [X, X, X, X],
        ]
    )
    assert_array_equal(neighbors.T, expected)