// Declaration de la taille des elements triangulaires
lc = 0.5;  // il est possible d'en spécifier plusieurs ...

// Declaration des points du contour : coord. X, coord. Y, coord. Z, taille des elements sur le contour
Point(1) = {0, 0, 0, lc};
Point(2) = {1., 0,  0, lc} ;
Point(3) = {1., 0.66, 0, lc};
Point(4) = {1., 1., 0, lc} ;
Point(5) = {0,  1., 0, lc} ;
Point(6) = {0,  0.33, 0, lc} ;

// Construction des lignes du contour en "reliant" les points
Line(1) = {1,2} ;
Line(2) = {2,3} ;
Line(3) = {3,4} ;  // sortie du fluide
Line(4) = {4,5} ;
Line(5) = {5,6} ;
Line(6) = {6,1} ;  // entree du fluide

// Connection des lignes
Line Loop(1) = {1,2,3,4,5,6} ; // Si domaine exterieur, lecture dans le sens trigo
// attention au sens de lecture sinon numerotation inversee !!!

// Creation de la surface a partir des "line loop" declarees
Plane Surface(1) = {1} ;

// Identification des elements de surface par un numero logique 31
Physical Surface(31) = {1} ;

