Changeset 2245 for palm/trunk/SOURCE/lpm_exchange_horiz.f90

Timestamp:

Jun 2, 2017 2:37:10 PM (7 years ago)

Author:

schwenkel

Message:

Bugfix in Add_particles_to_gridcell

File:

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

palm/trunk/SOURCE/lpm_exchange_horiz.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
-! Bugfix in lpm_move_particle that prevents
-! false production of additional particles
+! 
 ! 
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Bugfix in Add_particles_to_gridcell:
+! Apply boundary conditions also in y-direction
+! 
+! 2151 2017-02-15 10:42:16Z schwenkel
+! Bugfix in lpm_move_particle that prevents
+! false production of additional particles
 !
 ! 2000 2016-08-20 18:09:15Z knoop
@@ -880 +884 @@
 
              IF ( jp == -1 )  THEN
-                move_also_south(nr_move_south)%y =                             &
-                   move_also_south(nr_move_south)%y + ( ny + 1 ) * dy
-                move_also_south(nr_move_south)%origin_y =                      &
-                   move_also_south(nr_move_south)%origin_y + ( ny + 1 ) * dy
+!
+!--             Apply boundary condition along y
+                IF ( ibc_par_ns == 0 )  THEN
+                   move_also_south(nr_move_south)%y =                          &
+                      move_also_south(nr_move_south)%y + ( ny + 1 ) * dy
+                   move_also_south(nr_move_south)%origin_y =                   &
+                      move_also_south(nr_move_south)%origin_y + ( ny + 1 ) * dy
+                ELSEIF ( ibc_par_ns == 1 )  THEN
+!
+!--                Particle absorption
+                   move_also_south(nr_move_south)%particle_mask = .FALSE.
+                   deleted_particles = deleted_particles + 1
+
+                ELSEIF ( ibc_par_ns == 2 )  THEN
+!
+!--                Particle reflection
+                   move_also_south(nr_move_south)%y       =                    &
+                      -move_also_south(nr_move_south)%y
+                   move_also_south(nr_move_south)%speed_y =                    &
+                      -move_also_south(nr_move_south)%speed_y
+
+                ENDIF
              ENDIF
           ELSEIF ( jp >= nyn+1 )  THEN
@@ -906 +928 @@
              move_also_north(nr_move_north) = particle_array(n)
              IF ( jp == ny+1 )  THEN
-                move_also_north(nr_move_north)%y =                             &
-                   move_also_north(nr_move_north)%y - ( ny + 1 ) * dy
-                move_also_north(nr_move_north)%origin_y =                      &
-                   move_also_north(nr_move_north)%origin_y - ( ny + 1 ) * dy
+!
+!--             Apply boundary condition along y
+                IF ( ibc_par_ns == 0 )  THEN
+
+                   move_also_north(nr_move_north)%y =                          &
+                      move_also_north(nr_move_north)%y - ( ny + 1 ) * dy
+                   move_also_north(nr_move_north)%origin_y =                   &
+                      move_also_north(nr_move_north)%origin_y - ( ny + 1 ) * dy
+                ELSEIF ( ibc_par_ns == 1 )  THEN
+!
+!--                Particle absorption
+                   move_also_north(nr_move_north)%particle_mask = .FALSE.
+                   deleted_particles = deleted_particles + 1
+
+                ELSEIF ( ibc_par_ns == 2 )  THEN
+!
+!--                Particle reflection
+                   move_also_north(nr_move_north)%y       =                    &
+                      -move_also_north(nr_move_north)%y
+                   move_also_north(nr_move_north)%speed_y =                    &
+                      -move_also_north(nr_move_north)%speed_y
+
+                ENDIF
              ENDIF
           ELSE