Re: [firedrake] close to fixing my nonlinear QG solver

Hi Francis, sorry, I now saw your problem. You do: soln1.assign(soln0) psi1, u1, q1 = soln1.split() Gnon = F(psi1, u1, q1) You've been bitten by the "wrong split" bug!! ARGH. Only use function.split() for output. If you're going to use the split values in a form, you need to use "split(function)" If you replace your lines: psi1, u1, q1 = soln1.split() Gnon = ... with: psi1, u1, q1 = split(soln1) Gnon = ... psi1, u1, q1 = soln1.split() Everything works. Lawrence