  In the context of model order reduction of parametric elliptic problems, we
present a methodology to reconstruct a conforming flux from a given reduced
solution, that is locally conservative with respect to the underlying finite
element grid. All components of the procedure depend separably on the parameter
and allow for further use in offline/online decomposed computations, for
instance in the context of a posterior error estimation or flow problems.
