diff --git a/Antennas/dipole3d.geo b/Antennas/dipole3d.geo
index c2f94f227ece5056aae16926d8717bc260da90da..d75a7e7c63998af95f7484501a5c28b555510fd7 100644
--- a/Antennas/dipole3d.geo
+++ b/Antennas/dipole3d.geo
@@ -1,7 +1,7 @@
-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;
diff --git a/Antennas/dipole_data.geo b/Antennas/dipole_data.geo
index 78c6a30d5dd0695787be2fedfb8acb2eb6e9244e..3b8a1ef9e5f09ecd875f3afe0b73ff6ec62ea1fe 100644
--- a/Antennas/dipole_data.geo
+++ b/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);