Changeset 4616 for palm

Timestamp:

Jul 21, 2020 10:09:46 AM (4 years ago)

Author:

schwenkel

Message:

Bugfix in case of strechting: k-calculation limited lower bound of 1

File:

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

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Bugfix in case of strechting: k-calculation limited lower bound of 1
+! 
+! 4589 2020-07-06 12:34:09Z suehring
 ! remove unused variables
 ! 
@@ -1418 +1421 @@
                                IF ( dz_stretch_level /= -9999999.9_wp  .OR.           &
                                     dz_stretch_level_start(1) /= -9999999.9_wp )  THEN
-                                  kp = MINLOC( ABS( tmp_particle%z - zu ), DIM = 1 ) - 1
+                                  kp = MAX( MINLOC( ABS( tmp_particle%z - zu ), DIM = 1 ) - 1, 1 )
                                ELSE
                                   kp = INT( tmp_particle%z / dz(1) + 1 + offset_ocean_nzt )
@@ -7833 +7836 @@
        IF ( dz_stretch_level /= -9999999.9_wp  .OR.         &
             dz_stretch_level_start(1) /= -9999999.9_wp )  THEN
-          kp = MINLOC( ABS( particle_array(n)%z - zu ), DIM = 1 ) - 1
+          kp = MAX( MINLOC( ABS( particle_array(n)%z - zu ), DIM = 1 ) - 1, 1 )
        ELSE
           kp = INT( particle_array(n)%z / dz(1) + 1 + offset_ocean_nzt )