Hi all, I am trying to run a density driven flow problem where the computational domain has some initial perturbation in the concentration field. I have a unit cube domain and I want the concentrations at z >= 0.9 to be randomly perturbed by some random values. So I have a text file 'initialfield.txt' containing these perturbations which looks something like this: <xcoord> <ycoord> <zcoord> <value1> <xcoord> <ycoord> <zcoord> <value 2> ... Which is for a specific coarse grid resolution, say h-size = 1/20. How would I best go about projecting this data into a Firedrake/UFL Function? Does the initialfield.txt data set have to have the same connectivity/mapping as the Mesh function? Also, if I wanted to project this data onto a finer grid, like h-size = 1/40 or 1/60, would I create two meshes, one coarse mesh for the initialfield.txt data and one fine grid that the data shall be projected onto? Any guidelines or tips on how to best do this appreciated, thanks! Justin