Changeset 1929 for palm/trunk/SOURCE/lpm.f90

Timestamp:

Jun 9, 2016 4:25:25 PM (8 years ago)

Author:

suehring

Message:

several bugfixes in particle model and serial mode

File:

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

palm/trunk/SOURCE/lpm.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!
+! Call wall boundary conditions only if particles are in the vertical range of 
+! topography.
 ! 
 ! Former revisions:
@@ -111 +112 @@
 
     USE indices,                                                               &
-        ONLY: nxl, nxr, nys, nyn, nzb, nzt
+        ONLY: nxl, nxr, nys, nyn, nzb, nzb_max, nzt, nzb_w_inner
 
     USE kinds
@@ -300 +301 @@
                 CALL lpm_advec(i,j,k)
 !
-!--             Particle reflection from walls
-                IF ( topography /= 'flat' )  THEN
+!--             Particle reflection from walls. Only applied if the particles 
+!--             are in the vertical range of the topography. (Here, some
+!--             optimization is still possible.)
+                IF ( topography /= 'flat' .AND. k < nzb_max + 2 )  THEN
                    CALL lpm_boundary_conds( 'walls' )
                 ENDIF
@@ -334 +337 @@
        steps = steps + 1
        dt_3d_reached_l = ALL(grid_particles(:,:,:)%time_loop_done)
-
 !
 !--    Find out, if all particles on every PE have completed the LES timestep
@@ -359 +361 @@
 !--    Horizontal boundary conditions including exchange between subdmains
        CALL lpm_exchange_horiz
-
 !
 !--    Pack particles (eliminate those marked for deletion),