open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Everyone,

I am currently working on to build a script to convert a Selafin 3D file into a netcdf file on the Ugrid conventions. I notice that someone already did something closer to this, but only for 2D file (Coards convention), converting from unstructured into a regular grid.

My intention is to build a function which works also for 3D. I can load and access the data on Matlab and Python3, using the most recent python 3 scripts presented on the last TUC 2019 and can handle the interpolation for the first plane using IKLE2 (i.e. NPI,3) with triangulation function.

However, I am not understanding the IKLE3 variable, since it is something like N, 6...I was thinking to build a loop iterating on IKLE2 interval and reproduce it on the vertical layer, however the size of IKLE3 does not match IKLE2*NPLAN.

May anyone could help me on that? I mean, which would be the most straighforward way to do that?

Thanks in advance,

Regards,

Hi
IKLE3 should be (IKLE2*(NPLAN-1)) at least for prism structure.
IN 3D, between 2 planes, the element has 6 nodes. 3 for the triangle at the bottom and 3 for the triangle at the top...
As the node numbering is continuous, this means the numbering of the 3 top nodes is the numbering of the bottom ones plus NPOIN2...

Hope this helps