  We introduce interpolation operators with approximation and stability
properties suited for parabolic problems in primal and mixed formulations. We
derive localized error estimates for tensor product meshes (occurring in
classical time-marching schemes) as well as locally in space-time refined
meshes.
