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?