physical lines changed after boolean operators
Dear mister/miss,
I ran into the following problem: I am constructing a fairly large mesh of a river delta where I have 5 types of physical entities: 4 lines (TidalEntrance, River, closed, closedInner) and one physical surface (covering the entire domain). These will be boundaries with varying boundary conditions when I use the mesh in a hydrodynamic model (TELEMAC).
There are large islands in my domain. If I don't exclude them (no boolean operators) the physical lines remain perfectly in place where I have instructed them to be ( Physical Line(1) = {1, 2, 3} where 1, 2 and 3 are lines (Line(1) = {...}). However, after I apply the boolean operators, they shift and the physical labels change and other line segments now get the TidalEntrance label while the original ones now have the 'closed' label.
When I use the visibility or labelling tool in the gmsh gui, I see that the original curve labels have changed. Originally, the labels start at my node 1 where curve labeled 1 started and clockwise, all curves where labelled +1. After applying the boolean operators, the labeleling starts again at node 1 but goes 2, 4, 6, ... clockwise and 3,5,7,... counterclockwise. As a consequence, the physical line defined by curves 1, 2 and 3 is now located with node 1 in the middle.
All of this happens without meshing. I made the .geo file myself, based on a python program. I added the .geo file.
Kind regards, Ignace Pelckmans PhD-fellow, UAntwerp, department of BiologyGulf.geo