Changeset 3207 for palm

Timestamp:

Aug 27, 2018 12:55:33 PM (6 years ago)

Author:

schwenkel

Message:

Minor bugfix for sgs-velocities in case of cloud droplets

File:

• palm/trunk/SOURCE/lpm_advec.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_advec.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor bugfix for sgs-velocities in case of cloud droplets
+! 
+! 3189 2018-08-06 13:18:55Z raasch
 ! Bugfix: Index of the array dzw has to be k+1 during the interpolation. 
 ! Otherwise k=0 causes an abortion because dzw is allocated from 1 to nzt+1
@@ -977 +980 @@
                    IF ( use_sgs_for_particles )  THEN
                       lagr_timescale(n) = km(kp,jp,ip) / MAX( e(kp,jp,ip), 1.0E-20_wp )
-                      RL             = EXP( -1.0_wp * dt_3d / lagr_timescale(n) )
+                      RL             = EXP( -1.0_wp * dt_3d / MAX( lagr_timescale(n), &
+                                             1.0E-20_wp ) )
                       sigma          = SQRT( e(kp,jp,ip) )
 
@@ -1049 +1053 @@
                 IF ( use_sgs_for_particles )  THEN
                     lagr_timescale(n) = km(kp,jp,ip) / MAX( e(kp,jp,ip), 1.0E-20_wp )
-                    RL             = EXP( -1.0_wp * dt_3d / lagr_timescale(n) )
+                     RL             = EXP( -1.0_wp * dt_3d / MAX( lagr_timescale(n), &
+                                             1.0E-20_wp ) )
                     sigma          = SQRT( e(kp,jp,ip) )