neg_x=-1;
pos_x=2;
neg_y=-1;
pos_y=1;
z=0.6;
outer_scale=0.09;
refine_1 = 0.03;


//Below describes a random 2D geometry and links all those points
Point(1)= {-0.014548,-0.057256,0.0};
Point(2)= {0.04967,-0.01702,0.0};
Point(3)= {0.020572,0.042453000000000005,0.0};
Point(4)= {-0.03586400000000001,0.024861,0.0};
Point(5)= {-0.02998,-0.02998,0.0};
Line(1)={1,2};
Line(2)={2,3};
Line(3)={3,4};
Line(4)={4,5};
Line(5)={5,1};
Characteristic Length {1:5}=refine_1;

//An outer rectangular shape which describes the computation domain 
Point(1001)={neg_x,neg_y,0,outer_scale};
Point(1002)={pos_x,neg_y,0,outer_scale};
Point(1003)={pos_x,pos_y,0,outer_scale};
Point(1004)={neg_x,pos_y,0,outer_scale};
Line(1001)={1001,1002};
Line(1002)={1002,1003};
Line(1003)={1003,1004};
Line(1004)={1004,1001};
Curve Loop(1) = {1:5};
Curve Loop(2) = {1001:1004};
Plane Surface(1) = {1,2};

//Extrusion
Extrude {0,0,0.6}{Surface{1};}

//mesh refinement
Field[1] = Cylinder;
Field[1].Radius = 0.6;
Field[1].VIn = refine_1;
Field[1].VOut = outer_scale;
Field[1].XCenter = 0;
Field[1].YCenter = 0;
Field[1].ZAxis = 0.8;
Background Field = 1;

Physical Volume("fluid") = {1};
Physical Surface("east")={1196};
Physical Surface("west")={1188};
Physical Surface("top")={1192};
Physical Surface("bottom")={1200};
Physical Surface("periodic_0_l")={1};
Physical Surface("periodic_0_r")={1201};
Physical Surface("concaveTower")={1048,1052,1056,1060,1064,1068,1072,1076,1080,1084,1088,1092,1096,1100,1104,1108,1112,1116,1120,1124,1128,1132,1136,1140,1144,1148,1152,1156,1160,1164,1168,1172,1176,1180,1184};

