refine close to feed

Antennas/dipole3d.geo

-fac = 2 ;
+fac = 5 ;
 
-Mesh.CharacteristicLengthFactor = fac ;
-Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
+Mesh.MeshSizeFactor = fac ;
+Mesh.Algorithm3D = 10; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
 Mesh.Optimize = 1 ;
 
 // characteristic lengths & some transfinite number of divisions
@@ -243,3 +243,14 @@ If(CoefGeo==1 && !Flag_PML_Cyl)
   Physical Surface(SKINDIPOLEUP)  = skindipole1[];
   Physical Surface(SKINFEED)   = skindipole2[];
 EndIf
+
+
+Field[1] = Ball;
+Field[1].VIn = lcd/5;
+Field[1].VOut = lc;
+Field[1].XCenter = 0;
+Field[1].YCenter = 0;
+Field[1].ZCenter = 0;
+Field[1].Radius = 1*delta_gap;
+Field[1].Thickness = 4*delta_gap;
+Background Field = 1;

Antennas/dipole_data.geo

@@ -71,7 +71,7 @@ DefineConstant[
 rdipole = rdipole*mm; // in [m]
 AngleWedge = ((Flag_PML_Cyl==0) ? AngleWedge_deg : AngleWedgeCyl_deg) * deg2rad ;
 
-CoefGeo = (!Flag_3Dmodel) ? 2*Pi : 2*Pi/AngleWedge; // axisymmetry in 2D, 1/8 or 1/4 of the 3D model
+CoefGeo = (!Flag_3Dmodel) ? 2*Pi : (2*Pi/AngleWedge); // axisymmetry in 2D, 1/8 or 1/4 of the 3D model
 
 Printf("CoefGeo %g", CoefGeo);