The number of elements etc. seem to still be correct. Ive been looking through the thread that you sent, however, it seems to be about gmsh .msh files only. I have a .msh file from ansys icem. So I do not know if I can write the boundary parts like so:
meshio.write(“mesh.xdmf”, meshio.Mesh(points=msh.points, cells={“tetra”: msh.cells[“tetra”]}))
meshio.write(“mf.xdmf”, meshio.Mesh(points=msh.points, cells={“triangle”: msh.cells[“triangle”]},
cell_data={“triangle”: {“name_to_read”: msh.cell_data[“triangle”][“gmsh:physical”]}}))
when I try and run my code this part gives me an error .
boundaries = MeshFunction(“size_t”, mesh, 2) # Read the boundary indicators
xdmf.read(boundaries, ‘boundaries’)
ds = Measure(“ds”)(subdomain_data=boundaries)
ds_wall = ds(wall[0]) + ds(wall[1])
ds_veinwall = ds(wall[0])
this gives
File “/home/jack/anaconda2/envs/fenicsproject/lib/python3.7/site-packages/oasis-2018.1-py3.7.egg/oasis/problems/NSfracStep/simulation_Jack_3D.py”, line 183, in
xdmf.read(boundaries, ‘boundaries’)
RuntimeError:
*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at***
*** fenics-support@googlegroups.com***
*** Remember to include the error message listed below and, if possible,
*** include a minimal running example to reproduce the error.***
*** -------------------------------------------------------------------------
*** Error: Unable to open MeshFunction for reading.
*** Reason: Mesh Grid with data Attribute not found in XDMF.
*** Where: This error was encountered inside XDMFFile.cpp.
*** Process: 0*** *** DOLFIN version: 2019.1.0
*** Git changeset: f3de96341c54e434c48115de5a3975b482aee933
*** -------------------------------------------------------------------------
Any advice?
Cheers,
Jack