Changeset 828 for palm/trunk/SOURCE/init_particles.f90

Timestamp:

Feb 21, 2012 12:00:36 PM (12 years ago)

Author:

raasch

Message:

New:
---

Changed:

---

Optimization of collision kernels. Collision tables can be calculated once at
simulation start for defined radius (and dissipation) classes instead of
re-calculating them at every timestep and for the particle ensemble in
every gridbox.
For this purpose the particle feature color is renamed class.
New parpar parameters radius_classes and dissipation_classes.
(Makefile, advec_particles, check_parameters, data_output_dvrp, header, init_particles, lpm_collision_kernels, modules, package_parin, set_particle_attributes)

Lower limit for droplet radius changed from 1E-7 to 1E-8.
(advec_particles)

Complete re-formatting of collision code (including changes in names of
variables, modules and subroutines).
(advec_particles, lpm_collision_kernels)

Errors:

---

Bugfix: transformation factor for dissipation changed from 1E5 to 1E4
(advec_particles, lpm_collision_kernels)

File:

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

palm/trunk/SOURCE/init_particles.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! call of init_kernels, particle feature color renamed class
 !
 ! Former revisions:
@@ -72 +72 @@
     USE grid_variables
     USE indices
+    USE lpm_collision_kernels_mod
     USE particle_attributes
     USE pegrid
@@ -148 +149 @@
        IF ( pdz(j) == 9999999.9  .OR.  pdz(j) == 0.0 )  pdz(j) = pdz(j-1)
     ENDDO
+
+!
+!-- Initialize collision kernels
+    IF ( collision_kernel /= 'none' )  CALL init_kernels
 
 !
@@ -343 +348 @@
                             particles(n)%radius      = particle_groups(i)%radius
                             particles(n)%weight_factor =initial_weighting_factor
-                            particles(n)%color         = 1
+                            particles(n)%class         = 1
                             particles(n)%group         = i
                             particles(n)%tailpoints    = 0
@@ -470 +475 @@
              IF ( MINVAL( particles(1:n)%dvrp_psize  ) ==     &
                   MAXVAL( particles(1:n)%dvrp_psize  )  .AND. &
-                  MINVAL( particles(1:n)%color ) ==     &
-                  MAXVAL( particles(1:n)%color ) )  THEN
+                  MINVAL( particles(1:n)%class ) ==     &
+                  MAXVAL( particles(1:n)%class ) )  THEN
                 uniform_particles_l = .TRUE.
              ELSE
@@ -532 +537 @@
                 particle_tail_coordinates(1,2,nn) = particles(n)%y
                 particle_tail_coordinates(1,3,nn) = particles(n)%z
-                particle_tail_coordinates(1,4,nn) = particles(n)%color
+                particle_tail_coordinates(1,4,nn) = particles(n)%class
                 particles(n)%tailpoints = 1
                 IF ( minimum_tailpoint_distance /= 0.0 )  THEN
@@ -538 +543 @@
                    particle_tail_coordinates(2,2,nn) = particles(n)%y
                    particle_tail_coordinates(2,3,nn) = particles(n)%z
-                   particle_tail_coordinates(2,4,nn) = particles(n)%color
+                   particle_tail_coordinates(2,4,nn) = particles(n)%class
                    particle_tail_coordinates(1:2,5,nn) = 0.0
                    particles(n)%tailpoints = 2