
lc = 0.09;

Point(1) = {0, 0, 0, lc};
Point(2) = {.1, 0,  0, lc} ;
Point(3) = {.1, .3, 0, lc} ;
Point(4) = {0,  .3, 0, lc} ;

Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;

Line Loop(5) = {1,-2,3,4} ;

Plane Surface(6) = {5} ;
 
Transfinite Line{1} = 4;
Transfinite Line{2} = 4;
Transfinite Line{3} = 4;
Transfinite Line{4} = 4;

Transfinite Surface{6} = {1,2,3,4};

Recombine Surface{6};

// Volume 1 will be structured with HEX elements
extrudetovol1[] = Extrude {0.0, 0.0, 0.5} {Surface {6}; Layers {{1,1,2},{0.1,0.2,1}}; Recombine;};


surfOpposite = extrudetovol1[0];

// Extrude surfOpposite (= # 28) to Volume 2 filled mainly with TETs but with
// 2 layers of pyramids near surfOpposite
Extrude {0.0, 0.0, 0.7} {Surface {surfOpposite}; Layers {8}; QuadTriDbl; Recombine;}


// Note: You can achieve worse element quality with QuadTriDSngl instead of QuadTriDbl

