Changeset 4054 for palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

Timestamp:

Jun 27, 2019 7:42:18 AM (5 years ago)

Author:

raasch

Message:

bugfix for calculating the minimum particle time step

File:

• palm/trunk/SOURCE/lagrangian_particle_model_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for calculating the minimum particle time step
+! 
+! 4044 2019-06-19 12:28:27Z schwenkel
 ! Bugfix in case of grid strecting: corrected calculation of k-Index
 !
@@ -3738 +3741 @@
 !--          The particle timestep should not be too small in order to prevent
 !--          the number of particle timesteps of getting too large
-             IF ( dt_particle(n) < dt_min_part  .AND.  dt_min_part < dt_gap(n) )  THEN
-                dt_particle(n) = dt_min_part
+             IF ( dt_particle(n) < dt_min_part )  THEN
+                IF ( dt_min_part < dt_gap(n) )  THEN
+                   dt_particle(n) = dt_min_part
+                ELSE
+                   dt_particle(n) = dt_gap(n)
+                ENDIF
              ENDIF
              rvar1_temp(n) = particles(n)%rvar1