Changeset 835 for palm/trunk/SOURCE/lpm_collision_kernels.f90

Timestamp:

Feb 22, 2012 11:21:19 AM (12 years ago)

Author:

raasch

Message:

Bugfix: array diss can be used only in case of Wang kernel

File:

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

palm/trunk/SOURCE/lpm_collision_kernels.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: array diss can be used only in case of Wang kernel
 !
 ! Former revisions:
@@ -242 +242 @@
        radclass(1:radius_classes) = particles(pstart:pend)%radius * 100.0
 
-       epsilon = diss(k1,j1,i1) * 1.0E4   ! dissipation rate in cm**2/s**-3
+       IF ( wang_kernel )  THEN
+          epsilon = diss(k1,j1,i1) * 1.0E4   ! dissipation rate in cm**2/s**-3
+       ELSE
+          epsilon = 0.0
+       ENDIF
        urms    = 202.0 * ( epsilon / 400.0 )**( 0.33333333333 )