[Gmsh] Node label

[Philippe Ackerer]

Dear all,
Some troubles between node numbers in the screen 
(Tools->Option->Visibility->Mesh->Node Labels) and the node numbers in 
the .msh file.
I attached the .geo file.

On the screen, watch node number 41 (close to node 1). In the .msh file, 
it seems to be node number 39.

Could you tell me what I did wrong ?
Thanks

Best regards
Philippe

-- 
--------------------------------------------------------------
Philippe ACKERER
LHyGeS UMR 7517 UdS-CNRS-Engees
1, rue Blessig
67000 STRASBOURG
tel +33 (0)368 850 561
-------------------------------------------------------------

-------------- next part --------------

// As usual, we start by defining some variables:

r1=1.0;
r2=19.0;
Lcw=0.50;
Lc1 = 5.0;
Lc2 = 5.0;
dteta=2.0*3.14159/8.0;
teta=dteta;


// Then we define some points and some lines using these variables:

Point(1) = { 0.0, 0.0  , 0.0 , Lcw};
Point(2) = { r1, 0.0   , 0.0 , Lc1};
Point(3) = { r2 , 0.0   , 0.0 , Lc2};
x1=r2*Cos(teta);
y1=r2*Sin(teta);
Point(4) = { x1 , y1 , 0.0 , Lc2};
x1=r1*Cos(teta);
y1=r1*Sin(teta);
Point(5) = { x1 , y1 , 0.0 , Lc1};


// Circle(1) = {2,1,5}; starting point: 2, ending point: 5, center: 1.
Circle(1) = {2,1,5};
Circle(2) = {3,1,4};


Line(10) = {1,2};
Line(11) = {2,3};
Line(12) = {4,5};
Line(13) = {5,1};
 

Line Loop(21) = {10,1,13}; Plane Surface(31) = {21};
Line Loop(22) = {11,2,12,-1}; Plane Surface(32) = {22};

//Extrusion de la surface sur une distance h
h = 3.0;

out1[]=Extrude {0,0,h} {Surface{31,32}; Layers{ {3}, {1.00} };Recombine; };
Printf(" Top surface=%g", out1[0]);
Printf(" New volume=%g", out1[1]);
Printf(" Lateral surface=%g", out1[2]);
Printf(" Lateral surface=%g", out1[3]);
Printf(" Lateral surface=%g", out1[4]);
Printf(" Top surface=%g", out1[5]);
Printf(" New volume=%g", out1[6]);
Printf(" Lateral surface=%g", out1[7]);
Printf(" Lateral surface=%g", out1[8]);
Printf(" Lateral surface=%g", out1[9]);
Printf(" Lateral surface=%g", out1[10]);


Physical Surface(101) = {31};
Physical Surface(102) = {out1[0]};

Physical Volume(500) = {out1[2], out1[6]};