Ah, I see it's not the sphere bit that's the problem - it's the Interpolator. Just to confirm - I did follow the instructions on the website so I should have the right branches of tsfc and PyOP2. J ________________________________ From: firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of Shipton, Jemma <j.shipton@imperial.ac.uk> Sent: 04 June 2019 21:47 To: firedrake Subject: [firedrake] complex problem with sphere meshes Hi All, The following used to work on the complex branch but now throws an error.... mesh=IcosahedralSphereMesh(1, 3) x=SpatialCoordinate(mesh) mesh.init_cell_orientations(x) The error is: File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/mesh.py", line 1216, in init_cell_orientations f.interpolate(ufl.dot(expr, ufl.cross(ReferenceGrad(x)[:, 0], ReferenceGrad(x)[:, 1]))) File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/function.py", line 334, in interpolate return interpolation.interpolate(expression, self, subset=subset) File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/interpolation.py", line 37, in interpolate return Interpolator(expr, V, subset=subset).interpolate() File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/interpolation.py", line 57, in __init__ self.callable = make_interpolator(expr, V, subset) File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/interpolation.py", line 112, in make_interpolator loops.extend(_interpolator(V, f.dat, expr, subset)) File "/Users/jshipton/complex-firedrake/src/firedrake/firedrake/interpolation.py", line 171, in _interpolator ast, oriented, needs_cell_sizes, coefficients = compile_ufl_kernel(expr, to_pts, coords, parameters=parameters) ValueError: too many values to unpack (expected 4) Any help would be appreciated! Thanks! Jemma