FDS Source: Any zero or inflow non-HVAC BC has deposition.

Author: drjfloyd

Source/wall.f90

@@ -119,8 +119,8 @@ SUBROUTINE WALL_BC(T,DT,NM)
 
    IF (DEPOSITION .AND. .NOT.INITIALIZATION_PHASE .AND. CORRECTOR .AND. .NOT.SOLID_PHASE_ONLY) THEN
       IF (WC%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. &
-         (ANY(SF%LEAK_PATH>0) .OR. (B1%NODE_INDEX==0 .AND. ABS(SF%VEL)<TWO_EPSILON_EB .AND. &
-          ANY(ABS(SF%MASS_FLUX)<TWO_EPSILON_EB) .AND. ABS(SF%VOLUME_FLOW)<TWO_EPSILON_EB))) THEN
+         (ANY(SF%LEAK_PATH>0) .OR. &
+         (B1%NODE_INDEX==0 .AND. SF%VEL<TWO_EPSILON_EB .AND. SF%VOLUME_FLOW<TWO_EPSILON_EB))) THEN
          CALL CALC_DEPOSITION(DT,BC,B1,B2,WALL_INDEX=IW)
       ENDIF
    ENDIF
@@ -175,8 +175,8 @@ SUBROUTINE WALL_BC(T,DT,NM)
 
    IF (DEPOSITION .AND. .NOT.INITIALIZATION_PHASE .AND. CORRECTOR .AND. .NOT.SOLID_PHASE_ONLY) THEN
       IF (CFA%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. &
-         (ANY(SF%LEAK_PATH>0) .OR. (B1%NODE_INDEX==0 .AND. ABS(SF%VEL)<TWO_EPSILON_EB .AND. &
-          ANY(ABS(SF%MASS_FLUX)<TWO_EPSILON_EB) .AND. ABS(SF%VOLUME_FLOW)<TWO_EPSILON_EB))) THEN
+         (ANY(SF%LEAK_PATH>0) .OR. &
+         (B1%NODE_INDEX==0 .AND. SF%VEL<TWO_EPSILON_EB .AND. SF%VOLUME_FLOW<TWO_EPSILON_EB))) THEN
          CALL CALC_DEPOSITION(DT,BC,B1,B2,CFACE_INDEX=ICF)
       ENDIF
    ENDIF
@@ -1613,19 +1613,19 @@ SUBROUTINE CALC_DEPOSITION(DT,BC,B1,B2,WALL_INDEX,CFACE_INDEX)
 TYPE(BOUNDARY_COORD_TYPE), POINTER :: BC
 TYPE(BOUNDARY_PROP1_TYPE), POINTER :: B1
 TYPE(BOUNDARY_PROP2_TYPE), POINTER :: B2
-TYPE(SURFACE_TYPE), POINTER :: SF
+TYPE(SURFACE_TYPE), POINTER :: SFD
 
 IF (PRESENT(WALL_INDEX)) THEN
-   SF=>SURFACE(WALL(WALL_INDEX)%SURF_INDEX)
+   SFD=>SURFACE(WALL(WALL_INDEX)%SURF_INDEX)
    RSUM_G = RSUM(BC%IIG,BC%JJG,BC%KKG)
    ZZ_G(1:N_TRACKED_SPECIES) = ZZ(BC%IIG,BC%JJG,BC%KKG,1:N_TRACKED_SPECIES)
 ELSEIF (PRESENT(CFACE_INDEX)) THEN
-   SF=>SURFACE(CFACE(CFACE_INDEX)%SURF_INDEX)
+   SFD=>SURFACE(CFACE(CFACE_INDEX)%SURF_INDEX)
    RSUM_G= CFACE(CFACE_INDEX)%RSUM_G
    ZZ_G(1:N_TRACKED_SPECIES) = B1%ZZ_G(1:N_TRACKED_SPECIES)
 ENDIF
 
-IF (ANY(SF%LEAK_PATH>0)) THEN
+IF (ANY(SFD%LEAK_PATH>0)) THEN
    U_NORMAL = 0._EB
 ELSE
    U_NORMAL = B1%U_NORMAL
@@ -1673,7 +1673,6 @@ SUBROUTINE CALC_DEPOSITION(DT,BC,B1,B2,WALL_INDEX,CFACE_INDEX)
          U_TURB = B2%U_TAU * 0.17_EB
       ENDIF
    ENDIF
-
    B2%V_DEP = MAX(0._EB,U_THERM+U_TURB+U_GRAV+U_NORMAL)
    IF (B2%V_DEP <= TWO_EPSILON_EB) CYCLE SMIX_LOOP
    ZZ_GET = ZZ_GET * B1%RHO_G
@@ -1682,7 +1681,6 @@ SUBROUTINE CALC_DEPOSITION(DT,BC,B1,B2,WALL_INDEX,CFACE_INDEX)
    ZZ_GET(N) = Y_AEROSOL - YDEP
    IF (SM%AWM_INDEX > 0) B1%AWM_AEROSOL(SM%AWM_INDEX)= B1%AWM_AEROSOL(SM%AWM_INDEX)+YDEP/B1%RDN
    IF (SS%AWM_INDEX > 0) B1%AWM_AEROSOL(SS%AWM_INDEX)= B1%AWM_AEROSOL(SS%AWM_INDEX)+YDEP/B1%RDN
-
    !$OMP CRITICAL
    D_SOURCE(BC%IIG,BC%JJG,BC%KKG) = D_SOURCE(BC%IIG,BC%JJG,BC%KKG) - MW_RATIO*YDEP / B1%RHO_G / DT
    M_DOT_PPP(BC%IIG,BC%JJG,BC%KKG,N) = M_DOT_PPP(BC%IIG,BC%JJG,BC%KKG,N) - YDEP / DT