diff --git a/src/mhd/DirectInduction/mhd.h b/src/mhd/DirectInduction/mhd.h
index ef83e645851d100fe2c43ea108d5236a8bf43c7b..558eda981dc71cdbf1ac16eced4ed25519175cd4 100644
--- a/src/mhd/DirectInduction/mhd.h
+++ b/src/mhd/DirectInduction/mhd.h
@@ -348,7 +348,7 @@ __attribute__((always_inline)) INLINE static void mhd_end_gradient(
   B[2] = p->mhd_data.B_over_rho[2] * rho;
 
   float OW;
-  OW = 1.0f;
+  OW = 2.0f;
 
   float absB;
   absB = sqrtf(B[0]*B[0]+B[1]*B[1]+B[2]*B[2]);
diff --git a/src/mhd/DirectInduction/mhd_iact.h b/src/mhd/DirectInduction/mhd_iact.h
index bcc766a34b074a45f805f7a0fb48d06fe3680e4f..c22248179c1c2fc9455ba3922abd8216afdac687 100644
--- a/src/mhd/DirectInduction/mhd_iact.h
+++ b/src/mhd/DirectInduction/mhd_iact.h
@@ -160,7 +160,7 @@ __attribute__((always_inline)) INLINE static void runner_iact_mhd_gradient(
   /* Calculate OWAR */
 
   float OW;
-  OW = 1.0f;
+  OW = 2.0f;
 
   float absBi;
   absBi = sqrtf(Bi[0]*Bi[0]+Bi[1]*Bi[1]+Bi[2]*Bi[2]);
@@ -309,7 +309,7 @@ runner_iact_nonsym_mhd_gradient(const float r2, const float dx[3],
   /* Calculate OWAR */
 
   float OW;
-  OW = 1.0f;
+  OW = 2.0f;
 
   float absBj;
   absBj = sqrtf(Bj[0]*Bj[0]+Bj[1]*Bj[1]+Bj[2]*Bj[2]);