diff --git a/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro b/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
index 716a0453e500539014fdcb0739779bada164eb73..846d4103929ad4d84d241b6a3d4ae1b226b23b38 100644
--- a/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
+++ b/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
@@ -556,6 +556,10 @@ PostProcessing {
           Term { [ {d a} ]; In Vol_Mag; Jacobian Vol; }
         }
       }
+      { Name norm_of_b; Value {
+          Term { [ Norm[{d a}] ]; In Vol_Mag; Jacobian Vol; }
+        }
+      }
       { Name h; Value {
           Term { [ nu[] * {d a} ]; In Vol_Mag; Jacobian Vol; }
         }
diff --git a/Magnetodynamics/transfo.pro b/Magnetodynamics/transfo.pro
index 3ec7ac76be17beffb83d66a79e37b168c45207f4..d8d3ebb906b1fe92eab6d100e541040ec3170320 100644
--- a/Magnetodynamics/transfo.pro
+++ b/Magnetodynamics/transfo.pro
@@ -21,7 +21,10 @@ DefineConstant[
   // way as far as circuit-coupling is concerned
   ConductorType = {2, Choices{1 = "Massive", 2 = "Coil"}, Highlight "Blue",
     Name "Parameters/01Conductor type"}
-  Freq = {1, Min 0, Max 1e3, Step 1, Name "Parameters/Frequency"}
+  Freq = {1, Min 0, Max 1e3, Step 1,
+    Name "Parameters/Frequency"}
+  Flag_FrequencyDomain = {0, Choices{0, 1},
+    Name "Parameters/Frequency-domain?"}
   mur_Core = {1000, Min 1, Max 10000, Step 1,
     Name "Parameters/Core relative permeability"}
 ];
@@ -94,8 +97,6 @@ Function {
 // secondary.
 Flag_CircuitCoupling = 1;
 
-Flag_FrequencyDomain = 0;
-
 // Note that the voltage will not be equally distributed in the PLUS and MINUS
 // parts, which is the reason why we must apply the total voltage through a
 // circuit -- and we cannot simply use a current source like in Tutorial 7a.