/********************************************************************* 
 *
 *  Gmsh tutorial 1
 * 
 *  Variables, Elementary entities (Points, Lines, Surfaces), Physical
 *  entities (Points, Lines, Surfaces), Background mesh
 *
 *********************************************************************/

lc1= 2.0;
lc2 = 1.0 ;
pi = 3.1416;
ang1 = pi/20. ; DL= 100.0; DH = 30.0;  z1 = Tan(ang1)*DL; width = 5.0;     
Dx1 = DL; Dz1 = 3.0; 
Dx2 = 95; Dz2 = 10.0; 
Dx3 = 90.0; Dz3 = DH;
Dx4 = 70.0; Dz4 = DH;
Dx5 = Dx4;  Dz5 = DH - 3.0;
Dx6 = Dx4 - width; Dz6 = Dz5;
Dx7 = Dx6; Dz7 = DH;
Dx8 = 0.0; Dz8 = DH;
Dx9 = 0.0; Dz9 = z1 + 5;
Dx10 = 0.0 ;  Dz10 = z1;

Point(1) = {DL,  0.,  0, lc1} ;

// Three additional points are then defined:

Point(2) = {Dx1, Dz1, 0, lc1} ;
Point(3) = {Dx2, Dz2, 0, lc2} ;
Point(4) = {Dx3, Dz3, 0, lc2} ;
Point(5) = {Dx4, Dz4, 0, lc2} ;
Point(6) = {Dx5, Dz5, 0, lc2} ;
Point(7) = {Dx6, Dz6, 0, lc2} ;
Point(8) = {Dx7, Dz7,  0, lc2} ;
Point(9) = {Dx8, Dz8, 0, lc2} ;
Point(10)= {Dx9, Dz9, 0, lc1} ;
Point(11)= {Dx10, Dz10, 0, lc1} ;




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

Line Loop (13) = {2,3,4,5,6,7,8,9,-12};

Plane Surface (14) = {13};


Line Loop (15) = {1,12, 10, 11};
Plane Surface(16) = {13,15};

Physical Surface(1) = {13};
Physical Surface(2) = {15};

//Physical Point(1) = {1,2} ;

//Physical Line(10) = {1,2,4} ;
//MySurface = 100;
//Physical Surface(MySurface) = {6} ;

