diff --git a/Source/part.f90 b/Source/part.f90 index 5b1dbd156c4..a4273cd71c1 100644 --- a/Source/part.f90 +++ b/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