Changeset 828 for palm/trunk/SOURCE/modules.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/modules.f90 (modified) (5 diffs)

palm/trunk/SOURCE/modules.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! +dissipation_classes, radius_classes, use_kernel_tables,
+! particle feature color renamed class
 !
 ! Former revisions:
@@ -1195 +1196 @@
     INTEGER ::  mpi_particle_type
 #endif
-    INTEGER ::  ibc_par_lr, ibc_par_ns, ibc_par_b, ibc_par_t,                  &
-                iran_part = -1234567, maximum_number_of_particles = 1000,      &
+    INTEGER ::  dissipation_classes = 10, ibc_par_lr, ibc_par_ns, ibc_par_b,   &
+                ibc_par_t, iran_part = -1234567,                               &
+                maximum_number_of_particles = 1000,                            &
                 maximum_number_of_tailpoints = 100,                            &
                 maximum_number_of_tails = 0,                                   &
@@ -1203 +1205 @@
                 number_of_initial_tails = 0, offset_ocean_nzt = 0,             &
                 offset_ocean_nzt_m1 = 0, particles_per_point = 1,              &
-                particle_file_count = 0, skip_particles_for_tail = 100,        &
-                sort_count = 0, total_number_of_particles,                     &
-                total_number_of_tails = 0
+                particle_file_count = 0, radius_classes = 20,                  &
+                skip_particles_for_tail = 100, sort_count = 0,                 &
+                total_number_of_particles, total_number_of_tails = 0
 
     INTEGER, PARAMETER ::  max_number_of_particle_groups = 10
@@ -1215 +1217 @@
                 particle_advection = .FALSE., random_start_position = .FALSE., &
                 read_particles_from_restartfile = .TRUE.,                      &
-                uniform_particles = .TRUE., use_particle_tails = .FALSE.,      &
-                use_sgs_for_particles = .FALSE.,                               &
+                uniform_particles = .TRUE., use_kernel_tables = .FALSE.,       &
+                use_particle_tails = .FALSE., use_sgs_for_particles = .FALSE., &
                 wang_kernel = .FALSE., write_particle_statistics = .FALSE.
 
@@ -1247 +1249 @@
                    origin_z, radius, rvar1, rvar2, rvar3, speed_x, speed_y, &
                    speed_z, weight_factor, x, y, z
-       INTEGER ::  color, group, tailpoints, tail_id
+       INTEGER ::  class, group, tailpoints, tail_id
     END TYPE particle_type