//Abmessungen
//Radius des Behälters
r = 6E-2; // war mal 9 cm
//Radius der Achse
t = 0.1E-2;
//Länge eines Flügels
l = 2E-2;
//Dicke eines Flügels
d = 2.5E-4;
//Abstand der Flügel vom Ursprung
a = 1E-2;
//Länge des Übergangs der Auflösungen auf dem Flügel
b=0.005;

// Punkte an denen sich der Spilne orientiert 
e = d/8;
s = b/20;

//Elementauflösung der Punkte
/*
vHighPr = 0.00002;
highPr =  0.00004;
midPr =   0.00015;
lowPr =   0.0006;
vLowPr =  0.0006;
*/

vHighPr = 0.000005;
//highPr =  0.00004;
midPr =   0.00015;
lowPr =   0.0004;
vLowPr =  0.001;

Geometry.NumSubEdges = 5000; // erhöht die sample rate der splines

Point(0) = {0,0,0,lowPr}; // Ursprung

//Flügel rechts
Point(1) = {a,-0.5*d,0,vHighPr};
Point(2) = {a+e,-0.5*d,0,vHighPr};
Point(3) = {a+s,-0.5*d,0,vHighPr};
Point(4) = {a+b,-0.5*d,0,midPr};
Point(5) = {a+l-b,-0.5*d,0,midPr};
Point(6) = {a+l-s,-0.5*d,0,vHighPr};
Point(7) = {a+l-e,-0.5*d,0,vHighPr};
Point(8) = {a+l,-0.5*d,0,vHighPr};
Point(9) = {a+l,-0.5*d+e,0,vHighPr};
Point(10) = {a+l,0.5*d-e,0,vHighPr};
Point(11) = {a+l,0.5*d,0,vHighPr};
Point(12) = {a+l-e,0.5*d,0,vHighPr};
Point(13) = {a+l-s,0.5*d,0,vHighPr};

Point(14) = {a+l-b,0.5*d,0,midPr};

Point(15) = {a+b,0.5*d,0,midPr};
Point(16) = {a+s,0.5*d,0,vHighPr};
Point(17) = {a+e,0.5*d,0,vHighPr};
Point(18) = {a,0.5*d,0,vHighPr};
Point(19) = {a,0.5*d-e,0,vHighPr};
Point(20) = {a,-0.5*d+e,0,vHighPr};

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

BSpline(4)={6,7,9,10,12,13};
Physical Line (2) = { 4 };

Line(5) = {13,14};
Line(6) = {14,15};
Line(7) = {15,16};
Physical Line (3) = { 5,6,7 };

BSpline(8)={16,17,19,20,2,3};
Physical Line (4) = { 8 };

Line Loop(1001) = {1,2,3,4,5,6,7,8};

//Flügel links
Point(21) = {-a,0.5*d,0,vHighPr};
Point(22) = {-(a+e),0.5*d,0,vHighPr};
Point(23) = {-(a+s),0.5*d,0,vHighPr};
Point(24) = {-(a+b),0.5*d,0,midPr};
Point(25) = {-(a+l-b),0.5*d,0,midPr};
Point(26) = {-(a+l-s),0.5*d,0,vHighPr};
Point(27) = {-(a+l-e),0.5*d,0,vHighPr};
Point(28) = {-(a+l),0.5*d,0,vHighPr};
Point(29) = {-(a+l),0.5*d-e,0,vHighPr};
Point(30) = {-(a+l),-0.5*d+e,0,vHighPr};
Point(31) = {-(a+l),-0.5*d,0,vHighPr};
Point(32) = {-(a+l-e),-0.5*d,0,vHighPr};
Point(33) = {-(a+l-s),-0.5*d,0,vHighPr};

Point(34) = {-(a+l-b),-0.5*d,0,midPr};

Point(35) = {-(a+b),-0.5*d,0,midPr};
Point(36) = {-(a+s),-0.5*d,0,vHighPr};
Point(37) = {-(a+e),-0.5*d,0,vHighPr};
Point(38) = {-a,-0.5*d,0,vHighPr};
Point(39) = {-a,-(0.5*d-e),0,vHighPr};
Point(40) = {-a,0.5*d-e,0,vHighPr};

Line(11) = {23,24};
Line(12) = {24,25};
Line(13) = {25,26};
Physical Line (5) = { 11,12,13 };

BSpline(14)={26,27,29,30,32,33};
Physical Line (6) = { 14 };

Line(15) = {33,34};
Line(16) = {34,35};
Line(17) = {35,36};
Physical Line (7) = { 15,16,17 };

BSpline(18)={36,37,39,40,22,23};
Physical Line (8) = { 18 };

Line Loop(1002) = {11,12,13,14,15,16,17,18};

//Flügel oben
Point(41) = {0.5*d,a,0,vHighPr};
Point(42) = {0.5*d,a+e,0,vHighPr};
Point(43) = {0.5*d,a+s,0,vHighPr};
Point(44) = {0.5*d,a+b,0,midPr};
Point(45) = {0.5*d,a+l-b,0,midPr};
Point(46) = {0.5*d,a+l-s,0,vHighPr};
Point(47) = {0.5*d,a+l-e,0,vHighPr};
Point(48) = {0.5*d,a+l,0,vHighPr};
Point(49) = {0.5*d-e,a+l,0,vHighPr};
Point(50) = {-0.5*d+e,a+l,0,vHighPr};
Point(51) = {-0.5*d,a+l,0,vHighPr};
Point(52) = {-0.5*d,a+l-e,0,vHighPr};
Point(53) = {-0.5*d,a+l-s,0,vHighPr};

Point(54) = {-0.5*d,a+l-b,0,midPr};

Point(55) = {-0.5*d,a+b,0,midPr};
Point(56) = {-0.5*d,a+s,0,vHighPr};
Point(57) = {-0.5*d,a+e,0,vHighPr};
Point(58) = {-0.5*d,a,0,vHighPr};
Point(59) = {-0.5*d+e,a,0,vHighPr};
Point(60) = {0.5*d-e,a,0,vHighPr};

Line(21) = {43,44};
Line(22) = {44,45};
Line(23) = {45,46};
Physical Line (9) = { 21,22,23 };

BSpline(24)={46,47,49,50,52,53};
Physical Line (10) = { 24 };

Line(25) = {53,54};
Line(26) = {54,55};
Line(27) = {55,56};
Physical Line (11) = { 25,26,27 };

BSpline(28)={56,57,59,60,42,43};
Physical Line (12) = { 8 };

Line Loop(1003) = {21,22,23,24,25,26,27,28};



//Flügel unten
Point(61) = {-0.5*d,-a,0,vHighPr};
Point(62) = {-0.5*d,-(a+e),0,vHighPr};
Point(63) = {-0.5*d,-(a+s),0,vHighPr};
Point(64) = {-0.5*d,-(a+b),0,midPr};
Point(65) = {-0.5*d,-(a+l-b),0,midPr};
Point(66) = {-0.5*d,-(a+l-s),0,vHighPr};
Point(67) = {-0.5*d,-(a+l-e),0,vHighPr};
Point(68) = {-0.5*d,-(a+l),0,vHighPr};
Point(69) = {-0.5*d+e,-(a+l),0,vHighPr};
Point(70) = {0.5*d-e,-(a+l),0,vHighPr};
Point(71) = {0.5*d,-(a+l),0,vHighPr};
Point(72) = {0.5*d,-(a+l-e),0,vHighPr};
Point(73) = {0.5*d,-(a+l-s),0,vHighPr};

Point(74) = {0.5*d,-(a+l-b),0,midPr};

Point(75) = {0.5*d,-(a+b),0,midPr};
Point(76) = {0.5*d,-(a+s),0,vHighPr};
Point(77) = {0.5*d,-(a+e),0,vHighPr};
Point(78) = {0.5*d,-a,0,vHighPr};
Point(79) = {0.5*d-e,-a,0,vHighPr};
Point(80) = {-0.5*d+e,-a,0,vHighPr};

Line(31) = {63,64};
Line(32) = {64,65};
Line(33) = {65,66};
Physical Line (13) = { 31,32,33 };

BSpline(34)={66,67,69,70,72,73};
Physical Line (14) = { 34 };

Line(35) = {73,74};
Line(36) = {74,75};
Line(37) = {75,76};
Physical Line (15) = { 35,36,37 };

BSpline(38)={76,77,79,80,62,63};
Physical Line (16) = { 38 };

Line Loop(1004) = {31,32,33,34,35,36,37,38};


//Außenwand

Point(101) = {-r,0,0,vLowPr};
Point(102) = {r,0,0,vLowPr};

Circle(101) = {101,0,102};
Circle(102) = {102,0,101};
Physical Line (100) = {101,202};

Line Loop(2001) = {101,102};


Point(103) = {-t,0,0,lowPr};
Point(104) = {t,0,0,lowPr};

Circle(103) = {104,0,103};
Circle(104) = {103,0,104};

Line Loop(2002) = {103,104};

Physical Line (101) = {103,104};

//Flächen
Plane Surface(1) = {2001,2002,1001,1002,1003,1004};

Physical Surface ( 1 ) = { 1003 };



