ElmerGrid -increase generates elementtype 0 and elmersolver then fails

[tomatoma]

Hi,

I'd like to export a .unv mesh from Salome, to then convert it to Elmer format using ElmerGrid.
The mesh contains HEXA8, TETRA4 and PYRA5.
Salome cannot export PYRA5 elements to .unv mesh.
So I exported the mesh from Salome to Elmer by using the following Python script:
https://github.com/jvencels/salomeToElmer
as explained here:
https://www.elmerfem.org/forum/viewtopi ... =elmergrid

That script does not support second order elements, but I'd like to use them.
So I used the following command to increase the order of the elements:
ElmerGrid 2 2 Mesh_name -increase

The problem is that when ElmerGrid reaches the following step:
Saving boundary elements to mesh.boundary.
It writes 542 times the following message:
GetElementDimension: unknown elementtype 0
and in the file mesh.boundary there are indeed 542 rows with elements of type 0.
1 18 6429 0 0
2 18 6495 0 0
3 18 6498 0 0
...
542 18 7160 0 0
543 18 12428 0 203 9115 9117 56248

elmersolver then fails as follow:
...
ElmerAsciiMesh: Base mesh name: ./.
WARNING:: GetElementType: Element type code 0 not found. Ignoring element.
ERROR:: ReadBoundaryFile: Element of type 0could not be associated!
STOP 1

(no problems with the original first order mesh)

Can anyone please help?
Is there any way to find those type-0 elements in the mesh/geometry?

Here is the original first-order mesh obtained from Salome with the salomeToElmer script:

first_order.7z

(604.28 KiB) Downloaded 7 times

Elmergrid version:
Starting program Elmergrid, compiled on Apr 26 2024

Thanks

[kevinarden]

If you do
ElmerGrid 2 2 mesh -autoclean -increase

it works

elmermesh.png (32.53 KiB) Viewed 66 times

[tomatoma]

Thanks!
That was easy!