Hello again,

Thanks Colin, that helps a great deal.

q0 = Function(Vdg).interpolate(Expression("0.0"))

q0.dat.data[:] = 0.01*np.random.randn(*q0.dat.shape)

As I mentioned I want to put together an example that shows 2D turbulence, and this will help a great deal. I also plan to put together an example on the barotropic instability of a jet. If any of these are of interest I would be happy to try and put together a demo or something.

Francis