Well, trying to generate the mesh using the create_mesh.py code, generated the following error:
File "/Users/xyz/opt/anaconda3/envs/fenics2/lib/python3.8/site-packages/pygmsh/built_in/geometry.py", line 522, in <listcomp>
nodes = list (set ([p.id for and in edges_list for p in e.points]))
AttributeError: 'EllipseArc' object has no attribute 'points'
A branch with updated mesh generation can be found here. Note that it does no longer depend on pygmsh, but directly on the python API. See the following Pull request for updates on when this is in master.
Thanks, it worked, but the names of the facets marker values must be changed in the code stokes_problem.py
line 99: from create_mesh import inflow_marker, outflow_marker, wall_marker, obstacle_marker, c_x, c_y, L, H
and the other lines that call those names.