Merge pull request #12240 from mcgratta/master

FDS Source: Issue #12239. Correct sign of diffusion term in TOTAL MAS…
firemodels · Nov 17, 2023 · 035852d · 035852d
2 parents afc1ca6 + 08b8a01
commit 035852d
Showing 1 changed file with 6 additions and 6 deletions.
diff --git 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