newlines in output files
@@ -263,24 +263,25 @@ def write_buffer(f, mesh, point_format, cell_format):
print(
"""# Generated by meshio {}
mesh3d
+
dimension
3
+
geomtype
0
-
""".format(
__version__
),
file=f,
)
- print("# surfnr bcnr domin domout np p1 p2 p3", file=f)
+ print("\n# surfnr bcnr domin domout np p1 p2 p3", file=f)
print("surfaceelements", file=f)
print(cells_per_dim[2], file=f)
for block in mesh.cells:
if _topological_dimension[block.type] == 2:
_write_cells(f, block)
- print("# matnr np p1 p2 p3 p4", file=f)
+ print("\n# matnr np p1 p2 p3 p4", file=f)
print("volumeelements", file=f)
print(cells_per_dim[3], file=f)
for block in mesh.cells:
@@ -288,7 +289,7 @@ geomtype
_write_cells(f, block)
print(
- "# surfid 0 p1 p2 trignum1 trignum2 domin/surfnr1 domout/surfnr2 ednr1 dist1 ednr2 dist2",
+ "\n# surfid 0 p1 p2 trignum1 trignum2 domin/surfnr1 domout/surfnr2 ednr1 dist1 ednr2 dist2",
file=f,
)
print("edgesegmentsgi2", file=f)
@@ -297,7 +298,7 @@ geomtype
if _topological_dimension[block.type] == 1:
_write_cells(f, block)
- print("# X Y Z", file=f)
+ print("\n# X Y Z", file=f)
print("points", file=f)
print(len(mesh.points), file=f)
@@ -309,9 +310,8 @@ geomtype
else:
p = points
np.savetxt(f, p, "% 20.16f")
- f.write("\n\n")
- print("# pnum index", file=f)
+ print("\n# pnum index", file=f)
print("pointelements", file=f)
print("0", file=f)
@@ -321,7 +321,6 @@ materials
1
1 domain
-
bcnames
1
1 boundary
