PartialDate extraction bug

[bblay]

I think the last extract should return something too, but does not.

>>> import iris
>>> from iris.time import PartialDateTime
>>> iris.FUTURE.cell_datetime_objects = True
>>> cube = iris.load_cube("/data/local/itbb/git/iris-test-data/test_data/PP/globClim1/theta.pp")
>>> print cube
air_potential_temperature / (K)     (model_level_number: 38; latitude: 145; longitude: 192)
     Dimension coordinates:
          model_level_number                           x             -               -
          latitude                                     -             x               -
          longitude                                    -             -               x
     Auxiliary coordinates:
          level_height                                 x             -               -
          sigma                                        x             -               -
     Scalar coordinates:
          forecast_period: 26280 hours
          forecast_reference_time: 2005-01-01 00:00:00
          time: 2008-01-01 00:00:00, bound=(2007-12-01 00:00:00, 2008-01-01 00:00:00)
     Attributes:
          STASH: m01s00i004
          source: Data from Met Office Unified Model 6.06
     Cell methods:
          mean: time (1 hour)
>>> cube.extract(iris.Constraint(time=lambda cell: cell == PartialDateTime(year=2008)))
<iris 'Cube' of air_potential_temperature / (K) (model_level_number: 38; latitude: 145; longitude: 192)>
>>> cube.extract(iris.Constraint(time=lambda cell: cell == PartialDateTime(day=1)))
<iris 'Cube' of air_potential_temperature / (K) (model_level_number: 38; latitude: 145; longitude: 192)>
>>> cube.extract(iris.Constraint(time=lambda cell: cell == PartialDateTime(month=1)))
>>> 

[bblay]

I'm assuming this is because the lower bound month is 12, but is it correct to return nothing?

[esc24]

This is a real problem and one I've looked at in the dim and distant past. I just took a quick look and have the beginnings of a fix.

[pelson]

It essentially boils down to:

from iris.time import PartialDateTime
import datetime
cell = Cell(datetime.datetime(2012, 1, 1),
         [datetime.datetime(2011, 12, 1), datetime.datetime(2012, 3, 1)])
print cell.contains_point(PartialDateTime(month=1))

[ajdawson]

This looks to have been closed by #1016. Please reopen if that is incorrect.

[esc24]

That's right. Iris will now raise an exception rather than giving the wrong answer. It's frustrating that we can't do this though, at least in the case of a standard calendar.