def METHOD_NAME(self, grid, field, axis, coord):
    sl = [slice(None), slice(None), slice(None)]
    sl[axis] = slice(coord, coord + 1)
    if grid.ds.field_ordering == 1:
        sl.reverse()
    return self._read_data_set(grid, field)[tuple(sl)]