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

Timestamp:

Sep 6, 2017 3:22:27 PM (7 years ago)

Author:

suehring

Message:

Major bugfix in modeling SGS particle speeds.

File:

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

palm/trunk/SOURCE/lpm.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Major bugfixes in modeling SGS particle speeds (since revision 1359).
+! Particle sorting added to distinguish between already completed and 
+! non-completed particles.
+! 
+! 2263 2017-06-08 14:59:01Z schwenkel
 ! Implemented splitting and merging algorithm
 ! 
@@ -140 +145 @@
 
     USE lpm_pack_arrays_mod,                                                   &
-        ONLY:  lpm_pack_all_arrays
+        ONLY:  lpm_pack_all_arrays, lpm_sort
 
     USE particle_attributes,                                                   &
@@ -259 +264 @@
        CALL cpu_log( log_point_s(44), 'lpm_advec', 'pause' )
        
-       grid_particles(:,:,:)%time_loop_done = .TRUE.
 !
 !--    If particle advection includes SGS velocity components, calculate the
@@ -265 +269 @@
 !--    horizontally averaged profiles of the SGS TKE and the resolved-scale 
 !--    velocity variances)
-
        IF ( use_sgs_for_particles  .AND.  .NOT. cloud_droplets )  THEN
           CALL lpm_init_sgs_tke
        ENDIF
+
+!
+!--    In case SGS-particle speed is considered, particles may carry out 
+!--    several particle timesteps. In order to prevent unnecessary 
+!--    treatment of particles that already reached the final time level, 
+!--    particles are sorted into contiguous blocks of finished and 
+!--    not-finished particles, in addition to their already sorting
+!--    according to their sub-boxes. 
+       IF ( .NOT. first_loop_stride  .AND.  use_sgs_for_particles )            &
+          CALL lpm_sort
 
        DO  i = nxl, nxr
@@ -323 +336 @@
 !--             optimization is still possible.)
                 IF ( topography /= 'flat' .AND. k < nzb_max + 2 )  THEN
-                   CALL lpm_boundary_conds( 'walls' )
+                   CALL lpm_boundary_conds( 'walls' )!, i, j, k )
                 ENDIF
 !
@@ -333 +346 @@
 !--             the top or bottom boundary and delete those particles, which are
 !--             older than allowed
-                CALL lpm_boundary_conds( 'bottom/top' )
+                CALL lpm_boundary_conds( 'bottom/top' ) !, i, j, k )
 !
 !---            If not all particles of the actual grid cell have reached the 
-!--             LES timestep, this cell has to to another loop iteration. Due to
-!--             the fact that particles can move into neighboring grid cell, 
-!--             these neighbor cells also have to perform another loop iteration
+!--             LES timestep, this cell has to do another loop iteration. Due to
+!--             the fact that particles can move into neighboring grid cells, 
+!--             these neighbor cells also have to perform another loop iteration.
+!--             Please note, this realization does not work properly if 
+!--             particles move into another subdomain. 
                 IF ( .NOT. dt_3d_reached_l )  THEN
-                   ks = MAX(nzb+1,k)
-                   ke = MIN(nzt,k)
-                   js = MAX(nys,j)
-                   je = MIN(nyn,j)
-                   is = MAX(nxl,i)
-                   ie = MIN(nxr,i)
+                   ks = MAX(nzb+1,k-1)
+                   ke = MIN(nzt,k+1)
+                   js = MAX(nys,j-1)
+                   je = MIN(nyn,j+1)
+                   is = MAX(nxl,i-1)
+                   ie = MIN(nxr,i+1)
                    grid_particles(ks:ke,js:je,is:ie)%time_loop_done = .FALSE.
+                ELSE
+                   grid_particles(k,j,i)%time_loop_done = .TRUE.
                 ENDIF
 
@@ -392 +409 @@
 !--    determine new number of particles
        CALL lpm_pack_all_arrays
-
 !
 !--    Initialize variables for the next (sub-) timestep, i.e., for marking