Re: [firedrake] how to set up random initial conditions?

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