# High order gmsh support

**URL:** <https://pyfr.discourse.group/t/high-order-gmsh-support/413>\
**Category:** Just Starting\
**Created:** [8 June 2021 10:03 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413 "2021-06-08T10:03:51Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![Zhenyang](https://avatars.discourse-cdn.com/v4/letter/z/ec9cab/32.png) [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Post date:** [8 June 2021 10:03 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/1 "2021-06-08T10:03:51Z")

</div>

Hi,

I am trying to import high order mesh generated by Gmsh to Pyfr. However, some of the high order mesh element type (i.e type 28, fifth order edge) is not supported. Do you know how can I process this kind of mesh and import to pyfr? And if only lower order mesh is supported in the pyfr (i.e, order 3), can I use higher p order in the ini file without any interpolation pre-running?

Thanks ahead!

Best,  
Zhenyang

---

<div class="post-metadata">

**Author:** ![WillT](https://avatars.discourse-cdn.com/v4/letter/w/bc79bd/32.png) [@WillT](https://pyfr.discourse.group/u/WillT)\
**Post date:** [9 June 2021 14:05 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/2 "2021-06-09T14:05:23Z")

</div>

Are you trying to generate a 2D or 3D mesh?

---

<div class="post-metadata">

**Author:** ![Zhenyang](https://avatars.discourse-cdn.com/v4/letter/z/ec9cab/32.png) [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Post date:** [9 June 2021 14:16 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/3 "2021-06-09T14:16:40Z")

</div>

Hi,

I only tried to generate a 2D mesh.

---

<div class="post-metadata">

**Author:** ![WillT](https://avatars.discourse-cdn.com/v4/letter/w/bc79bd/32.png) [@WillT](https://pyfr.discourse.group/u/WillT)\
**Post date:** [9 June 2021 14:19 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/4 "2021-06-09T14:19:14Z")

</div>

When you say it is not supported, what is the error message and traceback?

---

<div class="post-metadata">

**Author:** ![fdw](https://avatars.discourse-cdn.com/v4/letter/f/a587f6/32.png) [@fdw](https://pyfr.discourse.group/u/fdw)\
**Post date:** [9 June 2021 15:07 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/5 "2021-06-09T15:07:54Z")

</div>

The error indicates that your mesh has 1D elements (edges). These should not be present in a 2D mesh (which should contain exclusively triangles and quadrilaterals).

Regards, Freddie.

---

<div class="post-metadata">

**Author:** ![Zhenyang](https://avatars.discourse-cdn.com/v4/letter/z/ec9cab/32.png) [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Post date:** [9 June 2021 15:13 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/6 "2021-06-09T15:13:47Z")

</div>

Hi,

For order 6 edge: “ValueError: Unsupported element type 62”

For order 5 edge: “ValueError: Unsupported element type 28”

element type can be found here: [Gmsh 4.8.4](https://gmsh.info/doc/texinfo/gmsh.html#MSH-file-format)

---

<div class="post-metadata">

**Author:** ![Zhenyang](https://avatars.discourse-cdn.com/v4/letter/z/ec9cab/32.png) [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Post date:** [9 June 2021 15:26 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/7 "2021-06-09T15:26:56Z")

</div>

Hi Freddie,

I guess from what generated by  
gmsh, they always list element from 1D to higher dimension in ‘$Element’ section. And for the lower order like order 3, it can be read by pyfr. I attach two test mesh files, only order is changed the rest is the same.

> **[z-yuan131/Pyfr\_test](https://github.com/z-yuan131/Pyfr_test)**
>
> Contribute to z-yuan131/Pyfr\_test development by creating an account on GitHub.

---

<div class="post-metadata">

**Author:** ![Zhenyang](https://avatars.discourse-cdn.com/v4/letter/z/ec9cab/32.png) [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Post date:** [9 June 2021 15:29 UTC](https://pyfr.discourse.group/t/high-order-gmsh-support/413/8 "2021-06-09T15:29:10Z")

</div>

Sorry forgot to attach the traceback:

Traceback (most recent call last):  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/bin/pyfr”, line 33, in   
sys.exit(load\_entry\_point(‘pyfr==1.11.0’, ‘console\_scripts’, ‘pyfr’)())  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/ **main**.py”, line 111, in main  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/ **main**.py”, line 122, in process\_import  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/readers/ **init**.py”, line 18, in get\_reader\_by\_extn  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/readers/gmsh.py”, line 89, in **init**  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/readers/gmsh.py”, line 222, in \_read\_eles  
File “/cfs/klemming/nobackup/z/zhenyang/conda-dirs/envs/pyfr-tegner/lib/python3.7/site-packages/pyfr-1.11.0-py3.7.egg/pyfr/readers/gmsh.py”, line 253, in \_read\_eles\_impl\_v41  
ValueError: Unsupported element type 28
