plot b and btarget on points of interest

Team25/shape.pro

@@ -313,6 +313,7 @@ PostProcessing {
       { Name w; Value{ Term{[(CompX[{d a}]-bx[])^2+(CompY[{d a}]-by[])^2 ]; In Domain; Jacobian Vol; } } }
       { Name az; Value{ Term{ [CompZ[{a}]]; In Domain; Jacobian Vol; } } }
       { Name b; Value{ Term{ [{d a}]; In Domain; Jacobian Vol; } } }
+      { Name bTarget; Value{ Term{ [ Vector[bx[],by[],0]  ]; In Domain; Jacobian Vol; } } }
       { Name bMag; Value{ Term{ [Norm[{d a}]]; In Domain; Jacobian Vol; } } }
       { Name js; Value{ Term{ [js[]]; In Coil; Jacobian Vol; } } }
       { Name mur; Value{ Term{ [1/(nu[{d a}]*mu0)]; In Domain; Jacobian Vol; } } }
@@ -330,8 +331,12 @@ PostOperation Get_ObjectiveConstraints UsingPost ObjectiveConstraints {
     Format SimpleTable, File StrCat[Opt_ResDir,"bMag.txt"]];
   Print[az, OnElementsOf Domain, Format NodeTable, File "",
     SendToServer StrCat[Opt_ResDir_Onelab,"a"], Hidden 1];
-  Print[bMag, OnElementsOf Domain, File StrCat[Opt_ResDir,"az.pos"]];
+  //Print[bMag, OnElementsOf Domain, File StrCat[Opt_ResDir,"az.pos"]];
   //Print[az, OnElementsOf Domain, File StrCat[Opt_ResDir,"az.pos"]];
+  Print[bTarget, OnGrid {(9.5e-3+2.25e-3)*Cos[$A],(9.5e-3+2.25e-3)*Sin[$A],0}{0:50*Pi/180:5*Pi/180,0,0},
+    File StrCat[Opt_ResDir,"bTarget.pos"]];
+  Print[b, OnGrid {(9.5e-3+2.25e-3)*Cos[$A],(9.5e-3+2.25e-3)*Sin[$A],0}{0:50*Pi/180:5*Pi/180,0,0},
+    File StrCat[Opt_ResDir,"b.pos"]];
   If(Flag_PrintLevel>5)
     Print[mur, OnElementsOf Domain, File StrCat[Opt_ResDir,"mur.pos"]];
     Print[az, OnElementsOf Domain, File StrCat[Opt_ResDir,"az.pos"]];