Hi everyone, Apologies for sending out so many e-mails in such a short amount of time. I guess I have a two part question: 1) Is it possible to concatenation two or more function spaces? In FEniCS I was able to do something like this: V = VectorFunctionSpace(mesh, "CG", 2) Q = FunctionSpace(mesh, "CG", 1) P = FunctionSpace(mesh, "DG", 0) W = V * (Q + P) But the above does not work with Firedrake. Is there a way to do that, or is there something better? This kind of leads me to the second question: 2) I would like to try Q2-P1disc elements for Darcy problem (I see that this element combination is commonly used for stokes equations, but I don't know if it will work for Darcy). Is this possible to do? If the answer to the above two is no, is there another "locally conservative" mixed formulation that's currently doable other than the classical forms like RT0? And preferably with a higher order rate of convergence? Thanks, Justin