Dear All, How do I instruct Firedrake about the dimension of the problem when using external mesh from gmsh? I had a working code with utility mesh. I changed: mesh = RectangleMesh(nx, nz, Lx, Lz) V = VectorFunctionSpace(mesh, "CG", 1) to: mesh = Mesh("mesh.msh") V = VectorFunctionSpace(mesh, "CG", 1, dim=2) yet I get the error at further expressions that worked previously: ... File "/home/mmtjs/work/programs/lin_coupled_fd_2d/lib/beam.py", line 73, in initialize_solvers + self.mu * ( fd.inner( delX, delv_B + fd.transpose(delv_B) ) ) File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/exproperators.py", line 177, in _add return Sum(self, o) File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/algebra.py", line 55, in __new__ error("Can't add expressions with different shapes.") File "/opt/firedrake/1510/lib/python2.7/site-packages/ufl/log.py", line 151, in error raise self._exception_type(self._format_raw(*message)) ufl.log.UFLException: Can't add expressions with different shapes. Thank you, Tomasz