fixing geo team32 -> Abs(CombinedBoundary{***})

Team32/t32.geo

@@ -271,17 +271,17 @@ If(Flag_3Dmodel==0)
   //=================================================
   // Physical regions for FE analysis (2D)
   //=================================================
-  Physical Surface(CORE)   = surf_ECore[];
+    Physical Surface("Core", CORE)   = surf_ECore[];
 
   nb_surf_coil = #surf_Coil[];
-  Physical Surface(COIL)   = surf_Coil[{0:nb_surf_coil-1:4}];
-  Physical Surface(COIL+1) = surf_Coil[{1:nb_surf_coil-1:4}];
-  Physical Surface(COIL+2) = surf_Coil[{2:nb_surf_coil-1:4}];
-  Physical Surface(COIL+3) = surf_Coil[{3:nb_surf_coil-1:4}];
+  Physical Surface("Coil right (left side)", COIL)   = surf_Coil[{0:nb_surf_coil-1:4}];
+  Physical Surface("Coil right (right side)", COIL+1) = surf_Coil[{1:nb_surf_coil-1:4}];
+  Physical Surface("Coil left (right side)", COIL+2) = surf_Coil[{2:nb_surf_coil-1:4}];
+  Physical Surface("Coil left (left side)", COIL+3) = surf_Coil[{3:nb_surf_coil-1:4}];
 
-  Physical Surface(AIR)    = surf_Air[];
-  Physical Surface(AIRINF) = surf_AirInf[];
-  Physical Line(SURF_AIROUT) = ln_bnd[];
+  Physical Surface("Air", AIR)    = surf_Air[];
+  Physical Surface("Air infinity shell", AIRINF) = surf_AirInf[];
+  Physical Line("Outer boundary", SURF_AIROUT) = ln_bnd[];
 
     //=================================================
   // Some colors... for aesthetics :-)
@@ -305,44 +305,44 @@ If(Flag_3Dmodel==1)
   //=================================================
   // Physical regions for FE analysis (3D)
   //=================================================
-  bnd_Core[]  = CombinedBoundary{Volume{vol_Core[]};};
+  bnd_Core[]  = Abs(CombinedBoundary{Volume{vol_Core[]};});
   bnd_Core[] -= surf_cut_xy[];
   If (Flag_Symmetry3D ==1)
     surf_cut_xz[] += bnd_Core[{5,9,16,20}];
     bnd_Core[] -= surf_cut_xz[];
   EndIf
   //Printf("",bnd_Core[]);
-  Physical Volume(CORE)      = vol_Core[];
-  Physical Surface(SKINCORE) = bnd_Core[];
+  Physical Volume("Core", CORE)      = vol_Core[];
+  Physical Surface("Skin Core", SKINCORE) = bnd_Core[];
 
-  Physical Volume(COIL+0)       = vol_CoilR[];
-  Physical Volume(COIL+1)       = vol_CoilL[];
+  Physical Volume("Coil right", COIL+0)       = vol_CoilR[];
+  Physical Volume("Coil left", COIL+1)       = vol_CoilL[];
 
-  Physical Volume(LEG_INCOIL+0) = vol_inCoilR[]; // for source computation
-  Physical Volume(LEG_INCOIL+1) = vol_inCoilL[];
+  Physical Volume("Leg in Coil right",LEG_INCOIL+0) = vol_inCoilR[]; // for source computation
+  Physical Volume("Leg in Coil left", LEG_INCOIL+1) = vol_inCoilL[];
 
   nb_surf_coil = #surf_Coil[];
-  Physical Surface(ELECCOIL+0) = surf_Coil[{1:nb_surf_coil-1:4}]; // e.g. the one of the right
-  Physical Surface(ELECCOIL+1) = surf_Coil[{2:nb_surf_coil-1:4}];
+  Physical Surface("Electrode Coil right",ELECCOIL+0) = surf_Coil[{1:nb_surf_coil-1:4}]; // e.g. the one of the right
+  Physical Surface("Electrode Coil left", ELECCOIL+1) = surf_Coil[{2:nb_surf_coil-1:4}];
 
   //Right
   // ======================================================
-  bnd_CoilR[]  = CombinedBoundary{Volume{vol_CoilR[]};};
+  bnd_CoilR[]  = Abs(CombinedBoundary{Volume{vol_CoilR[]};});
   bnd_CoilR[] -= surf_Coil[{0:nb_surf_coil-1:4,1:nb_surf_coil-1:4}];
   If (Flag_Symmetry3D ==1)
     surf_cut_xz[] += bnd_CoilR[{4,6,9,12,16}];
     bnd_CoilR[]   -= surf_cut_xz[];
   EndIf
-  Physical Surface(SKINCOIL+0) = bnd_CoilR[];
+  Physical Surface("Skin Coil right", SKINCOIL+0) = bnd_CoilR[];
 
-  bnd_CoilR_plus_hole[] = CombinedBoundary{Volume{vol_CoilR[], vol_inCoilR[]};};
-  bnd_inCoilR[]  = CombinedBoundary{Volume{vol_inCoilR[]};};
+  bnd_CoilR_plus_hole[] = Abs(CombinedBoundary{Volume{vol_CoilR[], vol_inCoilR[]};});
+    bnd_inCoilR[]  = Abs(CombinedBoundary{Volume{vol_inCoilR[]};});
   bnd_inCoilR[] -= bnd_CoilR_plus_hole[];
-  Physical Surface(SKINCOIL_CUT+0) = bnd_inCoilR[];
+Physical Surface("Skin Coil right- only around cut", SKINCOIL_CUT+0) = bnd_inCoilR[];
 
-  bnd_CoilR_plus_hole[] -= CombinedBoundary{Volume{vol_CoilR[]};};
+  bnd_CoilR_plus_hole[] -= Abs(CombinedBoundary{Volume{vol_CoilR[]};});
   bnd_CoilR_plus_hole[] -= surf_cut_xy[];
-  Physical Surface(COIL_CUT+0)= bnd_CoilR_plus_hole[1]; //cut down (top at 0)
+  Physical Surface("Cut Coil right", COIL_CUT+0)= bnd_CoilR_plus_hole[1]; //cut down (top at 0)
 
   If (Flag_Symmetry3D ==1)
     surf_cut_xz[] += bnd_CoilR_plus_hole[1];
@@ -350,37 +350,37 @@ If(Flag_3Dmodel==1)
 
   // Left
   // ======================================================
-  bnd_CoilL[]  = CombinedBoundary{Volume{vol_CoilL[]};};
+  bnd_CoilL[]  = Abs(CombinedBoundary{Volume{vol_CoilL[]};});
   bnd_CoilL[] -= surf_Coil[{2:nb_surf_coil-1:4,3:nb_surf_coil-1:4}];
   If (Flag_Symmetry3D ==1)
      surf_cut_xz[] += bnd_CoilL[{4,6,9,12,16}];
      bnd_CoilL[]   -= surf_cut_xz[];
    EndIf
-  Physical Surface(SKINCOIL+1) = bnd_CoilL[];
+   Physical Surface("Skin Coil left", SKINCOIL+1) = bnd_CoilL[];
 
 
-  bnd_CoilL_plus_hole[] = CombinedBoundary{Volume{vol_CoilL[], vol_inCoilL[]};};
-  bnd_inCoilL[]  = CombinedBoundary{Volume{vol_inCoilL[]};};
+  bnd_CoilL_plus_hole[] = Abs(CombinedBoundary{Volume{vol_CoilL[], vol_inCoilL[]};});
+    bnd_inCoilL[]  = Abs(CombinedBoundary{Volume{vol_inCoilL[]};});
   bnd_inCoilL[] -= bnd_CoilL_plus_hole[];
-  Physical Surface(SKINCOIL_CUT+1) = bnd_inCoilL[];
+  Physical Surface("Skin Coil left- only around cut", SKINCOIL_CUT+1) = bnd_inCoilL[];
 
-  bnd_CoilL_plus_hole[] -= CombinedBoundary{Volume{vol_CoilL[]};};
+  bnd_CoilL_plus_hole[] -= Abs(CombinedBoundary{Volume{vol_CoilL[]};});
   bnd_CoilL_plus_hole[] -= surf_cut_xy[];
-  Physical Surface(COIL_CUT+1)= bnd_CoilL_plus_hole[1]; //cut down (top at 0)
+  Physical Surface("Cut Coil left", COIL_CUT+1)= bnd_CoilL_plus_hole[1]; //cut down (top at 0)
 
   If (Flag_Symmetry3D ==1)
     surf_cut_xz[] += bnd_CoilL_plus_hole[1];
   EndIf
 
-  Physical Volume(AIR)         = vol_Air[];
-  Physical Surface(SURF_CUTXY) = surf_cut_xy[];
-  Physical Surface(SURF_CUTXZ) = surf_cut_xz[];
+  Physical Volume("Air", AIR)         = vol_Air[];
+  Physical Surface("Symmetry cut XY", SURF_CUTXY) = surf_cut_xy[];
+  Physical Surface("Symmetry cut XZ", SURF_CUTXZ) = surf_cut_xz[];
 
   all_vol[] = Volume '*';
-  bnd_all[] = CombinedBoundary{Volume{all_vol[]};};
+  bnd_all[] = Abs(CombinedBoundary{Volume{all_vol[]};});
   bnd_all[] -= surf_cut_xy[];
   bnd_all[] -= surf_cut_xz[];
-  Physical Surface(SURF_AIROUT) = bnd_all[];
+  Physical Surface("Outer boundary", SURF_AIROUT) = bnd_all[];
 
   //=================================================
   // Some colors... for aesthetics :-)
@@ -405,7 +405,3 @@ If(PostProcessing.NbViews==0)
     };
   EndFor
 EndIf
-
-
-
-