I was thinking of taking my Quasi-Geostrophic demo that I recently submitted and change the initial conditions to be a random field.  Then we can see the evolution of freely evolving 2D turbulence.  Unfortunately, I am not sure how to pick random initial conditions in firedrake.  

I found issue 514 that suggested doing the following:

import numpy as np 

perm_array = np.random.randn(Nx, Ny, Nz)

I guess this works if you have a regular mesh.  But if I have an irregular mesh then is there a way to define this for the mesh? 

Cheers, Francis