[Gmsh] problème extrusion ligne-physical surface

[Christophe Geuzaine]

SALANEUVE Laurent wrote:
> 
> J'effectue une extrusion de ligne et veux récupérer la surface créée afin de
> l'affecter à une physical surface.
> J'écris ainsi:
> 
> test_surf[p]=newreg;
> Extrude Line{l1, {0,0,8}} {
> Layers { {nz}, {test_surf[p]}, {1} };
> };
> 
> Physical Surface(22222)={test_surf[p]};
> 
> Or, lorsque l'on affiche la surface créée (test_surf[p]), en fait la
> référence correspond à une ligne.
> Ce qui m'étonne c'est que la physical surface est tout de même créée et
> renvoie la ligne en question.

OK. This is because, when you use the second argument in Layers to 
assign the surface number, you lose the 1-to-1 relation between surfaces 
and surface numbers (you directly assign a number to each element in the 
layer). This technique is actually only useful if you have more than one 
layer.

The best solution is to set the second argument to zero, or to remove it 
altogether. In that case, the surface number is generated automatically, 
and you can retrieve it with

p1[]=Extrude Line{l1, {0,0,8}} {
  Layers { {nz}, {1} };
};;


Best,

Christophe

-- 
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org