//parameter of the breast
//R = 0.06;
//H = 0.07;
//lc = 0.15*R;
s = 0.0;
r = 0.5 * (R*R + H*H) / H;
shift = r-H;


MINUS_THREE_DIM = 2147483648;
MINUS_TWO_DIM   = 1073741824;
MINUS_ONE_DIM   =  536870912;

bc_x = 0.0; bc_y = 0.0; bc_z = 0.0;
bup_x = 0.0; bup_y = 0.0; bup_z = R;
bdown_x = 0.0; bdown_y = 0.0; bdown_z = -R;
br_x = R; br_y = 0.0; br_z = 0.0;
bl_x = -R; bl_y = 0.0; bl_z = 0.0;
Point(1) = {bc_x, bc_y, bc_z, lc};
Point(2) = {bup_x, bup_y, bup_z, lc};
Point(3) = {bdown_x, bdown_y, bdown_z, lc};
Point(5) = {br_x, br_y, br_z, lc};
Point(6) = {bl_x, bl_y, bl_z, lc};

cc_x = 0.0; cc_y = shift; cc_z = 0.0;
Point(4) = {cc_x, cc_y, cc_z, lc}; //center of the breast

fc_x = 0.0; fc_z = s; fc_y = -Sqrt(r*r - s*s)+shift;
Point(7) = {fc_x, fc_y, fc_z, lc};


Circle(1) = {2,1,5};
Circle(2) = {5,1,3};
Circle(3) = {3,1,6};
Circle(4) = {6,1,2};
Circle(5) = {2,4,7};
Circle(6) = {7,4,3};
Circle(7) = {6,4,7};
Circle(8) = {7,4,5};

Line Loop(1) = {1,2,3,4};
Ruled Surface(1) = {1};
Line Loop(2) = {5,8,-1};
Ruled Surface(2) = {2};
Line Loop(3) = {6,-2,-8};
Ruled Surface(3) = {3};
Line Loop(4) = {5,-7,4};
Ruled Surface(4) = {4};
Line Loop(5) = {6,3,7};
Ruled Surface(5) = {5};
Surface Loop(1) = {1,2,3,4,5};
Volume(1) = {1};

Physical Surface(2) = {1};
Physical Surface(2+MINUS_TWO_DIM) = {1};

Physical Surface(3) = {2,3,4,5};
Physical Surface(3+MINUS_TWO_DIM) = {2,3,4,5};

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

Physical Volume(1) = {1};