2D mesh error: unable to recover the edge
Hello,
I want to generate a 2D mesh of simple planar surfaces geometry. For further handling of the mesh I need to insert 2 physical points so that I have a clear allocation between the triangles and the points. However, the mesh generation fails whenever the part for the physical points is added (last 6 lines, see below or in attached file GND_Antenna.geo). Without these lines the mesh can be generated. I do not see what is causing the error message "Unable to recover the edge...". I would very happy to get a hint. Here are the lines to create the geometry:
Content of geo-file GND_Antenna.geo
`SetFactory("OpenCASCADE");cl1=0.002999;
p1=newp; Point(p1) = {-0.0025,0.0005,0,cl1};
p2=newp; Point(p2) = {-0.0025,0.0025,0,cl1};
p3=newp; Point(p3) = {-0.025,0.0025,0,cl1};
p4=newp; Point(p4) = {-0.025,0.05,0,cl1};
p5=newp; Point(p5) = {0.025,0.05,0,cl1};
p6=newp; Point(p6) = {0.025,0.0025,0,cl1};
p7=newp; Point(p7) = {0.0025,0.0025,0,cl1};
p8=newp; Point(p8) = {0.0025,0.0005,0,cl1};
l1=newl; Line(l1) = {p1,p2};
l2=newl; Line(l2) = {p2,p3};
l3=newl; Line(l3) = {p3,p4};
l4=newl; Line(l4) = {p4,p5};
l5=newl; Line(l5) = {p5,p6};
l6=newl; Line(l6) = {p6,p7};
l7=newl; Line(l7) = {p7,p8};
l8=newl; Line(l8) = {p8,p1};
ll1 = newll ;
Curve Loop(ll1) = {l1,l2,l3,l4,l5,l6,l7,l8};
s1 = news;
Plane Surface(s1) = {ll1} ;
p1=newp; Point(p1) = {-0.0025,-0.0005,0,cl1};
p2=newp; Point(p2) = {-0.0025,-0.0025,0,cl1};
p3=newp; Point(p3) = {-0.025,-0.0025,0,cl1};
p4=newp; Point(p4) = {-0.025,-0.05,0,cl1};
p5=newp; Point(p5) = {0.025,-0.05,0,cl1};
p6=newp; Point(p6) = {0.025,-0.0025,0,cl1};
p7=newp; Point(p7) = {0.0025,-0.0025,0,cl1};
p8=newp; Point(p8) = {0.0025,-0.0005,0,cl1};
l1=newl; Line(l1) = {p1,p2};
l2=newl; Line(l2) = {p2,p3};
l3=newl; Line(l3) = {p3,p4};
l4=newl; Line(l4) = {p4,p5};
l5=newl; Line(l5) = {p5,p6};
l6=newl; Line(l6) = {p6,p7};
l7=newl; Line(l7) = {p7,p8};
l8=newl; Line(l8) = {p8,p1};
ll1 = newll ;
Curve Loop(ll1) = {l1,l2,l3,l4,l5,l6,l7,l8};
s2 = news;
Plane Surface(s2) = {ll1} ;
Physical Surface (1) = {s1};
Physical Surface (2) = {s2};
p1 = newp; Point(p1) = {0,-0.0005,0};
Point{p1} In Surface{s2};
Physical Point(1)={p1};
p2 = newp; Point(p2) = {0,0.0005,0};
Point{p2} In Surface{s1};
Physical Point(2)={p2};`
I thank you in advance for your support!