SetFactory("OpenCASCADE");

lref = 1.0;
np=newp;
line=newl;

Point(np+ 1)={   -1.5000,    0.0000,    0.0000,lref};
Point(np+ 2)={   -1.4808,    0.1951,    0.0000,lref};
Point(np+ 3)={   -1.4239,    0.3827,    0.0000,lref};
Point(np+ 4)={   -1.3315,    0.5556,    0.0000,lref};
Point(np+ 5)={   -1.2071,    0.7071,    0.0000,lref};
Point(np+ 6)={   -0.8725,    1.0097,    0.0000,lref};
Point(np+ 7)={   -0.5379,    1.2509,    0.0000,lref};
Point(np+ 8)={   -0.2033,    1.4350,    0.0000,lref};
Point(np+ 9)={    0.1313,    1.5662,    0.0000,lref};
Point(np+10)={    0.4659,    1.6487,    0.0000,lref};
Point(np+11)={    0.8005,    1.6867,    0.0000,lref};
Point(np+12)={    1.1351,    1.6843,    0.0000,lref};
Point(np+13)={    1.4697,    1.6458,    0.0000,lref};
Point(np+14)={    1.8043,    1.5754,    0.0000,lref};
Point(np+15)={    2.1389,    1.4772,    0.0000,lref};
Point(np+16)={    2.3507,    1.3482,    0.0000,lref};
Point(np+17)={    2.5353,    1.1147,    0.0000,lref};
Point(np+18)={    2.6784,    0.7949,    0.0000,lref};
Point(np+19)={    2.7690,    0.4135,    0.0000,lref};
Point(np+20)={    2.8000,    0.0000,    0.0000,lref};

Spline(line+1) = {
                   np+  1,np+  2,np+  3,np+  4,np+  5,np+  6,
                   np+  7,np+  8,np+  9,np+ 10,np+ 11,np+ 12,
                   np+ 13,np+ 14,np+ 15,np+ 16,np+ 17,np+ 18,
                   np+ 19,np+ 20};
Line(line+2) ={np+1,np+20};
Line Loop(1) = {line+1,line+2};
Plane Surface(1) = 1;
Extrude {{1, 0, 0}, {0, 0, 0}, 2*Pi} { Surface{1}; }
Delete {Surface{1};Line{3};}
