diff --git a/Source/divg.f90 b/Source/divg.f90
index 4948ce53d0e..1a84067dd0f 100644
--- a/Source/divg.f90
+++ b/Source/divg.f90
@@ -349,20 +349,20 @@ SUBROUTINE DIVERGENCE_PART_1(T,DT,NM)
          IF (STORE_SPECIES_FLUX) THEN
             IF (CORRECTOR) THEN
                SELECT CASE(IOR)
-                  CASE(-1) ; DIF_FXS(IIG  ,JJG,KKG,N) = -RHO_D_DZDN
+                  CASE(-1) ; DIF_FXS(IIG  ,JJG,KKG,N) =  RHO_D_DZDN
                   CASE( 1) ; DIF_FXS(IIG-1,JJG,KKG,N) = -RHO_D_DZDN
-                  CASE(-2) ; DIF_FYS(IIG,JJG  ,KKG,N) = -RHO_D_DZDN
+                  CASE(-2) ; DIF_FYS(IIG,JJG  ,KKG,N) =  RHO_D_DZDN
                   CASE( 2) ; DIF_FYS(IIG,JJG-1,KKG,N) = -RHO_D_DZDN
-                  CASE(-3) ; DIF_FZS(IIG,JJG,KKG  ,N) = -RHO_D_DZDN
+                  CASE(-3) ; DIF_FZS(IIG,JJG,KKG  ,N) =  RHO_D_DZDN
                   CASE( 3) ; DIF_FZS(IIG,JJG,KKG-1,N) = -RHO_D_DZDN
                END SELECT
             ELSE
                SELECT CASE(IOR)
-                  CASE(-1) ; DIF_FX(IIG  ,JJG,KKG,N) = 0.5_EB*(DIF_FXS(IIG , JJG,KKG,N)-RHO_D_DZDN)
+                  CASE(-1) ; DIF_FX(IIG  ,JJG,KKG,N) = 0.5_EB*(DIF_FXS(IIG , JJG,KKG,N)+RHO_D_DZDN)
                   CASE( 1) ; DIF_FX(IIG-1,JJG,KKG,N) = 0.5_EB*(DIF_FXS(IIG-1,JJG,KKG,N)-RHO_D_DZDN)
-                  CASE(-2) ; DIF_FY(IIG,JJG  ,KKG,N) = 0.5_EB*(DIF_FYS(IIG,JJG  ,KKG,N)-RHO_D_DZDN)
+                  CASE(-2) ; DIF_FY(IIG,JJG  ,KKG,N) = 0.5_EB*(DIF_FYS(IIG,JJG  ,KKG,N)+RHO_D_DZDN)
                   CASE( 2) ; DIF_FY(IIG,JJG-1,KKG,N) = 0.5_EB*(DIF_FYS(IIG,JJG-1,KKG,N)-RHO_D_DZDN)
-                  CASE(-3) ; DIF_FZ(IIG,JJG,KKG  ,N) = 0.5_EB*(DIF_FZS(IIG,JJG,KKG  ,N)-RHO_D_DZDN)
+                  CASE(-3) ; DIF_FZ(IIG,JJG,KKG  ,N) = 0.5_EB*(DIF_FZS(IIG,JJG,KKG  ,N)+RHO_D_DZDN)
                   CASE( 3) ; DIF_FZ(IIG,JJG,KKG-1,N) = 0.5_EB*(DIF_FZS(IIG,JJG,KKG-1,N)-RHO_D_DZDN)
                END SELECT
             ENDIF