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

Timestamp:

Jun 8, 2017 2:59:01 PM (7 years ago)

Author:

schwenkel

Message:

Implemented splitting and merging algorithm

File:

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

palm/trunk/SOURCE/lpm.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implemented splitting and merging algorithm
+! 
+! 2233 2017-05-30 18:08:54Z suehring
 !
 ! 2232 2017-05-30 17:47:52Z suehring
@@ -142 +145 @@
         ONLY:  collision_kernel, deleted_particles, deallocate_memory,         &
                dt_write_particle_data, dt_prel, end_time_prel,                 &
-               grid_particles, number_of_particles, number_of_particle_groups, &
-               particles, particle_groups, prt_count, step_dealloc,            &
-               time_prel, time_write_particle_data, trlp_count_sum,            &
-               trlp_count_recv_sum, trnp_count_sum,                            &
-               trnp_count_recv_sum, trrp_count_sum, trrp_count_recv_sum,       &
-               trsp_count_sum, trsp_count_recv_sum,                            &
+               grid_particles, merging,  number_of_particles,                  &
+               number_of_particle_groups, particles, particle_groups,          &
+               prt_count, splitting, step_dealloc, time_prel,                  &
+               time_write_particle_data, trlp_count_sum, trlp_count_recv_sum,  &
+               trnp_count_sum, trnp_count_recv_sum, trrp_count_sum,            & 
+               trrp_count_recv_sum, trsp_count_sum, trsp_count_recv_sum,       &
                use_sgs_for_particles, write_particle_statistics
-
+                              
     USE pegrid
 
@@ -369 +372 @@
        IF ( dt_3d_reached )  time_prel = time_prel + dt_3d
 !
+!--    Apply splitting and merging algorithm
+       IF ( cloud_droplets )  THEN
+          IF ( splitting ) THEN
+             CALL lpm_splitting
+          ENDIF
+          IF ( merging ) THEN
+             CALL lpm_merging
+          ENDIF
+       ENDIF
+!
 !--    Move Particles local to PE to a different grid cell
        CALL lpm_move_particle