Changeset 2417 for palm/trunk/SOURCE/user_lpm_advec.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/user_lpm_advec.f90 (modified) (3 diffs)

palm/trunk/SOURCE/user_lpm_advec.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Particle loops adapted for sub-box structure, i.e. for each sub-box the 
+! particle loop runs from start_index up to end_index instead from 1 to 
+! number_of_particles.
+! 
+! 2101 2017-01-05 16:42:31Z suehring
 !
 ! 2000 2016-08-20 18:09:15Z knoop
@@ -70 +75 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  ip   !<
-    INTEGER(iwp) ::  jp   !< 
-    INTEGER(iwp) ::  kp   !< 
-    INTEGER(iwp) ::  n    !< 
+    INTEGER(iwp) ::  ip   !< index of particle grid box, x-direction
+    INTEGER(iwp) ::  jp   !< index of particle grid box, y-direction
+    INTEGER(iwp) ::  kp   !< index of particle grid box, z-direction
+    INTEGER(iwp) ::  n    !< particle index
+    INTEGER(iwp) ::  nb   !< index of sub-box particles are sorted in
+
+    INTEGER(iwp), DIMENSION(0:7)  ::  start_index !< start particle index for current sub-box
+    INTEGER(iwp), DIMENSION(0:7)  ::  end_index   !< start particle index for current sub-box
 
 !
@@ -79 +88 @@
 !   number_of_particles = prt_count(kp,jp,ip)
 !   particles => grid_particles(kp,jp,ip)%particles(1:number_of_particles)
+!
+!   start_index = grid_particles(kp,jp,ip)%start_index
+!   end_index   = grid_particles(kp,jp,ip)%end_index
+!
 !   IF ( number_of_particles <= 0 )  CYCLE
 !   DO  n = 1, number_of_particles
-!
+!   DO  nb = 0, 7
+!      DO  n = start_index(nb), end_index(nb)
+!         particles(n)%xxx = 
+!      ENDDO
 !   ENDDO