NadaSulaiman

- The **orthogonality** is defined as the cosine of the angle **φ** between a **flowlink** (a line segment connecting two cell circumcenters) and a **netlink** (an edge of a cell,...

Hi all, many thanks for the efforts. I have adjusted the code to obtain identical orthogonality values as it was displayed inside the DELFT FM grid generation module ([RGFGRID](https://content.oss.deltares.nl/delft3d/manuals/RGFGRID_User_Manual.pdf)). Summary...

> > I am attempting to compute the orthogonality of the New Zealand msh obj (acquired by running the example 1), however I am unsure which edge of the triangle...

Alternatively, instead of saving the mesh as .14. and then converting it to _net.nc, we could write it directly as DELFT mesh by adding those three lines to the end...

Since `C = circumcenter(TR) ` does not equal ``` C(:,1)=ncread('delftexport_net.nc','mesh2d_face_x')'; C(:,2)=ncread('delftexport_net.nc','mesh2d_face_y')'; ``` Could anyone have a quick look at the definition of how the circumcenter of a triangle is computed...

Hi Keith. The grid is uploaded in the zip file here [[delftexport_net.zip](https://github.com/CHLNDDEV/OceanMesh2D/files/8598371/delftexport_net.zip)]. You should be able to run the code with it.

You are correct. I think it places it to the edge closest to the location of the circumcenter outside the element. Suppose we adjust the code above to get the...