diff --git a/Magnetodynamics/transfo.pro b/Magnetodynamics/transfo.pro
index d9df73dc26d05860dce0379be3899e51729a26ab..3d7cd43317f0dbb42de4c1f6a2719f5218a1da00 100644
--- a/Magnetodynamics/transfo.pro
+++ b/Magnetodynamics/transfo.pro
@@ -122,6 +122,8 @@ ElseIf (type_Source == 2) // voltage
     // Primary side
     E_in = Region[10001]; // arbitrary region number (not linked to the mesh)
     SourceV_Cir += Region[{E_in}];
+    R_in = Region[10002]; // arbitrary region number (not linked to the mesh)
+    Resistance_Cir += Region[{R_in}];
 
     // Secondary side
     R_out = Region[10101]; // arbitrary region number (not linked to the mesh)
@@ -134,9 +136,14 @@ ElseIf (type_Source == 2) // voltage
     // are defined)
     val_E_in = 1.;
     phase_E_in = 90 *deg; // Phase in radian (from phase in degree)
+
     // High value for an open-circuit test; Low value for a short-circuit test;
     // any value in-between for any charge
     Resistance[R_out] = 1e6;
+
+    // End-winding primary winding resistance for more realistic primary coil
+    // model
+    Resistance[R_in] = 1;
   }
 
   Constraint {
@@ -157,7 +164,8 @@ ElseIf (type_Source == 2) // voltage
         // PLUS and MINUS coil portions to be connected in series, together with
         // E_in (an additional resistor should be defined to represent the
         // Coil_1 end-winding (not considered in the 2D model))
-        { Region E_in; Branch {1,2}; }
+        { Region E_in; Branch {1,4}; }
+        { Region R_in; Branch {4,2}; }
 
         { Region Coil_1_P; Branch {2,3} ; }
         { Region Coil_1_M; Branch {3,1} ; }