lc = 0.009;
Point(1) = {0, 0, 0, lc};
Point(2) = {0.1, 0, 0, lc};
Point(3) = {0.1, 0.3, 0, lc};
Point(4) = {0, 0.3, 0, lc};
Line(1) = {1, 2};
Line(2) = {3, 2};
Line(3) = {3, 4};
Line(4) = {4, 1};
Line Loop(5) = {4,1,-2,3};
Plane Surface(6) = {5};

lc_c = 0.009; cx=0.05; cy=0.15; cz=-0.2;

Point(9)  = {cx, cy, 1*cz/5, lc_c};
Point(10) = {cx, cy, 0, lc_c};
Line(14) = {9,10};

Point{10} In Surface{6};


Field[1] = Attractor;
Field[1].NNodesByEdge = 50;
Field[1].EdgesList = {1, 2, 3, 4};

Field[2] = Threshold;
Field[2].IField = 1;
Field[2].Sigmoid = 1;
Field[2].LcMax = lc;
Field[2].DistMin = lc/3;
Field[2].LcMin = lc/3;
Field[2].DistMax = lc;

Field[3] = Threshold;
Field[3].IField = 1;
Field[3].Sigmoid = 1;
Field[3].LcMax = lc;
Field[3].DistMin = lc;
Field[3].LcMin = lc/2;
Field[3].DistMax = 1.5*lc;

Field[4] = Min;
Field[4].FieldsList = {2,3};
Background Field = 4;

// Don't extend the elements sizes from the boundary inside the domain
Mesh.CharacteristicLengthExtendFromBoundary = 0;

Mesh.ColorCarousel = 0;

