import iris import numpy as np iris.__version__ points = np.arange(11) bounds = np.array(([(x-2, x+2) for x in points])) example_coord = iris.coords.AuxCoord(points, bounds=bounds, long_name='example') print example_coord value = 5 for cell in example_coord.cells(): print('{!r} == {!r} is {}'.format(cell, value, cell==value)) for cell in example_coord.cells(): print('{!r} == {!r} is {}'.format(cell.point, value, cell.point==value)) example_coord = iris.coords.DimCoord(points, bounds=bounds, long_name='example') dummy_cube = iris.cube.Cube(np.random.random(11)) dummy_cube.add_dim_coord(example_coord, 0) cell_comparison_constraint = iris.Constraint(example=lambda cell: cell == 5) cell_point_comparison_constraint = iris.Constraint(example=lambda cell: cell.point == 5) print dummy_cube.extract(cell_comparison_constraint).coords() print dummy_cube.extract(cell_point_comparison_constraint).coords()