Changeset 1519 for palm

Timestamp:

Jan 8, 2015 10:20:42 AM (10 years ago)

Author:

hoffmann

Message:

bugfix: check if radius ratio < 1 for computation of collisions

File:

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

palm/trunk/SOURCE/lpm_collision_kernels.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: Using the new particle structure, particles are not sorted by size.
+! Hence, computation of collision efficiencies must ensure that the ratio of 
+! two colliding droplets is < 1.
 ! 
 ! Former revisions:
@@ -846 +848 @@
 
              ir = ira(j)
-             rq = radclass(i) / radclass(j)
+             rq = MIN( radclass(i) / radclass(j), radclass(j) / radclass(i) )
              iq = INT( rq * 20 ) + 1
              iq = MAX( iq , 2)
@@ -1003 +1005 @@
 
              ir = ira(j)
-             rq = radclass(i) / radclass(j)
+             rq = MIN( radclass(i) / radclass(j), radclass(j) / radclass(i) )
 
              DO  kk = 2, 11