modified model to have b along e_y, nicer pic

Magnetostatic/electromagnet.pro

@@ -88,7 +88,9 @@ Function {
   NbTurns = 1000 ;
   Current = DefineNumber[0.01, Name "Model parameters/Current",
 			 Help "Current injected in coil [A]"]; 
-  Js_fct[ Ind ] = NbTurns*Current/SurfaceArea[];
+  Js_fct[ Ind ] = -NbTurns*Current/SurfaceArea[];
+  /* The minus sign is to have the current in -e_z direction, 
+     so that the magnetic induction field is in +e_y direction */  
 }
 
 /* In the 2D approximation, the magnetic vector potential A and the current density Js
@@ -238,17 +240,20 @@ PostProcessing {
 }
 
 e = 1.e-5;
-p1 = {e,e,0};
-p2 = {0.25-e,e,0}; // horizontal cut through model, just above x-axis. 
+h = 0.02;
+p1 = {e,h,0};
+p2 = {0.25-e,h,0}; // horizontal cut through model, just above x-axis. 
 
 PostOperation {
 
   { Name Map_a; NameOfPostProcessing MagSta_a_2D;
     Operation {
-      Echo[ Str["l=PostProcessing.NbViews-1; View[l].IntervalsType = 3;"],
+      Echo[ Str["l=PostProcessing.NbViews-1;",
+		"View[l].IntervalsType = 1;",
+		"View[l].NbIso = 40;"],
 	    File "tmp.geo", LastTimeStepOnly] ;
       Print[ az, OnElementsOf Dom_Hcurl_a_Mag_2D, File "az.pos" ];
-      Print[ b, OnLine{{List[p1]}{List[p2]}} {1000}, File "by.pos" ];
+      Print[ b, OnLine{{List[p1]}{List[p2]}} {50}, File "by.pos" ];
     }
   }
 }