Changeset 3524 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Nov 14, 2018 1:36:44 PM (5 years ago)

Author:

raasch

Message:

unused variables removed, missing working precision added, missing preprocessor directives added, bugfix concerning allocation of t_surf_wall_v in nopointer case, declaration statements rearranged to avoid compile time errors, mpi_abort arguments replaced to avoid compile errors

File:

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

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added missing working precision
+! 
+! 3361 2018-10-16 20:39:37Z knoop
 ! ocean renamed ocean_mode
 ! 
@@ -804 +807 @@
 !
 !--                            Determine the grid indices of the particle position
-                               ip = tmp_particle%x * ddx
-                               jp = tmp_particle%y * ddy
-                               kp = tmp_particle%z / dz(1) + 1 + offset_ocean_nzt                               
+                               ip = INT( tmp_particle%x * ddx )
+                               jp = INT( tmp_particle%y * ddy )
+                               kp = INT( tmp_particle%z / dz(1) + 1 + offset_ocean_nzt )
                                DO WHILE( zw(kp) < tmp_particle%z ) 
                                   kp = kp + 1
@@ -1184 +1187 @@
                 IF ( particles(n)%weight_factor - FLOOR(particles(n)%weight_factor,KIND=wp) &
                      .GT. random_function( iran_part ) )  THEN
-                   particles(n)%weight_factor = FLOOR(particles(n)%weight_factor,KIND=wp) + 1.0
+                   particles(n)%weight_factor = FLOOR(particles(n)%weight_factor,KIND=wp) + 1.0_wp
                 ELSE
                    particles(n)%weight_factor = FLOOR(particles(n)%weight_factor,KIND=wp)