//

Mesh.SubdivisionAlgorithm=2;
          //From gmsh userguide:
	  //Mesh subdivision algorithm (0=none, 1=all quadrangles, 2=all hexahedra)
          //Default value: 0
          //Saved in: General.OptionsFileName


lc = 1;
lca = 0.01;

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

Point(5) = {0, 0, 0, lca} ;
Point(6) = {0.25, 0.01, 0, lca} ;
Point(7) = {0.5, 0, 0, lca} ;
Point(8) = {0.325, 0.075, 0, lca} ;

Circle (303) = {1, 5, 2 };
Circle (304) = {2, 5, 3 };
Circle (305) = {3, 5, 4 };
Circle (306) = {4, 5, 1 };



Spline(101) = {6, 7, 8, 5, 6};

Line Loop(5) = {303,304,305,306} ;
Line Loop(6) = {101} ;


Plane Surface(9) = {5,6} ;

Transfinite Line{101} = 40;
Transfinite Line{303,304,305,306} = 10;
Transfinite Surface{9} = {1,2,3,4};
Recombine Surface{9};

Extrude {0,0,0.001} {
  Surface{9}; Layers{1}; Recombine;
}


Physical Surface ("airfoil") = {332};
Physical Surface ("Farfield") = {316,320,324,328};

Physical Volume ("InternalField") = {1};

Physical Surface("FrontAndBack") = {6,333};
