Since the removal of this check now allows ```izip(tyx_cube, zyx_cube)```, which will use ```np.ndindex``` to give every (t,z) combination, the documentation should also be updated to explain what will happen: _(t0, z0), (t0, z1), ... (t0, zn), (t1, z0), (t1, z1), ..._