issue_comments

@shoyer - my experience is that the dummy netCDF4.datetime objects don't play nice with setting up a pandas time index, so a intermediate conversion step is necessary. I haven't looked into why this is exactly.

I just tried the new decoding and it seems to work.

OK, I just merged a fix into master.

Unfortunately, it's not terribly useful to be able to have arrays decoded as netCDF4.datetime objects, because it's not possible to use them with label based indexing, as they are not hashable (see Unidata/netcdf4-python#255).

Just out of curiosity, why do you usually convert netCDF4.datetime objects into real datetime objects? I'm guessing it's not because you want objects you can put in a dictionary?

If there is a better type than netCDF4.datetime to use for decoded dates, I'm definitely willing to consider it. As I said before, we don't make much use of non-standard dates.

It's also certainly possible (in principle) to keep around another array with the original, encoded dates. Right now all the decoding according to CF conventions is done in one large function with no options, but I would love for it to be more flexible and modular.

Thanks, the decode_cf keyword should get me around the problem for now.

I've made a habit of always directly converting my netCDF4.datetime to true datetime.datetime objects immediately, since netCDF4 only returns real datetime objects for the Gregorian calendars.

python f = netCDF4.Dataset('sample_for_xray.nc') decoded_times = netCDF4.num2date(f.variables['time'][:], f.variables['time'].units, f.variables['time'].calendar) for i, t in enumerate(decoded_times): decoded_times[i] = datetime.datetime(*t.timetuple()[:6])

The important piece to remember if this is done is that you have to be very picky about how you calculate timedeltas between these dates since they think they are on the Gregorian calendar. I usually just keep an ordinal based time array around for that reason.

Ouch! Thanks for filing the report and providing the sample file -- I will take a look.

For now, turn off automatic date decoding by calling xray.open_dataset('sample_for_xray.nc', decode_cf=False).

I'm guessing that part of the trouble might be that numpy and pandas provide poor support for alternative calendars (and honestly, I haven't tested them very much). I attempted to fall back on making arrays of python datetime objects, but in this case it looks like that didn't work -- somehow things got converted in a numpy native datetime64 array anyways.