stress that there is a single (fixed) dof for js

Magnetostatics/electromagnet.pro

@@ -218,10 +218,10 @@ Integration {
 
 /* The function space "Hregion_j_Mag_2D" provides one basis function, and hence
    one degree of freedom, per physical region in the abstract region
-   "Vol_S_Mag".  The constraint "SourceCurrentDensityZ" fixes all these dofs, so
-   the FunctionSpace "Hregion_j_Mag_2D" is fully fixed and has no FE unknowns.
-   One could thus have replaced it by a simple function and the Integral term
-   below in the weak formulation could have been written as
+   "Vol_S_Mag".  The constraint "SourceCurrentDensityZ" fixes all these dofs (a
+   single one!), so the FunctionSpace "Hregion_j_Mag_2D" is fully fixed and has
+   no FE unknowns.  One could thus have replaced it by a simple function and the
+   Integral term below in the weak formulation could have been written as
 
    Integral { [ Vector[ 0,0,-js_fct[] ] , {a} ];
      In Vol_S_Mag; Jacobian Vol; Integration Int; }