Defining boundary conditions in COMSOL Script

[razlin]

I'm using Comsol script to solve coefficient form PDE. I assign boundary conditions using e.g. fem.bnd.h={1 0 0 0} and fem.bnd.r={1 0 0 0}. However, depending on relative positions of various structures, the Comsol arbitrarily changes the numbering order of the boundaries in my geometry, which makes it impossible to automatically program geometry changes using the script. For instance, the 6 boundaries of the same cylinder could be numbered as 12,13,14,15,22,23 when it's located in one position but then they are changed to 8,9,10,11,12,13 when the cylinder moves to another position. Is there a way to make the boundary numbering static or at least know in what order it is assigned each time ?