Mark regions from mesh of Gmsh

mesh
4

I guess your are looking for something similar to: Dolfinx discontinous expression

As long as you can access the underlying PETScVector in the dolfin.Function, you can follow this thread