diff --git a/sources/sources.F90 b/sources/sources.F90
index 201a456..be5da5e 100644
--- a/sources/sources.F90
+++ b/sources/sources.F90
@@ -571,21 +571,16 @@ module sources
 !
         call divergence(dh(:), pdata%q(ibx:ibz,:,:,:), db(:,:,:))
 
-! update the momentum component increments (only in the adiabatic case for
-! KEPES), i.e.
+! update the momentum component increments, i.e.
 !     d/dt (ρv) + ∇.F = - (∇.B)B
 !
-        if (glm_type < 3 .or. ien > 0) then
-
-          pdata%du(imx,:,:,:) = pdata%du(imx,:,:,:) &
+        pdata%du(imx,:,:,:) = pdata%du(imx,:,:,:) &
                                                - db(:,:,:) * pdata%q(ibx,:,:,:)
-          pdata%du(imy,:,:,:) = pdata%du(imy,:,:,:) &
+        pdata%du(imy,:,:,:) = pdata%du(imy,:,:,:) &
                                                - db(:,:,:) * pdata%q(iby,:,:,:)
-          pdata%du(imz,:,:,:) = pdata%du(imz,:,:,:) &
+        pdata%du(imz,:,:,:) = pdata%du(imz,:,:,:) &
                                                - db(:,:,:) * pdata%q(ibz,:,:,:)
 
-        end if ! ien > 0
-
 ! add the HEGLM-MHD and KEPES source terms
 !
         if (glm_type > 1) then