So, here is my situation:

I am trying to introduce noise into my initial conditions for a Turbulent Channel DNS using awgn(). Currently, it is set up along the lines of u = 1.0 + 0.1*awgn(1.0). However, for what we are doing, we need this initial noise to be symmetric over the XY-plane at the middle of the channel (i.e. the bottom half is identical to the top half).

I have been trying various combinations of modules in FieldConvert, but I feel as though I’m missing something. (If you could offer and insight on this, that would be great.)

My other idea was to bring the initial field into a program such as paraview, and then make the necessary edits there. However, after reading and tinkering, I don’t believe that FieldConvert has the ability to read in a VTK file and convert it back to a FLD file. Is there any way to actually read in data from a source such as VTK file or some other non-Nektar++ source? I tried this with a tecplot file, but it does not seem to want to read it.