The dimension is non-standard in the sense that we only have one dof per face but require the problem for convenience to return a full dimension vector for each coordinate direction.
TODO: add similar check for
-
dirichlet and other stuff returning primary varaibles/numeqvector