diff --git a/Magnetodynamics/electromagnet.pro b/Magnetodynamics/electromagnet.pro
index b5e7138e09ad332ffd43db3f16ba3594bfc02bfd..a381e598b2fa389d36ee9badc349d12dc647c848 100644
--- a/Magnetodynamics/electromagnet.pro
+++ b/Magnetodynamics/electromagnet.pro
@@ -14,6 +14,8 @@
 
    To compute the solution interactively from the Gmsh GUI:
        File > Open > electromagnet.pro
+       You may choose the Resolution in the left panel: 
+       Magnetodynamics2D_av (default) or Magnetostatics2D_a
        Run (button at the bottom of the left panel)
    ------------------------------------------------------------------- */
 
@@ -32,8 +34,8 @@ Group {
   Vol_Mag = Region[{Air, Core, Ind, AirInf}]; // full magnetic domain
   Vol_C_Mag = Region[Core]; // massive conductors
   Vol_S_Mag = Region[Ind]; // stranded conductors (coils)
-  Vol_Inf_Mag = Region[AirInf]; // annulus for infinite shell transformation
-  Val_Rint = rInt; Val_Rext = rExt; // interior and exterior radii of annulus
+  Vol_Inf_Mag = Region[AirInf]; // ring-shaped shell for infinite transformation
+  Val_Rint = rInt; Val_Rext = rExt; // interior and exterior radii of ring
 }
 
 Function {
@@ -50,7 +52,7 @@ Function {
   sigma[ Ind ] = 5e7;
 
   Ns[ Ind ] = 1000 ; // number of turns in coil
-  Sc[ Ind ] = SurfaceArea[] ; // surface (cross section) of coil
+  Sc[ Ind ] = SurfaceArea[] ; // area of coil cross section
   // Current density in each coil portion for a unit current (will be multiplied
   // by the actual total current in the coil)
   js0[ Ind ] = Ns[]/Sc[] * Vector[0,0,-1];
@@ -66,7 +68,7 @@ Constraint {
   }
   { Name Current_2D;
     Case {
-      // represents the phasor amplitude for a dynamic analysis
+      // represents the phasor amplitude (peak to peak value) for a dynamic analysis
       { Region Ind; Value Current; }
     }
   }