Changeset 792 for palm/trunk/SOURCE/advec_particles.f90

Timestamp:

Dec 1, 2011 12:23:23 AM (12 years ago)

Author:

raasch

Message:

further adjustments for speedup of particle code

File:

• palm/trunk/SOURCE/advec_particles.f90 (modified) (4 diffs)

palm/trunk/SOURCE/advec_particles.f90

@@ -4 +4 @@
 ! Current revisions:
 ! ------------------
-! 
+! particle arrays (particles, particles_temp) implemented as pointers in
+! order to speed up sorting (see routine sort_particles)
 !
 ! Former revisions:
@@ -3995 +3996 @@
     INTEGER ::  i, ilow, j, k, n
 
-    TYPE(particle_type), DIMENSION(1:number_of_particles) ::  particles_temp
+    TYPE(particle_type), DIMENSION(:), POINTER ::  particles_temp
 
 
@@ -4001 +4002 @@
 
 !
-!-- Initialize the array used for counting and indexing the particles
-    prt_count = 0
+!-- Initialize counters and set pointer of the temporary array into which
+!-- particles are sorted to free memory
+    prt_count  = 0
+    sort_count = sort_count +1
+
+    SELECT CASE ( MOD( sort_count, 2 ) )
+
+       CASE ( 0 )
+
+          particles_temp => part_1
+          part_1 = particle_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+                                  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+                                  0.0, 0, 0, 0, 0 )
+
+       CASE ( 1 )
+
+          particles_temp => part_2
+          part_2 = particle_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+                                  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+                                  0.0, 0, 0, 0, 0 )
+
+    END SELECT
 
 !
@@ -4055 +4076 @@
     ENDDO
 
-    particles(1:number_of_particles) = particles_temp
+!
+!-- Redirect the particles pointer to the sorted array
+    SELECT CASE ( MOD( sort_count, 2 ) )
+
+       CASE ( 0 )
+
+          particles => part_1
+
+       CASE ( 1 )
+
+          particles => part_2
+
+    END SELECT
 
 !