def test_grid_layer_slice():
"""Test grid slice coordinates."""
grd = Grid()
grd.from_file(REEKFILE)
sarr1, _ibarr = grd.get_layer_slice(1)
sarrn, _ibarr = grd.get_layer_slice(grd.nlay, top=False)
cell1 = grd.get_xyz_cell_corners(ijk=(1, 1, 1))
celln = grd.get_xyz_cell_corners(ijk=(1, 1, grd.nlay))
celll = grd.get_xyz_cell_corners(ijk=(grd.ncol, grd.nrow, grd.nlay))
assert sarr1[0, 0, 0] == cell1[0]
assert sarr1[0, 0, 1] == cell1[1]
assert sarrn[0, 0, 0] == celln[12]
assert sarrn[0, 0, 1] == celln[13]
assert sarrn[-1, 0, 0] == celll[12]
assert sarrn[-1, 0, 1] == celll[13]
