//Abmessungen
//Radius des Behälters
r = 6E-2; // war mal 9 cm
//Länge eines Flügels
l = 2.32E-2;
//Dicke eines Flügels
d = 2.5E-4;
//Abstand der Flügel vom Ursprung
a = 0.94E-2;

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

vHighPr = 0.000018;
highPr =  0.00004;
midPr =   0.0002;
lowPr =   0.0006;
vLowPr =  0.001;

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

//rechter Flügel
Point(1) = {a,-0.5*d,0,midPr};
Point(2) = {a+l,-0.5*d,0,midPr};
Point(3) = {a+l,0.5*d,0,midPr};
Point(4) = {a,0.5*d,0,midPr};

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

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

//linker Flügel
Point(5) = {-a,0.5*d,0,midPr};
Point(6) = {-(a+l),0.5*d,0,midPr};
Point(7) = {-(a+l),-0.5*d,0,midPr};
Point(8) = {-a,-0.5*d,0,midPr};

Line(5) = {5,6};
Line(6) = {6,7};
Line(7) = {7,8};
Line(8) = {8,5};
Physical Line (5) = {5};
Physical Line (6) = {6};
Physical Line (7) = {7};
Physical Line (8) = {8};

Line Loop(102) = {5,6,7,8};

//oberer Flügel
Point(9) = {-0.5*d,a,0,midPr};
Point(10) = {-0.5*d,a+l,0,midPr};
Point(11) = {0.5*d,a+l,0,midPr};
Point(12) = {0.5*d,a,0,midPr};

Line(9) = {9,10};
Line(10) = {10,11};
Line(11) = {11,12};
Line(12) = {12,9};

Physical Line (9) = {9};
Physical Line (10) = {10};
Physical Line (11) = {11};
Physical Line (12) = {12};

Line Loop(103) = {9,10,11,12};

//unterer Flügel
Point(13) = {0.5*d,-a,0,midPr};
Point(14) = {0.5*d,-(a+l),0,midPr};
Point(15) = {-0.5*d,-(a+l),0,midPr};
Point(16) = {-0.5*d,-a,0,midPr};

Line(13) = {13,14};
Line(14) = {14,15};
Line(15) = {15,16};
Line(16) = {16,13};

Physical Line (13) = {13};
Physical Line (14) = {14};
Physical Line (15) = {15};
Physical Line (16) = {16};

Line Loop(104) = {13,14,15,16};

//Außenwand
Point(100) = {-r,0,0,vLowPr};
Point(200) = {r,0,0,vLowPr};

Circle(100) = {100,0,200};
Circle(200) = {200,0,100};

Physical Line (17) = {100,200};

Line Loop(105) = {100,200};

Plane Surface(106) = {101,102,103,104,105};

//lables for Elmer
Physical Surface ( 1 ) = { 106 };

//lokale Auflösung beeinflussen

//feine Auflösung für die Flügelecken
Field[1] = Attractor;
Field[1].EdgesList = {2,4,6,8,10,12,14,16};

Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = vHighPr;
Field[2].LcMax = vLowPr;
Field[2].DistMin = 0.025E-2;
Field[2].DistMax = 0.25E-2;

Field[5] = Attractor;
Field[5].EdgesList = {2,4,6,8,10,12,14,16};

Field[6] = Threshold;
Field[6].IField = 5;
Field[6].LcMin = highPr;
Field[6].LcMax = vLowPr;
Field[6].DistMin = 0.1E-2;
Field[6].DistMax = 1E-2;

// Auflösung der Flügel
Field[3] = Attractor;
Field[3].EdgesList = {1,3,5,7,9,11,13,15};

Field[4] = Threshold;
Field[4].IField = 3;
Field[4].LcMin = midPr;
Field[4].LcMax = vLowPr;
Field[4].DistMin = 0.5E-2;
Field[4].DistMax = 2E-2;


//Backround
Field[10] = Min;
Field[10].FieldsList = {2,4,6};

Background Field = 10;


