//
// 



Merge "Kugelschuttung.nas" ;

//Cylinder

h = 7;
x = 0.000000;
y = 0.000000;
z = -1.000000;
r = 1.2;

p1 = newp;  Point(p1) = {x,y,z,0.5};
p2 = newp;  Point(p2) = {x+r,y,z,0.5}; 
p3 = newp;  Point(p3) = {x,y+r,z,0.5}; 
p4 = newp;  Point(p4) = {-(x+r),y,z,0.5}; 
p5 = newp;  Point(p5) = {0,-(y+r),z,0.5}; 

c1 = newreg;  Circle(c1) = {p2,p1,p3}; 
c2 = newreg;  Circle(c2) = {p3,p1,p4}; 
c3 = newreg;  Circle(c3) = {p4,p1,p5}; 
c4 = newreg;  Circle(c4) = {p5,p1,p2}; 

l1  = newreg;  Line Loop(l1) = {c1,c2,c3,c4}; 
ps1 = newreg;  Plane Surface(ps1) = {l1}; 

cyl[]=Extrude {0,0,h} { Surface{ps1}; }; 

Physical Volume(99)={cyl[1]}; 


Surface Loop(1) = {l1};

Surface Loop(2) = { cyl[0]};

Volume(2)={2, 1};

Physical Volume("Fluid")  = {2};

Physical Surface("Bead  ") = {1};

//Physical Surface("Inlet") = {28};
//Physical Surface("Outlet") = {ps1};

Physical Surface("cylinder") = {2};



