IMPORTANT NOTE: The example implementation given here is obsolete and should 
not be considered as a model for implementing solvers based on face finite 
elements. For a model implementation that employs general utility subroutines
see instead the solver module 

   .../fem/src/modules/ModelMixedPoisson.F90
