Changeset 4114 for palm/trunk/SOURCE

Timestamp:

Jul 23, 2019 2:09:27 PM (5 years ago)

Author:

schwenkel

Message:

Bugfix: Added working precision

File:

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

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: Added working precision for if statement
+! 
+! 4054 2019-06-27 07:42:18Z raasch
 ! bugfix for calculating the minimum particle time step
 ! 
@@ -1338 +1341 @@
 !
 !--                            In case of stretching the actual k index is found iteratively
-                               IF ( dz_stretch_level .NE. -9999999.9  .OR.           &
-                                    dz_stretch_level_start(1) .NE. -9999999.9 ) THEN
+                               IF ( dz_stretch_level .NE. -9999999.9_wp  .OR.           &
+                                    dz_stretch_level_start(1) .NE. -9999999.9_wp ) THEN
                                   kp = MINLOC( ABS( tmp_particle%z - zu ), DIM = 1 ) - 1
                                ELSE
@@ -7352 +7355 @@
 !
 !--    In case of stretching the actual k index must be found
-       IF ( dz_stretch_level .NE. -9999999.9  .OR.         &
-            dz_stretch_level_start(1) .NE. -9999999.9 ) THEN
+       IF ( dz_stretch_level .NE. -9999999.9_wp  .OR.         &
+            dz_stretch_level_start(1) .NE. -9999999.9_wp ) THEN
           kp = MINLOC( ABS( particle_array(n)%z - zu ), DIM = 1 ) - 1
        ELSE
@@ -7523 +7526 @@
                 IF( zw(k)   < particles_before_move(n)%z ) k = k + 1
                 IF( zw(k-1) > particles_before_move(n)%z ) k = k - 1 
-                
+
 !--             For lpm_exchange_horiz to work properly particles need to be moved to the outermost gridboxes
 !--             of the respective processor. If the particle index is inside the processor the following lines
@@ -7531 +7534 @@
                 j = MIN ( j , nyn )
                 j = MAX ( j , nys )
-                
+
                 k = MIN ( k , nzt )
                 k = MAX ( k , nzb+1 )
-                
+
 !
 !--             Check, if particle has moved to another grid cell.