lc = 0.03;

// Indica il numero elementi in cui è suddiviso ogni 0.1 elementi.
// Il numero di mesh è dato da:
// parallelepipedi = 180 * nm^2

nm = 5;

Point(1) = {0,0,0,lc};

l1[] = Extrude {0,0.5,0} {
  Point{1}; Layers{nm*5,1};
};

l2[] = Extrude {0,0.3,0} {
  Point{2}; Layers{nm*3,1};
};

l3[] = Extrude {-0.5,0,0} {
  Point{3}; Layers{nm*5,1};
};

l3[] = Extrude {0,0.2,0} {
  Point{4}; Layers{nm*2,1};
};

s1[] = Extrude {0.5,0,0} {
  Line{l3[1]}; Layers{nm*5,1}; Recombine;
};

s2[] = Extrude {0.5,0,0} {
  Line{l2[1]}; Layers{nm*5,1}; Recombine;
};

s3[] = Extrude {0.5,0,0} {
  Line{l1[1]}; Layers{nm*5,1}; Recombine;
};

s4[] = Extrude {0.1,0,0} {
  Line{13}; Layers{nm*1,1}; Recombine;
};

s5[] = Extrude {0.5,0,0} {
  Line{17}; Layers{nm*5,1}; Recombine;
};

s6[] = Extrude {0,0.5,0} {
  Line{23}; Layers{nm*5,1}; Recombine;
};

s7[] = Extrude {0.1,0,0} {
  Line{27}; Layers{nm*1,1}; Recombine;
};

s8[] = Extrude {0.5,0,0} {
  Line{29}; Layers{nm*5,1}; Recombine;
};

s9[] = Extrude {0,-0.3,0} {
  Line{34}; Layers{nm*3,1}; Recombine;
};

s10[] = Extrude {0,-0.2,0} {
  Line{37}; Layers{nm*2,1}; Recombine;
};

s11[] = Extrude {0.5,0,0} {
  Line{43}; Layers{nm*5,1}; Recombine;
};

s12[] = Extrude {0.5,0,0} {
  Line{5}; Layers{nm*5,1}; Recombine;
};

// unisco tutte le superfici in un unico dominio
stot[] = {s1[1], s2[1], s3[1], s4[1], s5[1], s6[1], s7[1], s8[1], s9[1], s10[1], s11[1], s12[1]};

// Estrudo il dominio
Extrude {0,0,0.1} {
  Surface{stot[]}; Layers{1,1}; Recombine;
}

// INLET
Physical Surface(291) = {61};

// OUTLET
Physical Surface(292) = {289};

// WALL
Physical Surface(317) = {73,83,105,117,139,161,65,307,311,91,131,183,179,197,219,223,241,249,205,157,271,267,285,293};

// TOP/BOTTOM
Physical Surface(318) = {162,24,184,28,206,32,228,36,250,40,272,44,294,48,16,118,12,96,52,316,74,8,140,20};

Physical Volume(1000) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};