FDS Source: set DRAG_MAX consistently for all stationary particles

Author: rmcdermo

Source/part.f90

@@ -2955,9 +2955,6 @@ SUBROUTINE MOVE_IN_GAS
             ACCEL_X = -(K_TERM(1)+Y_TERM(1))*UBAR*DY(JJG_OLD)*DZ(KKG_OLD)*SFAC
             ACCEL_Y = -(K_TERM(2)+Y_TERM(2))*VBAR*DX(IIG_OLD)*DZ(KKG_OLD)*SFAC
             ACCEL_Z = -(K_TERM(3)+Y_TERM(3))*WBAR*DX(IIG_OLD)*DY(JJG_OLD)*SFAC
-            DRAG_MAX(1) = ACCEL_X/(-UBAR+TWO_EPSILON_EB)
-            DRAG_MAX(2) = ACCEL_Y/(-VBAR+TWO_EPSILON_EB)
-            DRAG_MAX(3) = ACCEL_Z/(-WBAR+TWO_EPSILON_EB)
          ELSE
             ACCEL_X = 0._EB
             ACCEL_Y = 0._EB
@@ -2973,10 +2970,10 @@ SUBROUTINE MOVE_IN_GAS
          ACCEL_X = -(MIN(DX(IIG_OLD),LP%DX)/DX(IIG_OLD))*(K_TERM(1)+Y_TERM(1))*UBAR*SFAC
          ACCEL_Y = -(MIN(DY(JJG_OLD),LP%DY)/DY(JJG_OLD))*(K_TERM(2)+Y_TERM(2))*VBAR*SFAC
          ACCEL_Z = -(MIN(DZ(KKG_OLD),LP%DZ)/DZ(KKG_OLD))*(K_TERM(3)+Y_TERM(3))*WBAR*SFAC
-         DRAG_MAX(1) = ACCEL_X/(-UBAR-SIGN(1._EB,UBAR)*TWO_EPSILON_EB)
-         DRAG_MAX(2) = ACCEL_Y/(-VBAR-SIGN(1._EB,VBAR)*TWO_EPSILON_EB)
-         DRAG_MAX(3) = ACCEL_Z/(-WBAR-SIGN(1._EB,WBAR)*TWO_EPSILON_EB)
    END SELECT
+   DRAG_MAX(1) = ACCEL_X/(-UBAR-SIGN(1._EB,UBAR)*TWO_EPSILON_EB)
+   DRAG_MAX(2) = ACCEL_Y/(-VBAR-SIGN(1._EB,VBAR)*TWO_EPSILON_EB)
+   DRAG_MAX(3) = ACCEL_Z/(-WBAR-SIGN(1._EB,WBAR)*TWO_EPSILON_EB)
    IF (TWO_D) THEN
       ACCEL_Y  = 0._EB
       DRAG_MAX(2) = 0._EB