How to set initial field value by every point's position

[dokken]

Your question is very general, and would depend on the structure of the data from the other software.

In general, you can

1. Create a relation between the location of each degree of freedom and its physical coordinate in space by using FunctionSpace.tabulate_dof_coordinates()
2. You can then directly assign data to the underlying vector in 'Function.vector()`

See for instance

How to modify the value of a Function object?

To me, it seems like your solution (c1) before trying to use the conditional, is highly unstable (spurious behavior). So this is something you have to fix regardless of how you would like to bound your function. The following alternative works for bounding the variables: ff = File("c1.pvd") for n in range(num_steps): print("Step:", str(n)) # Compute solution for diffusion solve(a == L, c1, bc) ff << c1 indices = np.argwhere(c1.vector()[:] > 3.75) c1.vector().vec().set…