implement new method to compute a posteriori error estimation (F.Casenave). Need to debug the execution step.