def METHOD_NAME():
    grid = RasterModelGrid((3, 4))
    assert np.all(
        grid.links_at_patch
        == np.array(
            [
                [4, 7, 3, 0],
                [5, 8, 4, 1],
                [6, 9, 5, 2],
                [11, 14, 10, 7],
                [12, 15, 11, 8],
                [13, 16, 12, 9],
            ],
            dtype=int,
        )
    )
    with raises(ValueError):
        grid.links_at_patch[0] = -1