Changeset 790 for palm/trunk/SOURCE/diffusion_e.f90

Timestamp:

Nov 29, 2011 3:11:20 AM (13 years ago)

Author:

raasch

Message:

Bugfix for output of mean particle radius + preliminary works for implementing the Wang collision kernel

File:

• palm/trunk/SOURCE/diffusion_e.f90 (modified) (5 diffs)

palm/trunk/SOURCE/diffusion_e.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! diss is also calculated in case that the Wang kernel is used
 !
 ! Former revisions:
@@ -160 +161 @@
 !--          Store dissipation if needed for calculating the sgs particle
 !--          velocities
-             IF ( use_sgs_for_particles )  THEN
+             IF ( use_sgs_for_particles  .OR.  wang_collision_kernel )  THEN
                 DO  j = nys, nyn
                    DO  k = nzb_s_inner(j,i)+1, nzt
@@ -250 +251 @@
 !--          Store dissipation if needed for calculating the sgs particle
 !--          velocities
-             IF ( use_sgs_for_particles )  THEN
+             IF ( use_sgs_for_particles  .OR.  wang_collision_kernel )  THEN
                 DO  j = nys, nyn
                    DO  k = nzb_s_inner(j,i)+1, nzt
@@ -264 +265 @@
 !
 !--    Boundary condition for dissipation
-       IF ( use_sgs_for_particles )  THEN
+       IF ( use_sgs_for_particles  .OR.  wang_collision_kernel )  THEN
           DO  i = nxl, nxr
              DO  j = nys, nyn
@@ -368 +369 @@
 !
 !--    Store dissipation if needed for calculating the sgs particle velocities
-       IF ( use_sgs_for_particles )  THEN
+       IF ( use_sgs_for_particles  .OR.  wang_collision_kernel )  THEN
           DO  k = nzb_s_inner(j,i)+1, nzt
              diss(k,j,i) = dissipation(k)