Changeset 1671 for palm

Timestamp:

Sep 25, 2015 3:29:37 AM (9 years ago)

Author:

raasch

Message:

bugfix: ghostpoint exchange for array diss in case that sgs velocities are used
for particles

File:

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

palm/trunk/SOURCE/time_integration.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! bugfix: ghostpoint exchange for array diss in case that sgs velocities are used
+! for particles
 ! 
 ! Former revisions:
@@ -214 +215 @@
 
     USE particle_attributes,                                                   &
-        ONLY:  particle_advection, particle_advection_start, wang_kernel
+        ONLY:  particle_advection, particle_advection_start,                   &
+               use_sgs_for_particles, wang_kernel
 
     USE pegrid
@@ -446 +448 @@
                    CALL exchange_horiz( ql_vp, nbgp )
                 ENDIF
-                IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
+                IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
+                THEN
+                   CALL exchange_horiz( diss, nbgp )
+                ENDIF
 
                 IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs
@@ -510 +515 @@
                    CALL exchange_horiz( ql_vp, nbgp )
                 ENDIF
-                IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
+                IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
+                THEN
+                   CALL exchange_horiz( diss, nbgp )
+                ENDIF
 
                 IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs
@@ -597 +605 @@
                 CALL exchange_horiz( ql_vp, nbgp )
              ENDIF
-             IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
+             IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
+             THEN
+                CALL exchange_horiz( diss, nbgp )
+             ENDIF
 
              IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs