Changeset 1880 for palm

Timestamp:

Apr 20, 2016 9:36:50 AM (9 years ago)

Author:

hoffmann

Message:

Bugfix in lpm_collision_kernels

File:

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

palm/trunk/SOURCE/lpm_collision_kernels.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: The index of the larger particle has to be chosen for interpolation.
 ! 
 ! Former revisions:
@@ -27 +27 @@
 ! 1873 2016-04-18 14:50:06Z maronga
 ! Module renamed (removed _mod)
-!
 ! 
 ! 1858 2016-04-13 13:12:11Z hoffmann
@@ -786 +785 @@
           DO  i = 1, j
 
-             ir = ira(j)
+             ir = MAX( ira(i), ira(j) )
              rq = MIN( radclass(i) / radclass(j), radclass(j) / radclass(i) )
              iq = INT( rq * 20 ) + 1
@@ -945 +944 @@
           DO  i = 1, j
 
-             ir = ira(j)
+             ir = MAX( ira(i), ira(j) )
              rq = MIN( radclass(i) / radclass(j), radclass(j) / radclass(i) )