diff --git a/applications/solvers/electromagnetics/magneticFoam/magnet.H b/applications/solvers/electromagnetics/magneticFoam/magnet.H
index 9be4e6e92711ad2c6dd03faaf65db9f7c9ada7a2..765d513a74b634bf2ddcf2eb8fc9a89a0c126bcb 100644
--- a/applications/solvers/electromagnetics/magneticFoam/magnet.H
+++ b/applications/solvers/electromagnetics/magneticFoam/magnet.H
@@ -60,7 +60,7 @@ class magnet
     // Private data
 
         word name_;
-        scalar relativPermeability_;
+        scalar relativePermeability_;
         dimensionedScalar remanence_;
         vector orientation_;
 
@@ -85,7 +85,7 @@ public:
         )
         :
             name_(name),
-            relativPermeability_(mur),
+            relativePermeability_(mur),
             remanence_("Mr", dimensionSet(0, -1, 0, 0, 0, 1, 0), Mr),
             orientation_(orientation)
         {}
@@ -111,7 +111,7 @@ public:
         //- Return relative permeability
         inline scalar mur() const
         {
-            return relativPermeability_;
+            return relativePermeability_;
         }
 
         //- Return remenance
@@ -133,7 +133,7 @@ public:
         {
             is.readBegin("magnet");
             is  >> m.name_
-                >> m.relativPermeability_
+                >> m.relativePermeability_
                 >> m.remanence_.value()
                 >> m.orientation_;
             is.readEnd("magnet");
@@ -148,7 +148,7 @@ public:
         {
             os  << token::BEGIN_LIST
                 << m.name_ << token::SPACE
-                << m.relativPermeability_ << token::SPACE
+                << m.relativePermeability_ << token::SPACE
                 << m.remanence_.value()
                 << m.orientation_
                 << token::END_LIST;
diff --git a/applications/solvers/electromagnetics/magneticFoam/magneticFoam.C b/applications/solvers/electromagnetics/magneticFoam/magneticFoam.C
index d8b23f1b89ede9c996bcab7b282bae18768a90f9..e55979061dfacc806277744ac87c6d8b0c8f144e 100644
--- a/applications/solvers/electromagnetics/magneticFoam/magneticFoam.C
+++ b/applications/solvers/electromagnetics/magneticFoam/magneticFoam.C
@@ -80,6 +80,23 @@ int main(int argc, char *argv[])
         );
 
         H.write();
+
+        Info<< nl
+            << "Creating field HdotGradH for time "
+            << runTime.timeName() << endl;
+
+        volVectorField HdotGradH
+        (
+            IOobject
+            (
+                "HdotGradH",
+                runTime.timeName(),
+                mesh
+            ),
+            H & fvc::grad(H)
+        );
+
+        HdotGradH.write();
     }
 
     if (!args.optionFound("noB"))